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Abstract 

The jet parton correlation has been studied for different cut scenarios and various 
jet algorithms. These migrations are getting important for measuring the strong 
couphng constant and the NLO gluon density g{^, Q^) from jet rates in DIS at 
HERA. 
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Kapitel 1 
Einleitung 



Ein wesentlicher Charakterzug der Menschen ist die Neugierde. Immer schon haben 
wir uns nicht mit der bloßen Existenz eines Phänomens zufriedengegeben, sondern 
haben versucht die Ursache dafür zu ergründen. 

Die älteste Frage ist die nach dem Woher der Welt und ihrer Zusammensetzung. 
Eine entscheidende Idee hatte hierbei schon der griechische Philosoph Demokrit, für 
den die Welt aus elementaren, unteilbaren Bausteinen ("atomos") aufgebaut war. 

Heute ist diese Frage immer noch nicht beantwortet, auch wenn die moderne 
Naturwissenschaft Schritt für Schritt zu immer kleineren Strukturen vordringt (siehe 
Abbildung Q. 

Eine Voraussetzung dafür war die Entwicklung des Modellbegriffes. Unter einem 
Modell verstehen wir in der Physik eine Theorie, die einen bestimmten Teil der Natur 
beschreibt, ohne daß Messungen den Vorhersagen des Modells widersprechen. Aus 
den meisten Theorien können wir zusätzliche Aussagen über bisher nicht beobachtbare 
Bereiche und neue Phänomene treffen. Die Aufgabe des Experimentalphysikers ist es, 
dieses zu überprüfen. Stimmen die Messungen nicht mit der Theorie überein, so sind 
die bestehenden Modelle in diesem Bereich nicht mehr gültig und die Theoretiker 
sind dazu aufgerufen, diese zu erweitern oder neue Modelle zu entwickeln. Dieses 
gegenseitige Wechselspiel ist sicher ein wesentlicher Grund für den Erfolg der heutigen 
Naturwissenschaft . 

Auch wenn die heute so häufig beschrieene Anwendbarkeit der Forschungsergebnisse 
bei der Grundlagenforschung nicht direkt erkennbar ist, so befriedigt sie dennoch die 
Neugier der Menschen und leistet einen Beitrag zum technischen und philosophischen 
Fortschritt. Als Beispiele seien hier nur die Relativitätstheorie oder das Heisenbergsche 
Unschärfeprinzip erwähnt, die das Weltbild und die Philosophie dieses Jahrhunderts 
wie kaum eine andere Entwicklung geprägt haben. 

Das zur Zeit gültige Standardmodell der Elementarteilchenphysik ist in der Lage 
alle heute bekannten Effekte der Hochenergiephysik durch drei Quark- und drei 
Lepton- Familien zu beschreiben, deren Mitglieder über Austauschteilchen miteinander 



wechselwirken. Die Partikel sind in Tabelle 1.1 aufgeführt. Es wird zwischen links- 



(L) und rechtshändigen (R) Teilchen unterschieden, da Neutrinos nur linkshändig in 



3 



4 



KAPITEL 1. EINLEITUNG 




+0,01 m 
Kristall 

1/10.000.000 

10"^m 
Molekül 

1/10 

10-^°m 
Atom 



1/10.000 

10-^'m 
Atomkern 



1/10 

10'^^m 
Proton 



1/1.000 

<10"^^m 
Elektron, 
Quark 



>- 

CO 

< 

X 



CO 
O 



O) 

c 



(0 
(0 

c 
o 



u 
c 



Abbildung 1.1: Vergleich der Größenordnungen verschiedener Strukturen 

(gedruckt mit freundlicher Genehmigung der DESY PR Abteilung) 
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Leptonen 
(Spin i) 



Quarks 
(Spin i) 




UR, du Cr, Sr tR, Ör 



Austauschteilchen (Spin 1) 


elektromagnetisch 


7 


schwach 




stark 


5fi,z = 1, . . . ,8 


Skalare (Spin 0) 


Higgsfeld 


h 



Tabelle 1.1: Elementarteilchen im Standardmodell 



der Natur vorkommen. Die Quarks sind alle bereits experimentell beobachtet worden, 
wobei die Existenz des Top-Quarks t nach ersten Hinweisen 1994 nun experimentell 
bestätigt wurde ||CDF95|| , ||DZe95|| . Bei den Lept onen ist nur das Tau-Neutrino Uj- 
bisher unbeobachtet. 



Die ebenfalls in Tabelle |LJ gezeigten Austauschteilchen übertragen die Wechselwirkung 
zwischen den Teilchen. Diese Partikel sind alle experimentell bestätigt. Ein wichtiges 
Ziel in der Hochenergiephysik ist, die verschiedenen Wechselwirkungen durch den 
gleichen Mechanismus mit einer einzigen Kopplungskonstante zu beschreiben. Diese 
sogenannte Vereinheitlichung ist bei der elektromagnetischen und der schwachen 
Kraft bereits gelungen. Um dies durch die Modelle der Großen Vereinheitlichung 
(GUT, grand unified theory) ebenfalls für die starke Kraft zu erreichen, ist eine 
genaue Kenntnis der starken Kopplungskonstante wesentlich. Die Vereinigung der 
Kopplungskonstanten wird bei einer heute nicht erreichbaren Energie erwartet. 

Ein Problem des Standardmodells ist die Masse der Teilchen zu erklären. Der 
Higgs-Mechanismus erreicht dies durch ein zusätzliches Teilchen, das Higgs-Boson 
h. Es besteht die Hoffnung, daß die Existenz dieses Teilchens bei den geplanten 
Beschleunigern LEP 200 oder LHC nachgewiesen werden kann. 

Das Standardmodell ist sicher nicht die letzte Theorie. Neben der großen vereinheitlichten 
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Theorie ist die hohe Anzahl an Eingabeparametern unbefriedigend. Ein Zusammenhang 
zwischen diesen Parametern wird in neueren Modellen wie der Supersymmetrie und 
der String-Theorie gesucht. Gleichzeitig stellen diese neue Vorhersagen auf, so daß 
sich in den nächsten Jahren die Möglichkeit weiterer Experimente und Messungen 
ergeben wird. 

Es bleiben somit noch viele Fragen offen, deren Beantwortung wir jedoch mit den 
heute bestehenden und zukünftigen Beschleunigern sowie durch engere Kooperation 
mit anderen Bereichen physikalischer Forschung (z.B. der Kosmologie) immer näherkommen 
werden. So werden wir hoffentlich irgendwann die Frage nach dem Woher und dem 
Woraus der Welt beantworten können. 



Kapitel 2 

Experimentelle Grundlagen 



2.1 Teilchenbeschleuniger 

In der Teilchenphysik werden Beschleuniger zur Bestimmung der Eigenschaften von 
Elementarteilchen eingesetzt. 

Zwei Teilchenreaktionen ermöglichen es, Daten zu gewinnen, Streuungen und 
Annihilationen . 

Zum einen kann man über Streuungen durch Beschuß des zu untersuchenden 
Teilchens dessen innere Struktur vermessen. Die erreichbare Auflösung Ax wird vom 
Viererimpulsübertrag ^/Q^ des Austauschteilchens über die Beziehung 

A. = ^ (2.1) 



bestimmt. Der Viererimpulsübertrag ist durch die zur Verfügung stehende Schwerpunktsenergie 

begrenzt. 

Zum anderen wird durch Annihilation von Teilchen und Antiteilchen eine definierte 
Schwerpunktsenergie erreicht, mit der neue Teilchen durch Resonanzen im Wirkungsquerschnitt 
nachgewiesen und untersucht werden können bzw. deren Masse und Lebensdauer 
bestimmt wird. Dazu muß jedoch die Stoßenergie in der Größenordnung der Masse 
des zu erzeugenden Teilchens sein. 

Beim Bau eines Beschleunigers, mit dem Teilchenreaktionen mit hoher Genauigkeit 
vermessen werden sollen, müssen verschiedene Designparameter beachtet werden. 
Von großer Bedeutung sind hier die Teilchensorte der Stoßpartner, die maximale 
Energie und die Luminosität. 

Eine möglichst hohe Schwerpunktsenergie erreicht man, indem zwei bewegte 
Partikel aneinander streuen bzw. annihilieren. Aus diesem Grunde sind sogenannte 
Fixed- Target-Experimente, d.h. Beschleuniger bei denen ein beschleunigtes Teilchen 
auf ein feststehendes, ruhendes Ziel geschossen wird, nur noch selten anzutreffen. Die 
Berechnung der Schwerpunktsenergie erfolgt unter Vernachlässigung der Teilchenmassen 
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bewegter Teilchen durchQ 



El 



Schwerpunkt 



{E^ + E^f - {pi + P2f 



(2.2) 
(2.3) 
(2.4) 



-^'Fixed-Target 




Die verbreitetste Bauart von Beschleunigern ist die eines Ringes, in dem die 
beiden Stoßpartner gegenläufig in einem Kreis umlaufen und nach und nach auf die 
Maximalenergie beschleunigt werden. Die Nachteile dieses Verfahrens liegen in der 
nötigen Synchronisation der Beschleunigungsvorgänge und dem Energieverlust durch 
Synchrotronstrahlung aufgrund der gekrümmten Bahn. Für die Zukunft sind jedoch 
auch einige Experimente geplant, bei denen die Stoßpartner durch zwei Linearbeschleuniger, 
die die Teilchen in einer geraden Röhre beschleunigen, aufeinander geschossen werden. 
Dazu muß jedoch die Stärke der beschleunigenden elektromagnetischen Felder vergrößert 
werden, da sonst die Beschleunigungsstrecken zu lang sind. Auch hier ist die Synchronisation 
der beiden Teilchenstrahlen kompliziert. Ein weiterer Nachteil ist, daß sich die 
Teilchenstrahlen nur einmal kreuzen und danach das Beschleunigersystem verlassen. 

Neben der Schwerpunktsenergie ist die spezifische Luminosität L einer der 
grundlegenden Parameter. Zusammen mit der Maschinenlaufzeit t läßt sich die 
integrierte Luminosität 



berechnen. Die zu erwartetende Ereigniszahl für einen Prozeß mit gegebenen 
Wirkungsquerschnitt a ergibt sich aus 



Eine hohe Luminosität erreicht man durch kleine Strahlquerschnitte und hohe 
Teilchenströme. 

Alle Colliderexperimente sind in eine der drei folgenden Kategorien nach Art 
der Stoßpartner einzuordnen. Die Vorteile bei der Beschleunigung von Elektronen 
entstehen aufgrund der geringen Masse. Schon bei einer Energie von 37 MeV erreichen 
die Elektronen eine Geschwindigkeit von 99.99% der Lichtgeschwindigkeit. Dadurch 
müssen bei der Beschleunigung im Ring nur kleine zeitliche Synchronisationen berücksichtigt 
werden. Die Probleme bei der Verwendung von Elektronen bestehen in den hohen 
Energieverlusten pro Umlauf durch Synchrotronstrahlung, die die maximale Strahlenergie 
begrenzen oder den Durchmesser des Ringes ansteigen lassen. Bei den schwereren 
Protonen hingegen ist ein Geschwindigkeitsbruchteil von 99.99% der Lichtgeschwindigkeit 
erst bei einer Teilchenenergie von 66.5 GeV erreicht, dadurch wird der Beschleunigungsvorgang 
komplizierter. Die erreichbaren Energien liegen höher, wobei berücksichtigt werden 
muß, daß bei tief- inelastischen Stößen nur ein Bruchteil xb des Protonimpules im 
harten Prozeß zur Verfügung steht. 

^Formel für CoUider gilt nur für Kollisionen von Teilchen aus entgegengesetzten Richtungen. 




N = C-a. 



2.1. TEILCHENBESCHLEUNIGER 
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maximale 


spezifische 


Name 


Ort 


Teilchensorte 


Strahlenergie 


Luminosität 








[GeV] 


[1030cm-2s^i] 


PETRAQ 


Hamburg 


e~^e~ 


23.4 


24 bei 17.5GeV 


PEPQ 


Stanford 


e~^e~ 


15 


60 


TRISTAN 


Tsukuba 


e~^e~ 


32 


37 


SLC 


Stanford 


e'^e~ 


50 


0.35 


LEP 


Genf 


e'^e~ 


55 


11 


LEP 200 


Genf 


e~^e~ 


700 


11 


SppS 


Genf 


ppb 


315 


6 


Tevatron I 


Chicago 


ppb 


1000 


10 


HERA 


Hamburg 


e^p 


30/820 


16 


LHcg 


CERN 


pp 


7000 


10000 



"Datennahme 1986 beendet 
''Datennahme 1990 beendet 

■^stufenweise Erhöhung der Energie von LEP auf eine Schwerpunktenergie von max. 193 GeV 
'^geplant, Fertigstellung voraussichtlich 2004 



Tabelle 2.1: Liste der Beschleuniger : Daten, außer LEP200, entnommen /{Far94\l 



Bei e~^e~ Beschleunigern werden Elektronen mit Positronen mit entgegengesetzt 
gleichem Impuls zur Kollision gebracht. Die dabei freiwerdende Energie kann zur 
Erzeugung von Teilchen verwendet werden, die die Quantenzahlen des Vakuums 
haben. Die Vorteile liegen in der einfachen Reaktionskinematik, da das Laborsystem 
gleichzeitig Schwerpunktsystem ist und bei Annihilation eine immer konstante, durch 
die Strahlenergie gegebene Schwerpunktsenergie vorliegt. 

Bei pp Maschinen streuen Protonen an deren Antiteilchen. Der Vorteil dieser 
Beschleuniger liegt in der höheren Energie der Strahlteilchen. Die Nachteile entstehen 
durch die innere Substruktur der Protonen, so daß nur ein Teil des Protonimpules 
in der Reaktion zur Verfügung steht und die Reste der Protonen zusätzliche Energie 
im Detektor hinterlassen. Dadurch entstehen hohe Jetmultiplizitäten. Ein weiteres 
Problem entsteht dadurch, daß die Strahlteilchen durch die starke Wechselwirkung 
im Anfangszustand zusätzliche Hintergrundereignisse hervorrufen können. Für pp 
CoUider, wie z.B. der inzwischen genehmigt LHC, gelten die gleichen Vor- und 
Nachteile, hinzu kommt jedoch, daß eine Annihilation der Protonen oder deren 
Konstituenten hier nicht möglich ist. 

Schließlich existiert noch die Möglichkeit Partikel unterschiedlicher Teilchensorte 
kollidieren zu lassen. Bei e^p Beschleunigern werden hochenergetische Protonen auf 
Elektronen oder Positronen geschossen. Die Vorteile liegen in der inneren Substruktur 
nur eines der Strahlteilchen, des Protons. Diese kann durch Beobachtung des gestreuten 
Elektrons und des hadronischen Endzustandes vermessen werden. Da der Energieübertrag 
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Abbildung 2.1: Der HERA Speicherring und die Vorbeschleuniger 



von Ereignis zu Ereignis unterschiedlich ist, können die Untersuchungen mit verschiedenen 
Schwerpunktsenergien ohne Veränderung der Strahlenergien durchgeführt werden. 
Die Nachteile entstehen durch die Verschiedenheit von Labor- und Schwerpunktsystem. 
In der Tabelle p.l| sind die wichtigsten Beschleuniger aufgeführt. 



2.2 HERA 

Der Beschleuniger HERA (Hadron Elektron Ring Anlage) ist der einzige e^p Beschleuniger 
und wurde nach siebenjähriger Bauzeit 1991 in Betrieb genommen. 

Die Beschleuniger, die auf dem DESY-Gelände ( Deutsches Elektronen Synchro- 
tron) bereits vorhanden waren, werden hierbei als Vorbeschleuniger benutzt. Eine 



Ubersicht über die einzelnen Geräte gibt die Abbildung |2?1 . 

Der Protonenring ist mit supraleitenden Magneten versehen. Bei der Injektion in 
HERA sind die Protonen bereits auf 40 GeV vorbeschleunigt. Die entgegengesetzt 
umlaufenden Elektronen^ haben eine Einschußenergie in HERA von 14 GeV. Dort 
werden herkömmliche Magnete verwendet, da die Energie und Masse der Elektronen 
wesentlich niedriger ist und so kleinere Magnetfelder benötigt werden. 



^In diesem Ring können sowohl Elektronen als auch deren Antiteilchen, die Positronen, 
beschleunigt werden. Nachdem bis Mitte 1994 Elektronen verwendet wurden, werden nun 
Positronen als Stol3partner der Protonen benutzt. Im folgenden wird trotzdem meist 
ausschließlich der Begriff Elektronen verwendet, da die Betrachtungen hier unabhängig von 
der Art der Strahllcptonen sind. Unterschiede ergeben sich aber z.B. in der Berechnung der 
Wirkungsquerschnitte elektroschwacher Prozesse. 



2.3. DER Hl DETEKTOR 



11 





Designparameter 




1994 




p-Ring 


e^-Rin 




p-Rin^ 


l e='=-Ring 


Maximalenergie [GeV] 


820 


30 




820 


27.55 


Anzahl Teilclienpakete 


210 


210 




153+170 153+15 


Stralilstrom/Paket [/lA] 


760 


290 




241 


101 


Einscliußenergie in HERA [GeV] 


40 


14 








magnetisches Ablenkfeld [T] 


4.65 


0.165 








Temperatur der Magnete 


4.4 K 


294 K 








Strahlungsverlust pro Umlauf 


125MeV 6eV 








Schwerpunktenergie [GeV] 




314 






300.6 


spez. Luminosität [10^''cm~^s~^] 




15 






1.40 


Umfang [m] 






6336 




Kollisionsrate [MHz] 






10.4 





"Die zweite Zahl gibt die Anzahl der Pilotbunche an, das sind die Strahlpakete, die nicht 
auf ein Paket des anderen Strahls treffen und so zur Messung des Strahl- Wand und Strahl-Gas 
Untergrundes dienen. 

''Der Wert von f994 entspricht der bei Hl gemessenen durchschnittlichen spezifischen 
Luminosität. Der höchste Wert war 4.3 • fO"^°cm~^s~^ 

Tabelle 2.2: Einige technische Parameter von HERA : Die Designwerte wurden 
hei der Datennahme 1994 teilweise noch nicht erreicht. 

Die Strahlteilchen kollidieren in den Wechselwirkungspunkten unter einem Winkel 
von 180°. Die Protonen und Elektronen werden zu sogenannten Bunches zusammengefaßt, 
d.h. zu Paketen aus Teilchen, die aufgrund der Beschleunigung durch elektrische 
Wechselfelder in Strahlrichtung eine begrenzte Ausdehnung haben. 

Es existieren vier Hallen, in denen Experimente aufgebaut werden können. Neben 
den beiden Vielzweckdetektoren Hl und ZEUS existiert seit 1995 mit HERMES ein 
drittes Experiment, das über ein Gastarget die Streuung an polarisierten Elektronen 
bzw. Positronen untersucht. Zur Zeit ist ein viertes Experiment HERA-B in Planung. 

Die technischen Daten von HERA sind in Tabelle ^ zusammengefaßt. 

2.3 Der Hl Detektor 

Der Hl Detektor ist einer der beiden Vielzweck 47r Detektoren am HERA Speicherring. 
Seine Aufgaben sind die Detektierung und Bestimmung von Reaktionsprodukten bei 
Stößen von Protonen und Elektronen. 

Dazu ist der nahezu komplette Raumbereich um den Wechselwirkungspunkt mit 
Spurkammern und Kalorimetern ausgestattet. Spurkammern dienen zur Bestimmung 
der Bahn geladener Teilchen. An der Krümmung dieser Bahnen in einem externen 
Magnetfeld läßt sich das Ladungsvorzeichen und in Verbindung mit der Stärke der 
Ionisation die Art des Teilchens bestimmen. Kalorimeter dienen zur Energiemessung, 
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indem sie die Menge der lonisationsladungen in den einzelnen Zellen bestimmen. 

Aus diesen Informationen kann man die Flugbahnen der Reaktionsprodukte und 
die Vertexposition, d.i. der Ort der Kollision, rekonstruieren. Außerdem sind noch 
Informationen über die Zeitpunkte der Ereignisse hilfreich. 

Da die Protonen eine höhere Strahlenergie haben und bei tief-inelastischen Streuungen 
nicht intakt bleiben, entspricht das Laborsystem nicht dem Schwerpunktsystem der 
Reaktion. Aus diesem Grunde sind die Anforderungen an die Meßgeräte in Proton- 
und in Elektronrichtung unterschiedlich und der Detektor selber nicht symmetrisch 
aufgebaut. 



Die Abbildung ^]2| zeigt den Hl Detektor in einer Rißzeichnung. Das Hl Koordinatensystem 
ist rechthändig und liegt so, daß die 2;-Achse entlang des Strahlrohres in Richtung 
der Protonen weist. Die x- Achse zeigt zum Speicherringmittelpunkt und die ?/- Achse 
somit nach oben. Der Nullpunkt liegt im nominalen Wechselwirkungspunkt. Meist 
erfolgt die Beschreibung jedoch in Polarkoordinaten, wobei 6 den Polarwinkel zur 
^-Achse darstellt. Die ungestreuten Protonen laufen unter dem Winkel 6 = 0° aus. 
Der Azimutal Winkel ip beschreibt die Lage in der x — y-Ebene. 

Die für diese Arbeit wichtigsten Komponenten werden hier nun näher erläutert, 
die Zahlen beziehen sich dabei auf Abbildung p.2| ||H1 93b|| : 



Spule (6) : Die supraleitende Spule erzeugt ein Magnetfeld, das eine Krümmung der 
Flugbahnen geladener Teilchen bewirkt. Dadurch kann das Ladungvorzeichen 
und der Impuls eines Teilchen bestimmt werden. Die Spule ist 5.16 m lang und 
hat einen Radius von 2.75 m. Bei den Spurkammern in der Mitte des Detektors 
erzeugt sie ein nahezu homogenes Magnetfeld von 1.2 T. 

Flüssig-Argon-Kalorimeter (LAr, 4 und 5) : Das mit flüssigem Argon gefüllte 
Kalorimeter deckt den Winkelbereich von vier bis 153 Grad ab. Durch seine 
Zweiteilung in einen hadronischen äußeren Teil und ein elektromagnetisches 
inneres Kalorimeter ist es in der Lage, nicht nur die Energie der Teilchen 
zu messen, sondern auf Grund des Verhältnisses der Ladungsdepositionen in 
den beiden Komponenten einen Hinweis auf die Art des Teilchens zu geben. 
Dies ist für die Berechnung der Energie wesentlich, da der Energieverlust pro 
Weglänge bei Hadronen und Leptonen unterschiedlich ist. Die charakterischen 
Größen zu den beiden Bereichen sind in Tabelle zusammengefaßt. Die 
Energieauflösung ist dabei über die Formel 



E/l GeV 




mit den Anteilen durch die Samplingfrequenz (A), das elektronische Rauschen 
(B) und die Kalibrationsunsicherheit (C) gegeben. 



Rückwärtiges elektromagnetisches Kalorimeter (BEMC, 12) : Das BEMC 
ist in Rückwärtsrichtung montiert und überdeckt den Winkelbereich von 151° 



3. DER Hl DETEKTOR 



HERA Experiment H1 




I 1 I Strahlrohr und Strahl mag nete 

[2~| Zentrale Spurkammern 

I 3 I Vorwärtsspurkammern und Übergangsstrahlungsmodul 

[4] Elektromagnetisches Kalorimeter (Blei) 

[s] Hadronisches Kalorimeter (Edelstahl) 

[g] Supraleitende Spule (1 .2T) 

Pt] Kompensationsmagnet 

[~8] Helium-Kälteanlage 



Flüssig-Argon 



[9] Myon-Kammern 

[TÖ| Instrumentiertes Eisen (Streamer-Röhren) 

[VT] Myon-Toroid-Magnet 

[12] warmes elektromagnetisches Kalorimeter 

[T3I Vorwärts-Kalorimeter 

[T4I Betonabschirmung 

[Tsl Flüssig-Argon-Kryostat 



Abbildung 2.2: Der Hl Detektor 
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hadronisch 


el.mag. 


Absorbermaterial 

Dicke der Absorberplatten [mm] 

Dicke der flüssigen Argonschichten [mm] 

Tiefe in charakteristischen Einheiten 


Stahl 

19 
2 ■ 2.4 
4.5 - 80 


Blei 
2.4 
2.35 
20 - 300 


A 
B 
C 
bei 


10-VGeV 

MeV] 

10-2] 

Energie [G( 


sV] 


46.1 ±0.7 
730 ± 30 
2.6 ±0.2 
3.7- 170 


11.16 ±0.05 
152 ±4 

0.64 ±0.07 
3.7- 80 



"in hadronischen Wechselwirkungslängen inklusive elektromagnetischem Kalorimeter 
''in Stahlungslängen 



Tabelle 2.3: Daten des Flüssig- Argon-Kalorimeters. Daten aus j \Hl 93<\ j 



bis 177°. Da man dort nur bei sehr kleinen Protonimpulsbruchteilen Xß hadronische 
Aktivität erwartet und das gestreute Elektron bei geringen Viererimpulstiberträgen 
(Q^ < lOOGeV^) genau in diesen Bereich fällt, ist das BEMC nur als elektromagnetisches 
Kalorimeter ausgelegt. Es besteht aus 2.5 mm dicken Bleiabsorberplatten und 
4 mm starken Szintillatorplatten, deren Licht über Photodioden ausgelesen 
wird. Die Energieaufiösung für Elektronen beträgt 



für Hadronen jedoch nur 



oe_ 
E 



E 



10% 



E/GeV 



had JE/GeV 



Instrumentiertes Eisen (Iron, 10) : Zur Rückführung des magnetischen Flusses 
befindet sich außerhalb der Spule ein Eisenjoch. Zwischen den zehn Eisenplatten 
der Dicke 75mm sind Streamerrohrkammern eingeschoben. Diese sind mit 
Pads bestückt, mit denen eine Energiemessung der nicht in den Kalorimetern 
absorbierten Teilchen möglich ist. 

Die anderen Detektorkomponenten werden bei dieser Analyse nicht direkt benutzt. 
Die Informationen sind jedoch in der Ereignisvorselektion meist über Veto- und 
Triggersignale bzw. in der Rekonstruktion als Zeit- oder Spurmessungen verwendet 
worden. Nähere Informationen zum Detektor oder zum Beschleuniger kann in der 
angegebenen Literatur gefunden werden ( ||H1 93 bp , ||ltt93| | , [[Rö95| | , | [Ros95|| ) . 



Kapitel 3 

Theoretische Grundlagen 



3.1 Tief-inelastische Lepton-Nukleon-Streuung 

Protonen haben im Gegensatz zu Elektronen eine innere Struktur, weshalb man 
zwischen elastischer und inelastischer Streuung unterscheiden muß. 

Bei der Lcpton-Nuklcon-Strcuung wird zwischen elastischer Streuung, d.h. das 
Proton bleibt intakt, und inelastischer unterschieden, d.h. das Proton zerbricht und 
die Streuung findet an einem Parton des Protons statt. 

Bei den in dieser Arbeit untersuchten tief-ineleistischen Stößen wird ein virtuelles 
Photon ausgetauscht, d.h. ein Photon dessen Viererimpulsquadrat negativ ist. 

Dieser Prozeß kann mit ep eX bezeichnet werden, wobei X den gesamten 
hadronischen Endzustand beschreibt. 

3.1.1 Kinematik 

Üblicherweise wird in der Hochenergiephysik ein Teilchen durch Vierervektoren beschrieben. 
Ein wichtiger Vektor ist der Impulsvierervektor, der definiert ist durch : 



Das Skalarprodukt eines Impulsvierervektors mit sich selbst ergibt die invariante 
Masse (mit Einsteinscher Summenkonvention) aus : 



Im folgenden werden die Ruhemassen vernachlässigt, da diese bei HERA für die 
einlaufenden Teilchen klein im Vergleich zu den jeweiligen Gesamtenergien sind. Für 





PP 



-p^ = {E + p) ■ {E -p) = m^. 



(3.1) 
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e(k) 



p(P) ^ 




i(Xi/pP). 



s X 



Abbildung 3.1: Schematische Darstehung einer Elektron Proton Streuung : 

Dargestellt werden hier die Bezeichnungen der wichtigsten Vektoren und Variablen. 



die Impulse der Strahlteilchen ergibt sich dann 

/ 820 GeV \ 

V 820 GeV / 



für die Protonen und für die Elektronen 



( 27.55 GeV ^^ 



V -27.55 GeV / 



In Abbildung ß.l| wird die tief-inelastische Streuung verdeutlicht. 

Das unter dem Winkel gestreute Elektron hat den Vierervektor k'. Das Austauschteilchen 
7 oder trägt den Impuls q = k — k'. 

In den nächsten Abschnitten betrachten wir die Mandelstam- Variablen s, t und 
u für diesen Prozeß. 



Sep = S = (P + k)' (3.2) 



ml + 2P-k + ml 



^EpE, (3.3) 
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Dies entspricht dem Quadrat der Schwerpunktsenergie -Ecms- Diese ist aufgrund 
der festen Strahlenergien konstant ^/s = 300.6 GeV. 
Für den Viererimpulsübertrag gilt 

t = q2 = (k-k')'. (3.4) 

Da t = jedoch negativ ist, wird meist 

Q' = -q'>0 

benutzt. 



Betrachten wir nur den harten Prozeß, der in der Abbildung 2?T| durch ein a 
gekennzeichnet ist, so erhalten wir nun als einlaufende Teilchen das Photon q und 
ein Parton i, d.h. ein Teil des Protons, mit dem Impuls Xt/pP. 

Dann ergibt sich die Mandelstam- Variable s zu 



s 



■yi 



(q + x,/pP)2. (3.5) 



Bei der Bestimmung s muß der Protonrest r von den Produkten des harten 
Streuprozesses getrennt werden. Rechnen wir mit dem gesamten hadronischen Endzustand 
X, so erhalten wir 

s^p = W' = (q + P)2. 

^ -Q^ + 2P-q (3.6) 

Führen wir die beiden Bjorken Skalenvariabein 

xb = X = — ^ — (3.7) 

und 

P q ^ 

Z/B^y = pT^ (3.8) 
ein, so ergeben sich die folgenden Beziehungen 

Q2 = xys (3.9) 

W ^ Q^^ = ys-Q^ (3.10) 

X 

^i/p = ^(1 + ^) (3-11) 

Ist s Null oder wesentlich kleiner als Q^, so entspricht Xi/p ~ x und x kann als 
der Impulsbruchteil des Protons betrachtet werden, den das Parton trägt. 

Analog ist y der Anteil der Elektronenenergie, der am harten Streuprozeß teilnimmt. 
Es gilt für beide Skalenvariablen x,y & [0, 1]. 
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3.1.2 Berechnung der kinematischen Größen 

Die Reaktionskinematik kann unter Berücksichtigung der Symmetrie bezüglich des 
Azimutalwinkels ip nach Formel ^]9| bei konstanter Schwerpunktsenergie s über zwei 
Größen beschrieben werden. Die Messung dieser Größen kann dabei über die Daten 
des gestreuten Elektrons und den hadronischen Endzustand erfolgen. Nun werden 
einige Methoden vorgestellt. 

Benutzt man nur das gestreute Elektronen e zur Bestimmung der Daten, erhält 
man die Größen und y über 

Ve = (3.12) 
Ql = 2EeE'^{l + cosee). (3.13) 

Die Berechnung mit den Kalorimeterinformationen über den hadronischen Endzustand 
(had) erfolgt durch 

l/had = ^^-'^ -jP-)^-'^ (3.14) 

QLd = (3.15) 

1 - Z/had 

Die Doppelwinkelmethode DA benutzt nur die Winkel des gestreuten Elektrons 
9e und den Winkel der Impulssumme über den hadronischen Endzustand öhad und 
ist somit von der Energiekalibration weitgehend unabhängig. 

, 6'had 

tan 

tan h tan — 

2 2 

Oe 

cot — 

Qda = 4E,2 2 ^317) 

, L/had , , "e 

tan h tan — 

2 2 

Andere Methoden können in der angegebenen Literatur ( ||BB94|| , ||Ros95|| ) nachgelesen 
werden. 



3.1.3 Quark Parton Modell 

Schreiben wir den Wirkungsquerschnitt der tief- inelastischen Streuung in Abhängigkeit 
der Lorentz Tensoren so ergibt sich 

da^\M\' = ^,L%W^^ (3.18) 
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mit dem Lorentz Tensor für punktförmige Teilchen gemittelt über alle Spineinstellungen 



HMs-l 



(3.19) 



e spms 



Wäre das Proton ebenfalls elementar, d.i. punktförmig, dann könnten wir für 
das Proton einen analogen Lorentz Tensor angeben. Da wir in der tief-inelastischen 
Streuung diese Annahme jedoch nicht machen können, setzen wir den allgemeinsten 
symmetrischen Lorentz Tensor an [[H1V184| 



l^M- = -W,g>''' + ^P''P'' + ^q^q'' + ^{P^q'' + q''Pn. (3.20) 

Eine Beziehung zwischen den Wi Funktionen ergibt sich durch Betrachtung der 
Stromerhaltung am Protonvertex. 



Wo 



Daraus folgt für |3.20 



P q 



^ Wo Wi 



(3.21) 
(3.22) 



w^ 

M2 



(3.23) 



Aus ß.l8| ergibt sich mit |3.19| und ß.23| 



L^/^-T^^, = 4m(k-k') + ^ [2(P-k)(P-k')-M2(k-k')l (3.24) 



Für die Lorentz Skalenvariablen Wi gilt bei punktförmigen Teilchen, d.i. wenn 



3.23| in |3.19| übergeht 



vW. 



Punkt T->Punkt / 



ÜJ] 



2m/ y - 
l^fl-i 



2Mu 



6 1 



1 



UJ 



(3.25) 
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Bei kleinem gilt, d.h. wenn wir das Proton als Teilchen mit Formfaktor G'^[Q'^) 
ausfassen : 

2MWf^'''''^ = -^G\Q^) 5(1 

^^elastisch ^ G\Q')6(^l-^j (3.26) 

Ein wesentlicher Unterschied ist, daß bei punktförmigen Teilchen der Wirkungsquerschnitt 
nur vom Verhältnis tu = bei strukturbehafteten jedoch sowohl von Q"^ als auch 
von u abhängt. 

Zerbricht nun das Proton in punktförmige Partonen, so können wir die Lepton 
Parton Wechselwirkung durch einen punktförmigen Lorentz Tensor beschreiben, wenn 
wir die Masse und den Impuls durch den jeweiligen Bruchteil x ersetzen. 

Dadurch können wir den Prozeß zweiteilen. Zum einen benötigen wir Partondichtefunktionen 
fi/p (particle density function, PDF), die die Wahrscheinlichkeit ein Parton des Typs 
i im Proton anzutreffen in Abhängigkeit der kinematischen Variablen beschreiben. 
Zum anderen bleibt dann noch die Streuung des Leptons am Parton übrig, welche. 



wie wir oben gesehen haben, einfach zu beschreiben ist (siehe auch Abbildung ^ 
Für die Partondichtefunktionen muß die Normierungsbedingung 



J2 j dxxfi^ 



li/pix) = 1 (3.27) 
gelten, da das ganze Proton aus Partonen aufgebaut ist. 



Mit der oben definierten Variable tu gilt dann aus 3.25 



Fi{uj) = xF[ = mW^f ""i^* = ^5 f 1 - — 

2XU \ XLü 

und daraus bei Summation über alle Partonen eines Protons 

F?(x) 



J dxe^Ji/pix)F;{u;) 

i 



F!{x) = ^F,{x) (3.29) 
2x 

Cj ist die Ladung des Partons i in Einheiten der Elementarladung. 

Die Funktionen Fi sind ausschließlich von einer Variablen abhängig, d.h. bei 
festem x ist F^ konstant bei verändertem Q^. Dieses Verhalten wird Bjorken Skalenverhalten 
(Bj Orken scaling) genannt. 

Der Aufbau des Protons aus nicht strukturbehafteten Teilchen kann dadurch 
nachgewiesen werden. 
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3.1.4 Quantenchromodynamik 

Bevor wir uns die Partondichtefunktionen genauer ansehen, haben wir noch ein 
Problem zu lösen : Wenn das Proton aus mehreren Partonen besteht, wieso sind 
diese Partonen nicht einzeln, sondern nur in einem Teilchen gebunden vorzufinden. 

Direkt damit gekoppelt ist das Problem, daß Teilchen existieren, die aus drei 
identischen Quarks mit gleichen Quantenzahlen aufgebaut sind (z.b. lA"*""^ >= 
\uuu > I tTT>)- Dies ist jedoch für Fermionen nach dem Pauli-Prinzip verboten 



Ber92 



Das Problem läßt sich durch die Einführung einer neuen Quantenzahl lösen. Diese 
wird in Analogie zu den drei Grundfarben Farbladung genannt, d.h. es existieren rote, 
grüne und blaue Quarks bzw. anti-rote, anti-grüne und anti-blaue Antiquarks. 

Das Proton ist, wie jedes andere Teilchen, das wir direkt beobachten können, 
farblos und somit entweder eine Kombination aus drei unterschiedlich farbigen Quarks 
(Baryonen) bzw. Antiquarks (Anti-Baryonen) oder ein Zusammenschluß aus einem 
Quark und einem Antiquark mit der zugehörigen Antifarbe (Mesonen). 

Dieses Verhalten bedingt die Existenz einer neuen Wechselwirkung, der starken 
Kraft, und der entsprechenden Austauschteilchen. Im Gegensatz zu der elektromag- 
netischen Wechselwirkung, die nur ein neutrales Austauschteilchen hat, existieren 
hier acht linear unabhängige Farb-Antifarb-Kombinationen und somit acht Austauschteilchen, 
die wir Gluonen nennen. 

Sehen wir uns nun die Verteilung von Xi/p an, so können wir verschiedene Szenarien 
unterscheiden. 

Ist das Proton ein punktförmiges, nicht strukturbehaftetes Teilchen, so nimmt an 
einer Streuung immer der gesammte Impuls P teil, d.h. Xi/p = 1 und es ergibt sich 



Bild |3^. 

Besteht das Proton hingegen aus Partonen, die sich in Richtung der Protonen 
bewegen und nicht miteinander wechselwirken, so haben diese Partonen die gleiche 
Geschwindigkeit und somit — unter Vernachlässigung des Massenunterschiedes — 
den gleichen Impuls. Für unsere Verteilung gilt dann Xi/p = ^. Bei einem Aufbau 
aus drei Quarks ergibt sich somit Bild |3.2|b. 



Wie wir jedoch bereits gesehen haben, sind bisher keine einzelnen Quarks beobachtet 
worden. Die Wahrscheinlichkeit ist also groß, daß es eine Partonwechselwirkung gibt 
und damit Impuls zwischen diesen übertragen wird. Die Xß- Verteilung hat in diesem 
Fall weiterhin ein Maximum bei |, sie ist jedoch verschmiert (Abbildung |3.2| c). 

In der Quantenelektrodynamik gibt es die Möglichkeit der Vakuumpolarisation, 
d.h. ein hochenergetisches Photon kann ein Teilchen-Antiteilchen-Paar erzeugen. 
Dies ist bei Gluonen ebenfalls möglich, wir erwarten zu den drei Quarks eine der 
Bremsstrahlung ähnliche Quark-See- Verteilung. Dadurch erhöht sich die Dichteverteilung 
bei kleinen Impulsbruchteilen (Abbildung ß.2| d). 

Der genaue Verlauf der Dichtefunktionen ist — wie bereits erwähnt — von der 
kinematischen Region abhängig und in machen Phasenraumb ereichen noch nicht 
gemessen worden. In Abbildung ß.2|e ist eine Parametrisierung der Dichtefunktionen, 
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Abbildung 3.2: 
Xi/p Verteilung für unterschiedliche Protonstrukturannahmen : a) Proton 
ist punktförmig, b) Proton besteht aus drei Quarks, c) zusätzliche Wechselwirkung 
zwischen Quarks, d) zusätzliche Seequarks, e) Parametrisierung der Messungen nach 
MRS(H) bei iJ,f = 10 GeV kPB9d^] ; siehe Erläuterungen im Text oder l\HM84\l 
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wie sie sich aus den Messungen ergibt, gezeigt. Wir erkennen, daß die Verteilungen 
ungefähr unserem Modell d entspricht, dies jedoch ein noch zu einfacher Ansatz ist. 

Integriert man über alle x und summiert über alle Valenz- und Seequarks, so 
ergibt sich beim Proton jedoch nicht 1, wie Gleichung |3.27| erwarten läßt. Es müssen 
zusätzlich noch die Gluonen berücksichtigt werden. Da die Gluonen in Seequarks 
zerfallen, ist deren Spektrum ähnlich der Seequarkverteilung (siehe auch Abbildung 
|3l^). Da die Gluonen jedoch nicht elektrisch geladen sind, koppeln sie nur stark und 
sind durch Streuung an Elektronen nicht direkt zu messen. 

Als Einführung in die Theorie der Quantenchromodynamik kann z.B. ||Ber92|| , 
HM84] oder |Nac92| empfohlen werden. 



3.1.5 Die laufende Kopplungskonstante ag 

In der Quantenelektrodynamik (QED) hat sich die Störungstheorie sinnvoll zur Berechnung 
der Wirkungsquerschnitte erwiesen. 

Mit ihrer Hilfe können wir diese in Terme in unterschiedlicher Ordnung in der 
Kopplungskonstanten a zerlegen : 

CT = Aoa° + Aia^ + ^2«^ + . . . + Aia' + ... (3.30) 



Die Grundvoraussetzung ist jedoch, daß die Kopplungskonstante wesentlich kleiner 
als eins ist und wir somit die Terme in höherer Ordnung vernachlässigen können. 

Die verschiedenen Terme Ai können dann mit Hilfe sogenannter Feynman- Gra- 
phen veranschaulicht und mit den zugehörigen Regeln berechnet werden. Näheres 
hierzu findet sich in der angegebenen Literatur ( |[HM84|| , ||1Z85|| ). 

Bei den heute möglichen genauen Messungen in der QED stellen wir fest, daß 

2 

die Kopplungskonstante a = ^ nicht konstant ist, sondern bei hohen Energien und 



damit nach ^]T| bei kleineren Abständen größer wird. So steigt der Wert von 



137 



bei Niedrigenergieexperimenten auf bei LEP, d.h. Energien in der Größe der 
Z-Boson-Masse Mz. Dies kann durch Vakuumpolarisation erklärt werden. Um die 
auftretenden Divergenzen zu umgehen wird anstelle der Ladung cq die Ladung e 






(3.31) 



verwendet. Dieser Vorgang wird Renormalisierung genannt. 

Der Wert der Kopplungskonstante ändert sich von a = a(0) auf 



IIMs4 



24 



KAPITEL 3. THEORETISCHE GRUNDLAGEN 



Wenden wir nun das gleiche Prinzip auf die QuantenchromodjTiamik (QCD) an. 
Schreiben wir den Wirkungsquerschnitt in Potenzen der starken Kopplungskonstante 
«s, so wird aus Gleichung ß.30| 



aoa° + aial + a2al + . . . + ajO;* + 



(3.33) 



Auch hier gibt es die Vakuumpolarisation, hinzu kommt jedoch noch die Kopplung 
von Austauschteilchen aneinander, die in der QED wegen der Ladungsneutralität der 
Photonen nicht auftritt. Somit renormalisieren wir die starke Kopplungskonstante 
durch 



I 



I 



I 
I 



I 



(3.34) 



T steht hier für transversale Polarisation, C für Coulomb Quarks, d.i. longitudinale 
bzw. skalare Polarisation ( |[HM84|| ). 

In Formeln ausgedrückt mit der Anzahl Quarkarten nj gilt dann 



127r 



(33 - 2nf) log 



A2 



(3.35) 



(3.36) 



wobei in der letzten Gleichung der freie QCD-Parameter 

/ -127r 



A' 



,(33-2n^)a,(/i2)^ 



(3.37) 



eingesetzt wurde. 

Wir erkennen hier, daß für die sechs uns bekannten Quarks der logarithmische 
Term im Gegensatz zur QED ein positives Vorzeichen hat. Dies bewirkt, daß die 
starke Kopplungskonstante bei kleineren Energien immer größer wird und wir somit 
nur farblose Teilchen sehen. Dieser Umstand wird color confinement genannt. 



3.1.6 Prozesse in führender und nächst-zu-führender Ordnung 

Im folgenden werden wir uns die Feynman-Graphen für einige der harten Streuprozesse 
ansehen. 

Der einfachste Graph ist in Abbildung |3.3| gezeigt. Die Streuung enthält keine 
starke Wechselwirkung, somit gehört sie in die Klasse der C(a°) Prozesse. Da Glu- 
onen elektrisch neutral sind, kann dieser Prozeß nur mit Quarks und Antiquarks als 
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Abbildung 3.4: 1+1 Ereignis in nächst zu führender Ordnung : a) Initial State 
Ahstrahlung, b) Final state Abstrahlung 
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Abbildung 3.5: virtuelle Korrekturen zu 1+1 Ereignis 



Partonen auftreten. Wir erhalten hier ein Parton zusätzlich zum Protonrest, deshalb 
wird dieser Prozeß als 1+1 Ereignis bezeichnet. 

Fügen wir in diesen Graphen eine starke Wechselwirkung ein, so existiert die 
Möglichkeit ein Gluon vom Quark abzustrahlen. Dies kann vor der Wechselwirkung 
mit dem Photon geschehen, dann wird der Prozeß initial State radiation ( ß.4| a) 
genannt. Erfolgt die Abstrahlung nach dem Stoß, so haben wir eine final state 
radiation ( |3.4| b). 

Ist das abgestrahlte Gluon weich, d.h. das Gluon hat eine geringe Energie, oder 
Quark und Gluon sind kollinear, d.h. die Impulsrichtung ist gleich, dann werden wir 
diesen Prozeß als 1+1 Ereignis klassifizieren. Da jedoch der Wirkungsquerschnitt 
für diesen Prozeß 0{as) ist, ist dies ein Beitrag zur nächst-zur-führenden Ordnung 
(NLO, Next-to-Leading Order). 

Einen weiteren Anteil erhalten wir durch virtuelle Korrekturen, d.h. Prozesse, die 
ein internes virtuelles Teilchen austauschen (Abbildung ß.5| ). Der Wirkungsquerschnitt 
dieser Prozesse selber ist zwar (9(a^), jedoch läßt sich der Endzustand nicht vom 
Prozeß in führender Ordnung unterscheiden. Durch die deshalb notwendige Addition 
der Matrixelemente entsteht bei der Berechnung des Wirkungsquerschnittes ein Interferenzterm 
in Oias). 



+ 



+ 




(3.38) 



Anzumerken ist hier noch, daß beim Prozeß ß.5| a beim Stoß zwischen Parton und 
Photon nicht der komplette Partonimpulsbruchteil x zur Verfügung steht, da das 
virtuelle Gluon ein Teil des Partonimpulses erhält. 

Gehen wir nun zur nächsten Prozeßklasse über. Wird entgegen unserer bisherigen 



Annahme ein hartes Gluon beim NLO 1+1 Ereignis (Abbildung p.4| ) abgestrahlt, so 
erhalten wir zwei harte Partonen zum Protonrest. Dieser Prozeß wird in Analogie zur 
Compton Streuung in der Quantenelektrodynamik QCD-Compton Prozeß genannt 
(Abbildung |3.6| a). Der Graph trägt somit zur führenden Ordnung der 2+1 Ereignisse 
bei und ist ö{a^. Die Entscheidung, wann ein Prozeß zur 2+1 oder zur 1+1 Ereignisklasse 



3.2. MONTE-CARLO-SIMULATIONEN 



27 




Abbildung 3.7: 2+1 Ereignis in nächst-zu- führender Ordnung 

gehört, muß durch einen Abschneideparameter getroffen werden. Diese Unterscheidung 
ist auch deshalb wichtig, da der Wirkungsquerschnitt für die Abstrahlung kollinearer 
bzw. weicher Partonen divergiert. Eine genauere Betrachung wird bei der Behandlung 
des analogen Problems in Monte- Carlo- Generatoren in Kapitel |3.2| durchgeführt. 

Bei der sogenannten Boson-Gluon-Fusion (Abbildung |3.6| b) wird ein Gluon vom 
Proton abgestrahlt. Das Gluon zerfällt dann in ein Quark- Antiquark-P aar. Die 
Streuung des Photons kann entweder am Quark oder am Antiquark erfolgen. Auch 
dieser Prozeß gehört zur Klasse der LO 2+1 Ereignisse. 

Die NLO 2+1 Prozesse entstehen analog den 1+1 Ereignissen durch Abstrahlung 
weicher Gluonen oder durch Interferenzterme bei den virtuellen Korrekturen (Abbildung |3.7| ). 
Dies sind Beiträge zur Störungsrechnung in ö{a^). 



3.2 Monte-Carlo-Simulationen 

Um sich ein Bild von den zu erwartenden Detektordaten zu machen, ist es notwendig 
Monte-Carlo-Simulationen einzusetzen. 

Die Simulation von Ereignissen geschieht hierbei in mehreren Schritten. Einen 
allgemeinen Uberblick über die unterschiedlichen Stufen eines Monte-Carlo-Generators 



gibt Abbildung ß.8| . Begonnen wird mit dem harten Streuprozeß, dessen Art und 



Kinematik den theoretischen Modellen entnommen wird. Nachdem zusätzliche Abstrahlungen 
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Partonen Hadronen Detektor Objekte 



p 




harte QCD-Streuung Fragmentation Simulation 

Partonschauer Teilchenzerfall Rekonstruktion 



pertubativ nicht-pertubativ 

Abbildung 3.8: Übersicht über die Vorgänge bei Monte-Carlo-Generationen 

von Partonen berücksichtigt worden sind, wird die Hadronisierung durchgeführt, d.h. 
die farbigen Quarks werden zu farblosen Hadronen zusammengefaßt. Der Zerfall 
kurzlebiger Teilchen muß dann ebenso beachtet werden wie Detektoreffekte, die 
durch unsensitve Detektorregionen oder die unterschiedlichen Detektoreffizienzen 
entstehen. 

Im folgenden werden die Schritte kurz erläutert. Genaueres ist der Literatur über 
Monte-Carlo-Generatoren bzw. den Programmanleitungen zu entnehmen ( ||Sey95|| , 



Ham93| 


1 


lng93 


1 


MKSW93 


1 


LÖ94 



3.2.1 Beschreibung auf Parton-Niveau 

Wir haben bereits im vorigen Abschnitt die einzelnen, für uns interessanten Prozesse 
betrachtet. Anhand dieser Wirkungsquerschnitte wird dann bei den Simulationen 
ein Ereignis gewürfelt. 

Leider gibt es keinen Ereignisgenerator, der Ereignisse auf der Basis von NLO 
Berechnungen erzeugen kann. Alle vorhandenen Programme benutzen Matrixelemente 
in führender Ordnung (LO) für die Erzeugung von 2+1 Ereignissen. 

Um dennoch die Korrekturen berücksichtigen zu können, die durch Gluonabstrahlung 
in höheren Ordnungen entstehen, werden als Ausweg QCD-inspirierte Partonschauer 
(PS, parton shower) verwendet. Diese Partonschauer müssen von den Abstrahlungen 
in höheren Ordnungen der QCD unterschieden werden, da die Schauer unabhängig 
vom harten Subprozeß angewandt werden. Es treten somit keine Interferenzen oder 
virtuelle Korrekturen auf. Partonschauer haben sich zusammen mit der im nächsten 
Abschnitt erläuterten Hadronisierung als nützliches Hilfsmittel bei der Beschreibung 
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der Daten erwiesen. 

Wir müssen hierbei zwischen zwei Fällen unterscheiden : Zum einen kann die 
Gluonabstrahlung vor der Streuung erfolgen (initial State parton shower), zum anderen 
danach (final state parton shower). 

Bei den final state Partonschauern treten zwei Divergenzen auf, d.h. an diesen 
Stellen im Phasenraum wird die Abstrahlwahrscheinlichkeit eins bzw. der Wirkungsquerschnitt 
für diesen Prozeß unendlich. Eine Divergenz erhalten wir bei Kollinear ität, d.h. 
wenn der Impuls des abgestrahlten Gluons in die gleiche Richtung zeigt wie der des 
abstrahlenden Partons. Dieses Problem können wir durch einen Abschneideparameter 
beseitigen, der z.B. eine minimale invariante Masse des Gluon Parton Paares verlangt. 
Der andere Fall tritt bei weicher Abstrahlung auf, d.i. wenn das Gluon eine sehr 
geringe Energie erhält. Auch dieses Problem läßt sich software-technisch umgehen, 
wobei es jedoch mehrere Vorgehensweisen gibt, deren Erklärung hier zu weit führen 
würde. 

Bei initial state Partonschauer ist die Situation noch wesentlich komplizierter. 
Das größte zusätzhche Hindernis ist die Veränderung der Energie des Partons durch 
die Abstrahlung und damit der Kinematik des harten Prozesses. Dies können wir 
durch Rückwärtsentwicklung umgehen, d.h. es wird nicht die Energie des Partons 
nach Abstrahlung verringert, sondern die vor der Abstrahlung vergrößert. 

Es ergeben sich weitere Probleme, wenn wir mehrfache Gluonabstrahlung berücksichtigen. 
Die verschiedenen Modelle haben unterschiedliche Algorithmen, die z.B. die Reihenfolge, 
in der die Gluonen entstehen, und den Abschneidepunkt, wann dieses Verfahren 
abzubrechen ist, betreffen. 

Das Ariadne Monte-Carlo-Programm benutzt mit dem Colour Dipole Model ein 
völlig anderes Schema, das Partonschauer komplett umgeht. Die einzelnen Farb- 
Antifarb-Paare werden als Dipole verstanden, die analog der Elektrodynamik durch 
Polarisation neue Dipole hervorrufen. 

3.2.2 Die Hadronisierung 

Bei der Hadronisierung werden aus den farbigen Partonen farblose Hadronen gebildet. 
Es ist bisher noch nicht gelungen diesen Vorgang theoretisch zu beschreiben, da hier 
der pertubative Bereich verlassen wird und die Störungsrechnung somit nicht mehr 
funktioniert. 

Die beiden meistbenutzten Modelle sind die Cluster Fragmentation und die String 
Fragmentation. 

Beim Cluster Modell werden jeweils zwei im Phasenraum benachbarte Partonen 
entsprechend ihren Farbladungen zu einem sogenannten Cluster kombiniert, der dann 
entsprechend seiner Energie, Ladung und der vorhandenen Quarkarten in Hadronen 
zerfällt. 

Das String Modell hingegen benutzt ein sehr anschauliches Analogon. Es zieht 
zwischen der Farbladung und der zugehörigen Antifarbladung eine "Saite", die analog 
den Feldlinien in der Elektrodynamik Energie enthält. Aufgrund des color confinements 
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ist die Energie hier jedoch proportional zur Länge der Feldhnien. Wird die Saite zu 
lang und damit die Spannung zu groß, so reißt diese und es entsteht ein neues Quark 
Antiquark Paar. Die Saiten, die am Ende des Verfahrens vorhanden sind, bilden 
dann die Hadronen. 

Die Probleme in beiden Modellen liegen in der Erzeugung von Baryonen, da diese 
nicht aus Quark Antiquark Paaren, sondern aus jeweils drei Quarks bestehen. Leider 
versagen die Modelle hier und nur durch einige Tricks und eine große Anzahl freier 
Parameter ist es im String Modell gelungen, eine Beschreibung der existierenden 
Daten zu erreichen. 

Im abschließenden Schritt auf Hadronniveau müssen die Teilchenbahnen nun im 
Magnetfeld berechnet und der Zerfall kurzlebiger Komponenten in Betracht gezogen 
werden. Dies ist durch die Kenntnisse der Zerfallszeiten und der verschiedenen Kanäle 
jedoch verglichen mit den vorhergehenden Schritten einfach und durch die Theorie 
beschrieben. 



3.2.3 Beschreibung auf Detektor-Niveau 



Im letzten Schritt werden die Detektoreffekte simuliert. Die danach vorhandenen 
Daten entsprechen denen gemessener Ereignisse. Wir haben nur noch solche Werte, 
die z.B. der Ladungsverteilung in den Kalorimetern oder den Treffern in den Spurkammern 
entsprechen. 

Diese Stufe ist unabhängig von den Monte- Carlo- Generatoren und wird von GE- 
ANT ||GEA93 | mit speziellen Hl Anpassungen durchgeführt [|H1 95a . 



Die Daten lassen sich nun mit den wirklichen Meßergebnissen des Detektors 
vergleichen. Um dies einfacher zu gestalten werden im abschließenden Rekonstruktionsschritt 
Teilchenidentifikationen, Energieumrechnungen, Klassifizierungen und andere Berechnungen 



durchgeführt [[Hl 95c| • 



Kapitel 4 

Messungen mit Jet-Raten 



4.1 Jets 

Wir werden nun Messungen durchführen und dann mit den dabei gewonnenen Detektordaten 
die Theorie testen. Dazu ist es notwendig eine Verbindung zwischen dem Detektorniveau 
und dem Partonniveau zu suchen. 

Einen Uberbhck über die einzelnen Niveaus gibt Abbildung ^?T| (siehe auch Kapitel 



Wenn wir erst den Zwischenschritt über das Hadronniveau machen, so können wir 
erkennen, daß aus einigen wenigen Partonen eine wesentlich größere Anzahl Hadronen 
entstanden ist. Ein einfaches Modell, das eine Lokale Parton Hadron Dualität 
(LPHD, local parton hadron duality) vorhersagt, ist dabei erstaunlich erfolgreich. 
Die Grundidee ist hierbei, daß die Anzahl der Hadronen in einem Phasenraumbereich 
proportional zur Anzahl der Partonen ist Pey95|| . Die Aussagekraft des Modells ist 



zwar stark begrenzt, aber es zeigt uns, daß wir mehrere Hadronen zu einem Gebilde, 
das wir Jet nennen, zusammenfassen müssen. 

Auf Detektorniveau tritt das gleiche Phänomen auf, da die Hadronen auf ihrem 
Weg durch den Detektor an unterschiedlichen Stellen Energie deponieren und nicht 
immer sofort gestoppt werden. Aufgrund der begrenzten Raumaufiösung der Detektorzellen 
können wir jedoch die genaue Anzahl der Ladungsdepositionen kaum erfassen. Es 
ist deshalb sinnvoller Bereiche mit hohem Ladungsinhalt zu sogenannten Clustern 
zusammenzufassen . 

Einen Uberblick über die Anzahl Partonen, Hadronen und Cluster pro Ereignis 



gibt Abbildung 4.2 



Wir müssen somit einen Algorithmus entwerfen, der Hadronen oder Cluster zu 
Jets zusammenfaßt. Sind zwei Partonen eng benachbart, so kann es vorkommen, daß 
beide Partonen auf Hadronniveau nur einen Jet ergeben. Deshalb ist es sinnvoll den 
gleichen Algorithmus auch auf die Partonen anzuwenden und diese Fehlinterpretationen 
im Wirkungsquerschnitt zu berücksichtigen. Dies haben wir im Kapitel p.l.6| jedoch 
schon getan. Der Begriff Objekt beschreibt im folgenden somit Partonen, Hadronen 
bzw. Cluster in Abhängigkeit des betrachteten Niveaus. 
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Abbildung 4.1: Graphische Darstellung der unterschiedlichen Niveaus :Zu 

sehen sind im äußeren Bereich die beiden zentralen Kalorimeter (LAr, BEMC) mit 
der Energieverteilung in den Clustern (dunkle Punkte, Detektorlevel). Im inneren 
ist ein Graph eines Prozesses gezeichnet, der diese Messung hervorgerufen haben 
könnte. Die Strahlteilchen laufen von den Seiten in den Detektor ein und streuen 
aneinander. Es ergeben sich der Protonrest und zwei harte Partonen (Partonlevel, 
nach Partonschauern). Diese bilden dann Hadronen, die durch die Auffächerung am 
Ende der Partonlinien verdeutlicht werden (Hadronlevel) . 
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Abbildung 4.2: Anzahl Partikel pro Ereignis auf den folgenden Niveaus : 

a) Parton, b) Hadron, c) Cluster. Dargestellt sind alle Ereignisse des MEHQ 



Datensatzes (siehe Tabelle \571\ ). 



Einige Jetalgorithmen werden nun vorgestellt. Zur Berechnung der NLO Wirkungsquerschnitte 
existieren nur die beiden Programme PROJET ||Gra93|| und DISJET |[BM94|] . Da 
in diesen nur der JADE Algorithmus eingebaut ist, sind wir gezwungen diesen zu 
verwenden. Für die Zukunft sind verschiedene neue Programme angekündigt bzw. 
befinden sich zur Zeit in der Testphase, die mehrere verschiedene Algorithmen enthalten. 
Daher beschränken wir uns hier nicht auf den JADE Algorithmus. 



4.1.1 Cone Algorithmus 

Der Cone Algorithmus ist sicher der intuitivste, da er räumlich benachbarte Partikel 
zu einem Jet zusammenfaßt. Die Abstandsberechnung geschieht hierbei über den 
Azimutalwinkel (p und die Pseudorapidität rj, die mit dem Polarwinkel über 



(4.1) 



zusammenhängt. Der Algorithmus startet mit einer Anzahl Jet Kandidaten, die 
Partikeln mit einer transversalen Energie größer als das einstellbare Minimum entsprechen. 
Dann werden alle Partikel hinzugefügt, die in einem Radius 




R 



ÖT]'^ + öip'^ < R„ 



(4.2) 



um den Kandidaten mit dem höchsten transversalen Impuls liegen. Dies wird dann 
der Reihe nach mit allen Kandidaten wiederholt. Am Ende werden alle Kandidaten, 
die eine bestimmte transversale Mindestenergie besitzen, als Jet aufgefaßt. 
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4.1.2 JADE und verwandte Algorithmen 

Bei den Clusteralgorithmen wird nicht der räumliche Abstand als Maß verwendet, 
sondern die invariante Masse zweier Partikel. Meist wird die invariante Masse jedoch 
durch eine Skala geteilt, um unabhängig von der Energie des stattfindenden Prozesses 
zu werden. Bei e^e~ Collidern wird die Schwerpunktenergie E"^^ = verwendet, 
bei HERA verwenden wir die invariante Masse des hadronischen Endzustandes W"^. 
Eine Begründung hierzu kann in |[Nis94b|| Kapitel 5.1 gefunden werden. Es ergibt 
sich dann als Schnittgröße 

» - (4.3) 

Der Algorithmus besteht aus drei Schritten, die bis zur Erfüllung der Abbruchbedingung 
wiederholt werden. 

1. Für alle möglichen Objektpaare wird die invariante Masse rriij berechnet. 
Es wird der niedrigste Wert rriki bestimmt. 

2. Ist — ^ > Heut, so wird das Verfahren abgebrochen und die übrig gebliebenen 
Objekte entsprechen den Jets. Ansonsten wird mit der Berechnung fortgefahren. 

3. Die Teilchen k und /, werden zusammengefaßt und das Summenteilchen wird 
an ihrer Stelle in die Liste aufgenommen. Der Algorithmus fängt nun wieder 
mit Schritt 1, der Berechnung der invarianten Massen, von vorne an. 

Für die Rekombination zweier Teilchen gibt es keine vorgeschriebene Regel, da 
in der Theorie mit masselosen Objekten gerechnet wird. In der Praxis haben die 
Partikel jedoch eine nicht verschwindende Masse. Zusätzlich ergibt sich bei Addition 
der Vierervektoren eine Masse, die nicht der Summe der Massen der Ausgangsteilchen 
entspricht. Die unterschiedlichen Ideen haben zu verschiedenen Algorithmen geführt. 



die in Tabelle aufgelistet sind. 

Im Gegensatz zum Cone Algorithmus wird hier keine minimale transversale Energie 
für die Jets verlangt. Somit sollte ein Jet, der durch den Protonrest (remnant) 
entsteht, gefunden werden. Dies geschieht auch auf Parton- und Hadronniveau. Auf 
dem Detektorniveau hingegen verläßt der größte Teil des Remnants den Detektor 
durch das Strahlrohr. Aus diesem Grunde werden die Cluster um ein Pseudoteilchen 
erweitert, dessen transversaler Impuls zu Null gesetzt wird. 

Pps.,z = Pz + k^- Phz (4.4) 

Cluster 

Px=Py = 

E = b1. 



wobei die Summe über alle Cluster einschließlich derer des gestreuten Elektron läuft. 
Zusätzlich ist die Kenntnis über den gesamten hadronischen Endzustand, inklusive 
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Schema 


Berechnung von mf^ 


Rekombination 


Anmerkungen 


JADE 


2E,E,(l-cos(%)) 


Pfc = Pi + Pj 


Massen werden bei m'f, 
Berechnung vernachlässigt 


E 


(Pi + P,)' 


Pk = Pi + Pj 


lorentz-invariant 


EO 


(Pi + P,)' 


Ek = (Ei + Ej) 
Pk= + 


p ist nicht erhalten 


P 


(Pi + Pj? 


Pk=Pi+ Pj 
Ek = \pk\ 


E nicht erhalten 



Tabelle 4.1: Rekombinationsschemata von JADE und verwandten 
Algorithmen 



q 



^ 1 ^ V xaaaXww^ 



Abbildung 4.3: Schematische Darstellung der Kinematik im Breitsystem 



Protonrestes, zur W"^ Berechnung erforderlich. Weiteres zum Pseudoteilchenansatz 
ist in | GM91 | und [ Nis94b | zu finden, der JADE Algorithmus ist in [ JAD86|| näher 
beschrieben. 



4.1.3 kr Algorithmus 



Der oder Durham Algorithmus ( ||CDB9^ , Pes93|| ) gehört auch in die Klasse der 
Clusteralgorithmen. 

Der Algorithmus benutzt die Partikel nicht im Labor- sondern im Breitsystem, in 
dem das an der Streuung teilnehmende Parton durch das Photon wie an einer Wand 
zurückgeworfen wird (siehe auch Abbildung [4.3|) . 



Wir können uns das Verfahren zweigeteilt vorstellen. Im ersten Schritt wird der 
Protonrest von den übrigen sogenannten Macro Jets separiert und diese Macro Jets 
werden dann im abschließenden Schritt aufgetrennt. 

Der erste Schritt kann dabei folgendermaßen beschrieben werden : 

1. Zuerst wird eine Skala Et definiert, die die harte Streuung umschreibt. 



Q2 > ^2 > ^2 
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2. Dann wird der Schnittwert für alle Objekte i relativ zum einlaufenden Proton 

_ 2(l-cos(g,,)) 

und für alle Objektpaare 

_ 2(1 -C0S(%)) ■ p2N 

^ min(i.^ , i/j. j 

bestimmt. 

3. Es wird das Minimum aller {yip, yij} berechnet und in Abhängigkeit des Ergebnisses 
einer der folgenden Schritte eingeleitet : 

• Sind alle Schnittwerte größer 1, so wird das Verfahren abgebrochen und 
mit dem zweiten Teil des Algorithmus begonnen. 

• Ist ein Uip kleiner als alle anderen Schnittwerte, so wird das Objekt i in 
den Beamjet aufgenommen, d.h. für die weitere Berechnung verworfen. 

• Ist ein Uij kleiner als alle anderen Schnittwerte, so werden die beiden 
Objekte i und j zu einem neuen Partikel zusammengefaßt. 

4. Die Schritte 2 und 3 werden solange wiederholt, bis die Abbruchbedingung 
mm{{yip, Uij}) > 1 erfüllt ist. Die übrigbleibenden Partikel werden als Macro 
Jets bezeichnet. 

Im zweiten Teil des Verfahrens werden nun nur noch die nicht zum Beamjet 
gehörenden Objekte verwendet, d.h. die im ersten Teil entstandenen Macro Jets 
werden wieder in die Ursprungscluster zerlegt. 



1. Wir definieren nun einen Schnittwert 



ycut = < 1. 



2. Für alle Okjektpaare (i, j) wird nun analog zum ersten Teil t/ij berechnet. 

3. Ist das Minimum aller yij kleiner als unser Auflösungsparameter ycut, so wird 
das Verfahren fortgesetzt, ansonsten ist der Jetalgorithmus hier beendet und 
die verbhebenen Partikel ergeben die gesuchten Jets. 

4. Die beiden Objekte, die das minimale yij bilden, werden rekombiniert. Die 
verschiedenen Schemata sind analog denen des JADE Algorithmus definiert. 
Danach wird das Verfahren ab Schritt 2 wiederholt. 



4.2. DIE MESSUNG DER STARKEN KOPPLUNGSKONSTANTE 
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4.2 Die Messung der starken Kopplungskonstante 

Wenden wir uns nun der Messung eines der entscheidenden QCD-Parameter zu. Die 
Kopplung der starken Wechselwirkung wird durch die laufende Kopplungskonstante 
beschrieben. Wie wir schon in Kapitel ß.l.5| gesehen haben, können wir den 



Wirkungsquerschnitt durch die Störungsrechnung in Potenzen der Kopplungskonstante 
beschreiben. 

Tun wir dies nun bis zur nächst-zur-führenden Ordnung in Abhängigkeit der 
Jetanzahlen, so erhalten wir : 

rfcTi+i((5^a,,?/cut) = a2o(Q^) + A2)a2l(Q^?/cut) (4.5) 

da2+l(Q^a„?/cut) = a,(Q^A2)a3i(Q^?/cut) + a.(Q^A^)^a32(Q^^/cut) (4.6) 

Hierbei lassen sich die CL-ij clU.S den Wirkungsquerschnitten berechnen. Sie beschreiben 
Ereignisse mit i Jets, d.h. {i — 1) harten Jets und dem Protonrest, und sind Beiträge 
in Ordnung a^. 

Da wir die Wirkungsquerschnitte nicht direkt messen, sondern nur die Anzahl an 
n + 1 Jetereignissen bestimmen können, ist es sinnvoll mit der Jetrate zu rechnen. 
Außerdem erreichen wir dann Unabhängigkeit von Fehlern in der Luminositäts- und 
Effizienzenbestimmung. 

Die 2+1 Jetrate ist definiert durch : 

R... - §±i._Jy-^ (4.7) 
a,(g2, A2)a3i(Q2, yeut) + a'(Q', A')a32(Q', Vcnt) 



a20(<5^) + «s(<5^ A2)[a2l(g2,ycut) + ä3l(Q2,?/cut)] 



(4.^ 



N2+1 und äsi trägt der Möglichkeit Rechnung, daß manche 2+1 Ereignisse zwar 
zum totalen Wirkungsquerschnitt gerechnet werden, aufgrund von Phasenraumschnitten, 
z.B. im Jetwinkel, jedoch nicht zur 2+1 Jetanzahl. 

Wir können diese Gleichung nun nach aufiösen und erhalten dann 

asiQ^) = ^[ai,•(Q^ycut),Ä2+l(Q^l/cut)]. (4.9) 

Damit können wir durch Messung der Jetrate und Berechnung der Koeffizienten 
ttij bei verschiedenen Q^, die Stärke der starken Kopplung und deren Laufen mit 
bestimmen. 

Einsichtig ist, daß wir die Faktoren a^j nur von der Theorie, d.h. auf Partonniveau, 
berechnen können. Daher müssen wir dafür sorgen, daß unsere auf Detektorniveau 
gemessene Jetrate der auf Partonniveau entspricht. Ansonsten müssen wir einen 
Korrekturfaktor 

R = (4.10) 

-n^2+l,detektor 
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einführen. 

Mit Hilfe eines Fitprogrammes läßt sich somit nun unser Meßwert bestimmen. 
Dies wurde zum ersten Mal von der Hl CoUaboration durchgeführt ([ Nis94b|| ) und ist 
inzwischen auch von anderen Experimenten erfolgreich angewandt worden |[ZEU95|] . 



4.3 Die Messung der Gluondichtefunktion 

Nahezu analog gestaltet sich die Berechnung der Gluondichtefunktion fg/p{x). 

Hier beruht die Messung auf der Aufteilung des Wirkungsquerschnittes in einen 
quark-induzierten und einen gluon-induzierten Teil. 



er = / dxi/p 



fq/p{xi/p, fi'^)ag{xi/p, ^'^) + fg/p{xi/p, ^^)ag{xi/p, (4.11) 



Um die Messung direkt auszuführen ist es notwendig, QCD-Compton Ereignisse 
(Abbildung ß.6| a) von Boson-Gluon- Fusion Ereignissen (Abbildung |3.6| b) zu trennen, 
da wir so die Integrale über die beiden Terme einzeln messen können. Um dann 
die Dichtefunktion zu bestimmen, müssen wir die einzelnen Wirkungsquerschnitte 
berechnen und eine Parametrisierung der Quarkdichtefunktionen einsetzen. Die 
Dichtefunktionen der Quarks sind wesentlich genauer bekannt als die der Gluonen. 

Auch hier kann man durch Teilen der Gleichung [4.11| durch die Effizienzen und die 



Luminosität zu unseren Meßgrößen, den Jetanzahlen gelangen. Für dieses Verfahren 
ist also eine verläßliche Jeterkennung und die richtige Behandlung der Partonjet- 
Detektorjet-Korrelation wichtig. 

Bei einer NLO Berechnung können wir die oben beschriebene Aufteilung nicht 
mehr vornehmen. Dies ist an Abbildung ß.7|b zu erkennen. Wir können den Prozeß als 



Boson-Gluon-Prozeß rechnen, wenn wir, wie in dem Feynman-Graphen angedeutet, 
das Gluon zum Quark fassen, oder wir können das Ereignis zu den QCD-Compton 
Ereignissen zählen, wenn das aus dem Paarbildungsprozeß entstandene Antiquark 
kollinear zum Protonrest ist und wir es somit zum Remnant packen. 

Die einzige Möglichkeit ist also, das komplette Integral zu fitten. Um den Fit 
durchzuführen, müssen wir jedoch die Wirkungsquerschnitte für jeden Fitschritt neu 
berechnen. Dies ist mit der heutigen Rechnertechnologie nicht möglich. 

Ein Ausweg bietet die Methode der Meilin Transformation. Gehen wir durch eine 
solche Transformation in den Momentenraum über, so läßt sich dort das Moment 
des Wirkungsquerschnittes als Produkt der Momente der Partondichtefunktionen fi/p 
und der Momente der Einzelwirkungsquerschnitte cij für den aktuellen kinematischen 
Bereich berechnen. Führen wir nun den Fit nur in bestimmten vorgegebenen Phasenraumintervallen 
durch, so müssen die Momente der Einzelwirkungsquerschnitte an den Stützstellen 
nur einmal berechnet werden. Da dies der zeitaufwendige Teil ist, die Transformation 
der Dichtefunktionen und die Rücktransformation der Momente aber schnell erfolgen 
kann, haben wir hiermit eine Möglichkeit gefunden den Fit in nächst-zu-führender 
Ordnung durchzuführen [|GHVB9^|. 



Kapitel 5 

Ergebnisse der Untersuchungen 



5.1 Datensätze 



Wir wollen nun die Korrelation zwischen den Jets auf Parton- und denen auf Detektorniveau 
genauer untersuchen. Dazu benötigen wir die schon angesprochenen Monte-Carlo- 
Simulationen. In dieser Arbeit verwenden wir Datensätze, die mit den Programmen 

IKSW93[| und Ariadne 4.3 ||Lö94|| erzeugt wurden. 



Lepto 6.1 |lng92 



Herwig 5.8 

Meist wird eine Unterteilung in einen niedrigen und einen hohen Bereich vorgenommen, 
da das Elektron dann ins rückwärtige Kalorimeter (Q^ < 100 GeV^) bzw. ins Flüssig- 
Argon-Kalorimeter (Q^ > 100 GeV^) gestreut wird. Die einzelnen Datensätze sind 



in Tabelle p.l| aufgeführt 



Zum Vergleich mit den Daten wurden die vorselektiertenQ Meßwerte aus der 
Datennahme des Jahres 1994 verwendet. Um eine Verfälschung der Daten gegenüber 
den zur Verfügung stehenden Monte-Carlos zu verhindern, wurden nur Positron Runs 
mit nominellem Vertex verwendet. Es gibt sich somit eine zur Verfügung stehende 
integrierte Luminosität von C = 2.74 pb~^. Die Daten sind ebenfalls in die beiden 

Bereiche aufgeteilt. Genaueres ist der Tabelle |5.1| zu entnehmen. 



5.2 Standardschnitte 

Die in diesem Kapitel erläuterten Schnitte werden benötigt, um zum einen möglichst 
viele Nicht-DIS Ereignisse zu verwerfen und zum anderen spezielle Eigenarten der 
Hard- und Software zu berücksichtigen. 



^Die Vorselektion wurde von der HI-ELAN Gruppe vorgenommen, eine genauere Selektion auf 
DIS Ereignisse von der Jet-Gruppe. 
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Generator 


Anzahl Ereignisse 


Simulationsversion 


Kürzel 


Datei 


Format 


Rekonstruktionsversion 


Bereich 


Lepto 6.1 
MEPS3621 


29998 
DSTQ 


3.06/06 
6.00/20 


MEPSHQ 
g2 > 90GeV2 


Lepto 6.1 


99031 


3.06/06 


MEPSLQ 


MEPS3620 


DST 


6.00/20 


g2 > OGeV^ 


Ariadne 


22715 


3.05/04 


ARHQ 
g2 > lOOGeV^ 


MEAR3652 


V2g 


5.02/18 


Ariadne 


87248 


3.05/03 


ARLQ 


MEAR3641 


DST 


6.00/12 


Q2 > 6GeV^ 


Herwig 5.8 


34998 


3.06/27 


HEHQ 


MEPS3604 


DST 


6.01/08 


> 90GeV2 


Herwig 5.8 


63039 


3.06/14 


HELQ 


MEPS3602 


DST 


6.01/05 


Q2 > eGeV^ 


Generator 


Anzahl Ereignisse 


Format 


Kürzel 


Daten 


21651 


DST6 


DAHQ 


Daten 


162892 


DST6 


DALQ 



°DST bedeutet Data Summary Tape 
''Version 2 der Files, Hl Standardformat 



Tabelle 5.1: Liste der verwendeten Datensätze 



5.2.1 Kinematische Variablen 



Wie in Kapitel ß . 1 . 2| erläutert , benötigen wir für die Berechnung der Reaktionskinematik 



Kenntnis über nur zwei unabhängige Variablen. Wir wollen nun die verschiedenen 
Berechnungsmethoden genauer betrachten. 



In Abbildung 5.1 ist die Qualität der Q Bestimmung für vier Methoden gezeigt. 



Auf der Abszisse ist das vom Generator zur Erzeugung des Ereignisses benutzte 
negative Quadrat des Viererimpulsübertrages Ql^c aufgetragen, auf der Ordinate 
jeweils der auf Detektorniveau rekonstruierte Wert. Verwendet wurden die Elektron- 
(Qg), die Doppelwinkel- {Qda)^ die Sigma- (Ql) [[BB94|| und die Jacquet-Blondel- 



Methode (Qjb). Gezeigt sind die Daten der Lepto- Monte-Carlos (MEPSHQ, MEPSLQ). 
Die Größe der Boxen in diesen Diagrammen gibt jeweils die Anzahl der Ereignisse 
in diesem Bereich an. Der Maßstab ist hierbei logarithmisch, so daß auch einzelne 
Einträge noch eine kleine Box erzeugen. Insgesamt erkennen wir, daß die Rekonstruktion 
mit allen Methoden gut funktioniert. Die Jacquet-Blondel Methode ist von der 
Auswertung des hadronischen Endzustandes besonders abhängig. Dadurch haben 
Detektorverluste und die hadronische Energiekalibration einen großen Einfluß. Deshalb 
sind hier auch die größten Abweichungen erkennbar. Die Sigma Methode ist der 
Jacquet-Blondel Methode sehr ähnlich, durch Verwendung der Größe E — des 
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totalen Ereignisses im Gegensatz zur zweifachen Strahlelektronenergie ist die Abhängigkeit 
von der Bestimmung des hadronischen Endzustandes geringer. Außerdem berücksichtigt 
sie damit Photonabstrahlungen am einlaufenden Elektronast. Diese Abstrahlungen 
sind in unseren Datensätzen nicht vorhanden und die Sigma Methode kann ihren 
Vorteil nicht voll ausspielen. Die Elektronmethode, die vollkommen unabhängig vom 
hadronischen Endzustand ist, ist über den gesamten kinematischen Bereich gesehen 
die geeigneteste. 

Die Bestimmung der zweiten, von unabhängigen Variable yß kann ebenfalls 
über die obengenannten Methoden erfolgen. Die entsprechenden Verteilungen sind 
für die hohen in Abbildung |5.2| , für die niedrigen in Abbildung |5.3| gezeigt. 

Bei hohen zeigt die i/b Verteilung ebenfalls eine gute Ubereinstimmung bei 
allen Methoden. Bei i/b größer 0.7 werden jedoch auch die Abweichungen größer. 
In diesem Bereich ist der Untergrund durch Photoproduktionsereignisse groß, da 
das gestreute Elektron den Detektor unbeobachtet durch das Strahlrohr verläßt 
und niederenergetische Hadronen ein Elektron mit hohem yB vertäuschen können. 
Deshalb betrachten wir im weiteren nur Ereignisse mit yß < 0.7. Vergleichen wir 
die unterschiedlichen Methoden, so ist auch hier die Elektronmethode die genaueste. 
Nur in Bereich yß < 0.2 erhalten wir mit der Doppelwinkelmethode die höchste 
Genauigkeit. 

Für kleiner 100 GeV^ fällt auf, daß kaum Einträge mit yß > 0.6 vorhanden 
sind. Dies ist durch den Zusammenhang zwischen yß, der Energie des gestreuten 
Elektrons E'^ und dessen Winkel 6e durch die Gleichung O erklärbar. Der Rand 
des BEMC bei 151° und die minimale Clusterenergie des Elektronkandidaten von 
10 GeV ergeben dann eine Abschätzung für den Grenzwert. 



P ■ (k - k') 

?/iJ,max 



P k 



2Ee-^:(l-COs(öe)) 



2E 

max '^-^^ 



-ge,mm(l " COs(ge,min)) _^ 10(1 - COs(15r)) 

2E, ~ 2 ■ 27.5 



0.66 (5.1) 



Die geringe Ereignisanzahl mit yßC > 0.6 entsteht ausnahmslos durch Ereignisse 
mit generiertem > 100 GeV, die durch das Fehlen eines oberen Schnittes 
im MEPSLQ Monte-Carlo entstehen. Diese Elektronen fallen deshalb nicht in den 
Bereich des BEMC und die oben angegebene Schranke gilt nicht. Eine Aufteilung und 
Untersuchung der gestreuten Elektronen in Abhängigkeit der Detektorkomponente 
wird in Kapitel p.'I.'S^ vorgenommen. Ansonsten erkennen wir auch hier ein Verhalten 
analog dem bei hohen Werten. 

Zusammenfassend können wir also aus der Rekonstruktion der kinematischen 
Größen schließen : 

• Die Berechnung von erfolgt im gesamten Phasenraum am besten durch die 
Elektronmethode. 
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Abbildung 5.1: Rekonstruktion mit unterschiedlichen Methoden : 

a) Elektronmethode, b) Sigma (E), c) Doppelwinkelmethode (DA), d) Jacquet-Blondel 
(JE) 
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Abbildung 5.2: ys Rekonstruktion bei > 100 GeV mit unterschiedlichen 
Methoden : a) Elektronmethode, b) Sigma (E), c) Doppelwinkelmethode (DA), 
d) Jacquet-Blondel (JB) 
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Abbildung 5.3: ys Rekonstruktion bei < 100 GeV mit unterschiedlichen 
Methoden : a) Elektronmethode, b) Sigma (E), c) Doppelwinkelmethode (DA), 
d) Jacquet-Blondel (JE) 
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Abbildung 5.4: a) ISO-W^^ Linien, b) W'^ Verteilung und c) ^^vertex Verteilung. 

Die Verteilungen b) und c) sind auf gleiche Flächen normiert. 

• Die Berechnung von i/b erfolgt für yß > 0.2 mit der Elektronmethode, ansons- 
ten mit der Doppelwinkelmethode. 

• Ein Schnitt auf < 0.7 ist aufgrund von Problemen bei der Rekonstruktion 
und aufgrund des Photoproduktionsuntergrundes sinnvoll. 



5.2.2 Kinematische Schnitte 

Ein weiterer Schnitt ergibt sich aus den Strahlungskorrekturen. Diese sind bei 
kleinem t/ß bzw. bei kleinem W"^ besonders groß. 

Nach Gleichung ß.lO| gilt = ys — Q"^ . Dies ist in Abbildung |5.4| a verdeutlicht. 
Wir erkennen, daß ein Schnitt von W"^ > 5000 GeV^ ungefähr einem yß > 0.05 
Schnitt entspricht. Dies gilt nahezu unabhängig von Q^. Die Verteilung für alle 



MEPS Ereignisse ist in Abbildung p.4| 3 gezeigt. 

Wie wir in Kapitel ^]2| gesehen haben, sind unsere Jetraten unabhängig von den 
Effizienzen und der Luminosität. Dies gilt natürlich nur dann, wenn wir keinen 
Effizienzunterschied in den einzelnen Ereignissen berücksichtigen müssen. Es ist 
sinnvoll den Bereich des primären Vertex einzuschränken, da so bei jedem Ereignis 
die Detektorkomponenten unter den gleichen Winkeln erscheinen und wir damit für 
alle Ereignisse gleiche Bedingungen erhalten. Die einlaufenden Teilchenstrahlen sind 
in der x — Ebene stark gebündelt. Daher ist nur ein Schnitt in Strahlrichtung 
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sinnvoll. Die ^Vertex Verteilung für die beiden Lepto-Monte-Carlo sind in Abbildung 
|5.4| c gezeigt. Der nominale Wecliselwirkungspunkt liegt bei z = 4 cm. Einen Schnitt 
von —21 cm < zvertex < 29 cm erfüllen 2% der Ereignisse nicht. 
Wir führen somit folgende kinematischen Schnitte ein : 

• W^ = yE^ s-Ql> 5000 GeV^ 

• —21 cm < -^Vertex < 29 cm 



5.2.3 Methoden zur Identifizierung des gestreuten Elektrons 

Wir haben gesehen, wie wichtig es für die kinematischen Berechnungen bei Verwendung 
der Elektronmethode ist, das gestreute Elektron richtig im Detektor zu identifizieren. 

Da das Elektron in zwei verschiedene Detektorkomponenten gestreut werden 
kann, benutzen wir für jede eine eigene Routine, sogenannte Elektronfinder. Im 
Flüssig- Argon-Kalorimeter wird die Routine QFSELH, im rückwärtigen elektromag- 
netischen Kalorimeter die Routine QFSELM verwendet. Beide sind in der offiziellen 
Hl Funktionsbibliothek HlPHAN fHl 95b[ | enthalten und wurden mit Ausnahme des 



Hb Schnittes mit den voreingestellten Steuerkarten benutzt. Der i/b Schnitt wurde 

nicht in der QFSELH Routine durchgeführt, sondern erst nach der Suche der Elektronkandidaten. 

Wurden keine Elektronkandidaten gefunden, so wird das Ereignis auf Detektorniveau 
verworfen. Ergibt die Suche in nur einer Routine einen Kandidaten, wird dieser für 
die weitere Betrachtung verwendet. Finden jedoch beide Routinen einen Kandidaten, 
so wird der Kandidat im BEMC als gestreutes Elektron identifiziert, wenn dessen 
Energie größer als 10 GeV ist. Ansonsten wird für die nachfolgenden Schnitte der 
LAr Kandidat verwendet. Dadurch verwerfen wir Elektronkandidaten, die durch 
Hadronen im BEMC vorgetäuscht werden. 

Nachdem somit ein Kandidat festgelegt ist, werden in Abhängigkeit von der 
registrierenden Detektorkomponente die in den folgenden Abschnitten erläuterten 
Bedingungen überprüft. 

Auf Partonniveau werden die Schnitte anhand des Bereiches ausgewählt, 
d.h. für < 100 GeV werden die BEMC Schnitte, ansonsten die LAr Schnitte 
durchgeführt. 



Elektronkandidaten im BEMC 

Für Elektronkandidaten im BEMC werden einige zusätzliche Schnitte eingeführt. 

Im Bereich kleiner Viererimpulsüberträge ist eine Kontamination des Datensatzes 
durch Photoproduktionsereignisse zu erwarten. Diese Ereignisse verwerfen wir durch 
ein minimales Q von 10 GeV^ und eine minimale Energie des Kandidaten von 
14 GeV. 

Ebenso sollte durch eine Winkelüberprüfung sichergestellt werden, daß das gestreute 
Elektron auch wirklich in den Bereich des BEMC fällt. 
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Da im BEMC in der Mitte aufgrund des Strahlrohres ein quadratischer Bereich 
unbestückt ist, müssen wir Kandidaten aus diesen Bereich durch einen Schnitt auf 
die X und y Koordinate des Elektronclusterschwerpunktes {center of gravity) aus 
unserem Datensatz entfernen, damit das gestreute Elektron vollständig im Detektor 
enthalten ist. 

Ein weiterer Schnitt ergibt sich aus der Tatsache, daß außer dem Kalorimeter 
noch eine Rückwärtige Proportional Kammer (BPC, hackward proportional Chamber) 
vorhanden ist. Hier fordern wir einen Treffer im Umkreis von 15 cm um unseren 
Kandidaten. 

Zusammengefaßt ergibt sich somit : 

• 10 GeV^ < < 100 GeV^ 

• 160° <9e< 172.5° 

• Ee>U GeV 

• \xcg\ > 14 cm oder \ycG\ > 14 cm 

• dßpc < 15 cm 

Elektronkandidaten im LAr 

Auch hier machen wir Schnitte in und im Elektronwinkel ög. 

Da unter kleineren Winkeln auch vermehrt Hadronen zu finden sind, müssen wir 
hier die Energieverteilung in den beiden unterschiedlichen Detektorteilen getrennt 
betrachten. Im hadronischen Kalorimeter sollte nur wenig Energie vorhanden sein, da 
Elektronen im vorgelagerten elektromagnetischen Teil ihre Energie komplett abgeben. 
Im elektromagnetischen Teil selber sollte eine relativ kleine, im Gegensatz zu hadronischen 
Jets kompaktere Ladungsverteilung zu erkennen sein. 

Zusätzlich wird noch ein Schnitt in der Größe E — für das gesamte Ereignis 
durchgeführt. Ein Wert in Höhe der zweifachen Strahlelektronenergie stellt sicher, 
daß der hadronische Endzustand vollständig erkannt wurde. Dies ist nicht der Fall, 
wenn der Elektronkandidat durch ein Hadron vorgetäuscht wurde und das wirkliche 
gestreute Elektron undetektiert bleibt. Näheres zur Rekonstruktion des hadronischen 
Endzustandes wird im nächsten Abschnitt erläutert. 

Zusammengefaßt ergeben sich für LAr Kalorimeter Elektronen die Schnitte : 

• > 100 GeV^ 

• 10° <9e< 148° 

• Eemc < 1-2 GeV in einem Konus 15 cm < r < 30 cm im elektromagnetischen 
Kalorimeter um den Elektronkandidaten. 
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• Eh AD < 0.5 GeV in einem Konus r < 30 cm im hadronischen Kalorimeter um 
den Elektronkandidaten. 

• 38 GeV <J:E-p,<70 GeV 

5.2.4 Rekonstruktionsqualität des hadronischen Endzustandes 

Wir müssen zunächst die Partikel bestimmen, die wir in unserem Jetalgorithmus 
zu Jets zusammenfassen wollen. Dies sind auf Partonniveau die vom Generator 
erzeugten Partonen und auf Detektorniveau die vom Rekonstruktionsprogramm Hl- 



REC|H1 95c| bestimmten hadronischen Cluster mit dem für die Algorithmen benötigten 



Pseudoteilchen (siehe hierzu Kapitel [4.1.2| ). 



Die richtige Partikelidentifizierung können wir überprüfen, indem wir uns die 
Energie- und Impulsverteilungen ansehen. 

Für die Energie erwarten wir die Summe der Energien der Strahlteilchen 

Etot = 820.0 GeV + 27.5 GeV = 847.5 GeV. (5.2) 

Die Verteilungen für die Ereignisse, die mit den drei Monte-Carlo-Generatoren erzeugt 



wurden, sehen wir in Abbildung |5.5| a,e und i. 

Wir erkennen, daß auf Partonniveau alle Ereignisse bei der erwarteten Energie 
auftreten. Bei den Clustern ist die Verteilung wesentlich breiter, der Peak ist jedoch 
auch hier — im Rahmen der Ungenauigkeit durch die Verluste — an der erwarteten 
Stelle. Die Form der Verteilung ist zum Teil durch die globale Kalibration der 
Kalorimeter zu erklären. Untersuchungen haben Unterschiede in der Empfindlichkeit 
der Detektoroktanten und der Detektorscheiben von einigen Prozent ergeben ||KM95| 



Der größte Teil wird aber durch Teilchen verursacht, die durch totes Material aufgehalten 
worden sind oder den Detektor unbeobachtet verlassen haben. Durch die Einführung 
des Pseudoteilchens ist die Aussagekraft der hadronischen Gesamtenergieverteilung 
jedoch eingeschränkt. 

Die Qualität der Rekonstruktion des Protonrestes durch das Pseudoteilchen läßt 
sich an der Verteilung der zur Strahlachse transversalen Impulse erkennen. Da beide 
einlaufenden Strahlteilchen keine transversalen Impulse besitzen, erwarten wir auch 
für die Impulssumme über alle Endzustandspartikel eine Ausrichtung entlang der 
Strahlachse. Die Verteilungen im Partonniveau bestätigen dies (Abbildung p75| b, c, f, 
g, j, k). Auf Detektorniveau ist ebenfalls ein Maximum bei Null zu sehen, jedoch ist 
auch hier die Verteilung breiter. Die Werte liegen im Rahmen der Meßungenauigkeit. 
Da das Pseudoteilchen keinen transversalen Impuls erhält, ist die Rekonstruktion des 
hadronischen Endzustandes somit gelungen und der Pseudoteilchenansatz gerechtfertigt. 

Die Impulssumme in Richtung der Strahlachse stimmt auf Detektorniveau durch 
die Definition des Pseudoteilchens mit dem Erwartungswert überein. Dies ist hier 
deshalb nicht gezeigt. 

In Abbildung [5.5|d, h und 1 ist die Verteilung der Summe von E — p^ gezeigt. Aus 



den obengenannten folgt, daß wir hier eine Kopie der Energieverteilung sehen, jetzt 
jedoch um den Erwartungswert 2Ee. 
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Abbildung 5.5: Verteilung der Energie- und Impulssummen. Dargestellt für 
Lepto (a-d), Ariadne (e-h) und Herwig (i-l) mit Detektor- (durchgezogen) und 
Partonniveau (gestrichelt). 
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Abbildung 5.6: Verteilung der Energie- und Impulssummen. Dargestellt für 
Daten (durchgezogen) und Lepto- Monte- Carlo (gestrichelt), nur Detektorniveau. 



In Abbildung |5^ sind die gleichen Verteilungen für die Daten gezeigt. Die 
Übereinstimmung zwischen den Daten und dem Monte-Carlo-Detektorniveau ist gut, 
die Datenverteilung der Energie ist jedoch etwas breiter. 

Aufgrund der Verteilung ist es sinnvoll Schnitte einzuführen, die Ereignisse mit 
zu großen Abweichungen verwerfen. Wir werden deshalb in dieser Arbeit folgenden 
Grenzen verwenden : 

• Die Energiesumme muß auf Partonniveau besser als 1%, auf Detektorniveau 
besser als 4% mit dem Erwartungswert übereinstimmen. 

• Die transversalen Impulse px und Py müssen jeweils kleiner als 1 GeV auf 
Parton- und 5 GeV auf Detektorniveau sein. 



5.2.5 Q2 Einteilung 



Für die Messung der starken Kopplungskonstanten in verschiedenen Bereichen, 
müssen wir eine Einteilung in sogenannte Eins vornehmen. 

Bei der ersten Messung der starken Kopplungskonstante mit den Daten von 
1993 wurde eine Fünfteilung in drei niedrige und zwei hohe Bins vorgenommen 
Nis94b|| , | [Nis94a|| , ||Nis95a|| . Mit der höheren Statistik aus dem Jahr 1994 ist eine 



Einteilung in acht Bins, fünf niedrige und drei hohe, möglich ||Nis95b|| , ||Eis95|| . 

Die Einteilung wurde so gewählt, daß die Anzahl der 2+1 Ereignisse in den Bins 
bei den Daten ungefähr gleich groß ist, gleichzeitig der Bereich aber nicht zu weit 
wird. Diese Wahl ist sinnvoll, da der statistische Fehler der Jetrate fast ausschließlich 
durch den Fehler der 2+1 Jetanzahl bestimmt wird. Ich behalte diese Einteilung 
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Bin 


Bereich 


mittlerer Wert 


Anzahl Ereignisse 


Anzahl 2+1 Ereignisse 






GeV^ 






GeV^ 




nach Stdschnitten 


nach Stdschnitten 


1 


10 — 14 


11.9 


14810 


14951 


744 


854 


2 


14 — 1 


8 


15.9 


10788 


10893 


582 


651 


3 


1 


8 — 25 


21.4 


10239 


10329 


650 


674 


4 


25 — 40 


32.1 


9451 


9578 


746 


701 


5 


40 — 100 


66.6 


6604 


6547 


628 


604 


6 


100 — 300 


186.6 


7917 


7749 


1142 


961 


7 


300 — 700 


479.1 


5047 


4660 


1049 


896 


8 


700 — 4000 


2011.7 


2236 


1848 


521 


396 



Tabelle 5.2: Einteilung in Bins. Neben den Grenzen ist das mittlere 
im jeweiligen Bereich eingetragen. Die Einteilung wurde so gewählt, daß die 2+1 
Jetanzahlen in allen Bins ungefähr gleich groß ist. Mehr hierzu im Text. Bei den 
Anzahlen bezieht sich die erste Zahl jeweils auf das Parton-, die zweite auf das 
Detektorniveau. 



deshalb bei. Die genauen Angaben über die einzelnen Bins sind in Tabelle |5.2| 
aufgeführt. 

Die Anzahl der Ereignisse mit 3 oder mehr harten Jets liegt auf beiden Niveaus 
summiert über alle Bins unter 3%0. Der Einfluß dieser Ereignisse kann somit vernachlässigt 
werden. 



5.3 Fehlerbestimmung 

Der statistische Fehler der Jetanzahlen ist gegeben durch die Wurzel 



A^. 

Die Jetraten berechnen sich ähnlich Effizienzen durch 

A^2+l 



R 



n+l 



N 



(5.3) 



(5.4) 



tot 



Wir müssen daher einen Effizienzenfehler benutzen, da die beiden Anzahlen korreliert 
sind : 



Ae 



'e(l-e) 



N 



(5.5) 



^Weniger als 200 Ereignisse mit mehr als zwei harten Jets pro Niveau. 
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oder in Jetanzahlen geschrieben 

Ai?2+1 



N2+1 


(1 


R2+1 (1 — R2+1) 



(5.6) 



tot 



Den Fehler für den Korrekturfaktor 



j^Parton 

~ ^Detektor {^■'^) 

hingegen berechnen wir durch lineare Fehlerfortpflanzung. Der Fehler für einen 
Quotienten zweier Größen ist dann 

R - ^ (5.8) 



m = ^^■^- + ^^-^' (5.9) 



Auch hier haben wir zwar eine Korrelation, da die gleichen Ereignisse verwendet 
wurden, um die Raten auf den unterschiedlichen Niveaus zu bestimmen. Durch 
Benutzung verschiedener Monte-Carlos mit gleichen Einstellungen für Parton- und 
Detektorniveau wäre es mögliche diese Korrelation vollständig auszuschalten. Wir 
können sie jedoch bei den hohen Ereigniszahlen vernachlässigen. 



5.4 Motivation zur Untersuchung der Jet— Parton— 
Korrelation 

Betrachten wir die Korrelation der (n + 1) Jet Ereignisanzahlen zwischen Parton- 
und Detektorniveau, so ergeben sich Matrizen der Art 



Bin 


Parton 


1 


1 + 1 


2+1 


Clus- 


1+1 


an 


021 


ter 


2+1 


012 


Ö22 



Die Hauptdiagonalelemente an und beschreiben die auf beiden Niveaus gleich 
klassifizierten Ereignisse und die Nebendiagonalelemente a2i und ai2 die Migrationen, 
d.i. der Fall, in dem durch Fehl-Erkennung ein Ereignis beim Ubergang von Parton- 
auf Detektorniveau die Jetklasse wechselt. Diese Elemente sollten möghchst niedrige 
Werte annehmen. 



13 MOTIVATION 



53 



Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


1 


1 + 1 


2+1 


2 


1 + 1 


2+1 


3 


1 + 1 


2+1 


Clus- 


1+1 


12458 


326 


Clus- 


1 + 1 


9419 


247 


Clus- 


1+1 


8970 


291 


ter 


2+1 


434 


343 


ter 


2+1 


336 


267 


ter 


2+1 


336 


287 



Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


4 


1 + 1 


2+1 


5 


1 + 1 


2+1 


6 


1 + 1 


2+1 


Clus- 


1+1 


8108 


345 


Clus- 


1+1 


5464 


256 


Clus- 


1+1 


5934 


448 


ter 


2+1 


305 


345 


ter 


2+1 


259 


299 


ter 


2+1 


274 


586 



Bin 


Parton 


Bin 


Parton 


7 


1 + 1 


2+1 






1+1 


2+1 


Clus- 


1+1 


3324 


341 


Clus- 


1+1 


1287 


154 


ter 


2+1 


234 


605 


ter 


2+1 


99 


286 



Bin 


1 


2 


3 


4 


5 


6 


7 


8 


Rate Parton 


4.9 


5.0 


5.8 


7.6 


8.8 


14.3 


21.0 


24.1 


in % 


±0.2 


±0.3 


±0.3 


±0.3 


±0.4 


±0.5 


±0.7 


±1.1 


Rate Detektor 


5.7 


5.9 


6.3 


7.1 


8.9 


11.9 


18.6 


21.1 


in % 


±0.2 


±0.3 


±0.3 


±0.3 


±0.4 


±0.4 


±0.6 


±1.0 


Korrektur- 


0.86 


0.85 


0.93 


1.06 


0.99 


1.20 


1.13 


1.14 


faktor 


±0.07 


±0.08 


±0.08 


±0.08 


±0.09 


±0.08 


±0.07 


±0.10 


Reinheit in % 


44.1 


44.3 


46.1 


53.1 


53.6 


68.1 


72.1 


74.3 



Tabelle 5.3: Migrationen und Korrekturfaktoren bei Lepto mit 
Standardschnitten, in Bins 
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Benutzen wir nur die Standardschnitte, so ergeben sich die in Tabelle |5^ verzeichneten 
Migrationen für die MEPS Datensätze in den verschiedenen Bins. Es wurde hier 
der JADE Algorithmus verwendet^. 

Aus den Migrationen, beschrieben durch die Elemente aparton,ciustor lassen sich die 
Raten einfach berechnen 

^Detektor ^ 0^2 + 0^ ^^^^^^ 
Y,i,j O'ij 

RllT = ^^^^ (5.11) 



Neben den Raten ist in der Tabelle 5.3 auch der Korrekturfaktor 



R = ^^^±^ (5.12) 

022 + 0-12 

berechnet worden. Hier ist jedoch anzumerken, daß bei den Migrationstabellen 
nur solche Ereignisse verwendet werden konnten, die auf beiden Niveaus die Standardschnitte 
erfüllen. Im sechsten Bin haben wir z.B. insgesamt 7242 Ereignisse, obwohl in diesem 
Bin 7749 Ereignisse die Standardschnitte auf Detektorniveau erfüllen und sogar 7917 
die auf Partonniveau (siehe Tabelle p.2| ). 

Für unsere späteren Betrachtungen benutzen wir die Raten, die sich aus den 
kompletten Sätzen ergeben. Diese können sich somit leicht von den hier angegebenen 
unterscheiden. 

Bei den Raten fällt uns auf, daß diese mit anwachsendem ebenfalls steigen. 
Wir erwarten für die starke Kopplungskonstante Us jedoch genau einen umgekehrten 
Verlauf. Die Abhängigkeit der Vorfaktoren in Gleichung |4.8| ist jedoch stärker, so 
daß der Abfall der starken Kopplungskonstanten mit überkompensiert wird. 

Die Korrekturfaktoren liegen bei eins und zeigen keine großen Probleme bei der 
Jet-Parton-Korrelation an. 

Sehen wir uns die Migrationen genauer an, so fällt jedoch auf, daß die 2+1 
Jetklassen auf den beiden Niveaus nicht gut übereinstimmen. In den unteren Bins 
liegt die Reinheit der 2+1 Jetereignisse auf Detektorniveau 

rParton und Detektor 



ArFart- 

^''2+1 

/\ /"Detektor 
^^2+1 

(5.13) 



022 + O12 

bei nur 45 bis 50%. In den hohen Bins ist sie mit 65 bis 75% zwar besser, aber 
sicherlich nicht zufriedenstellend. 



•^Da der JADE Algorithmus mit dem Auflösungsparameter j/cut = 0.02 unsere Standardwahl 
ist, wird im folgenden der Algorithmus nur bei Verwendung eines anderen Algorithmuses explizit 
angegeben. 
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Zur Bekämpfung von Partonschauern wurde schon bei der Analyse der im Jahr 
1993 genommenen Daten ein Jetwinkelschnitt eingeführt 

10° < öjet < 145°. (5.14) 



Eine genauere Untersuchung dieses Schnittes wird in Kapitel p.8.1| durchgeführt. 
In Tabelle ß.4| sind die Migrationen für diesen Schnitt notiert. 



Wir erkennen hier die gleiche Systematik. Für die Abhängigkeit gilt auch 
hier das obengesagte. Der Korrekturfaktor liegt nun jedoch mit durchschnittlich 1.4 
weiter von 1 entfernt als ohne 6 Schnitt. Um dies zu erklären, vergleichen wir die 
Migrationen näher. 

Betrachten wir zuerst die niedrigen Bins. Durch den zusätzlichen Schnitt hat 
sich die Anzahl der Ereignisse, die in beiden Niveaus zwei harte Jets haben (2+1 
Jet, 022) auf ca. 40% reduziert. Die Ereignisse ai2, die auf Partonniveau einen harten 
Jet haben, auf Detektorniveau jedoch zwei, sind zu über 70% unterdrückt worden. 
Im Gegensatz dazu werden nur etwa 42% der Ereignisse, die auf Detektorniveau 
einen harten Jet, auf Partonniveau jedoch zwei haben, verworfen. Die Konsequenzen 
abgesehen vom vergrößerten Korrekturfaktor sind kleinere Raten auf beiden Niveaus 
und bei kleinen eine leicht verbesserte Reinheit unserer Detektorprobe. 

Dieser Schnitt deutet auf Unterschiede in der Parton- und der Detektor- Jetwin- 
kelverteilung hin. Hierauf werden wir in Kapitel ß.8.1| zurückkommen. 



Bei den hohen Q"^ Bins ist die Unterdrückung um 55% für 022, 51% für 012 und 
29% für 021 nicht ganz so extrem, aber trotzdem sehr unterschiedlich. Insbesondere 
ist der Effekt auf die auf beiden Niveaus gleich erkannten Ereignisse groß, was die 
höheren relativen Fehler bewirkt. 

Um die Stärke der Migrationen zu charakterisieren genügt also nicht nur der 
Korrekturfaktor, sondern die Reinheit unserer Jetklassifizierung ist ebenfalls wichtig. 

Unser Problem bei der Jetanalyse ist deutlich geworden. Wir suchen einen oder 
mehrere Schnitte, die die Migrationen, d.h. die Elemente 012 und 021 der Migrationsmatrix 
unterdrücken. Außerdem sollte die Reduktion möglichst auf beide Elemente gleich 
stark ausfallen, da sonst der Korrekturfaktor zu großQ wird. Wünschenswert wäre 
zusätzlich, daß die auf beiden Niveaus gleicherkannten Ereignisse, insbesondere die, 
die beide Male als 2+1 Jetereignisse klassifiziert wurden, nicht zu sehr reduziert 
werden. Dies würde sonst den statistischen Fehler stark erhöhen. 



5.5 Verschiedene Monte-Carlo-Generatoren 

Bevor wir nun weitere Schnitte testen, ist es sinnvoll die Abhängigkeit vom Monte- 
Carlo- Generator und dem Jetalgorithmus zu betrachten. Ersteres untersuchen wir 
in diesem Abschnitt, letzteres im nächsten. 

■^Unter einem großen Korrekturfaktor verstehen wir einen Wert, der eine große Abweichung vom 
idealen Wert 1 bedeutet. Es ist dabei nebensächHch, ob die Abweichung nach unten, d.i. ein kleiner 
Korrekturfaktorwert, oder nach oben, d.i. ein großer Korrekturfaktorwert, ausfällt. 
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Bin 
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Bin 


Parton 


Bin 


Parton 


1 


1+1 


2+1 


2 


1+1 


2+1 


3 


1+1 


2+1 


Clus- 


1+1 


6830 


219 


Clus- 


1+1 


6500 


135 


Clus- 


1+1 


6712 


182 


ter 


2+1 


136 


147 


ter 


2+1 


71 


97 


ter 


2+1 


95 


109 



Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


4 


1+1 


2+1 


5 


1+1 


2+1 


6 


1+1 


2+1 


Clus- 


1+1 


6707 


183 


Clus- 


1+1 


5333 


138 


Clus- 


1+1 


5197 


323 


ter 


2+1 


91 


122 


ter 


2+1 


67 


114 


ter 


2+1 


141 


274 



Bin 


Parton 


Bin 


Parton 


7 


1 + 1 


2+1 




i 


1+1 


2+1 


Clus- 


1+1 


3804 


246 


Clus- 


1+1 


1506 


105 


ter 


2+1 


109 


249 


ter 


2+1 


48 


141 



Bin 


1 


2 


3 


4 


5 


6 


7 


8 


Rate Parton 


5.0 


3.4 


4.1 


4.3 


4.5 


10.1 


11.2 


13.7 


in% 


±0.3 


±0.3 


±0.3 


±0.3 


±0.3 


±0.4 


±0.5 


±0.9 


Rate Detektor 


3.9 


2.5 


2.9 


3.0 


3.2 


7.0 


8.1 


10.5 


in 'X 


±0.3 


±0.2 


±0.2 


±0.3 


±0.3 


±0.1 


±0.5 


±0.8 


Korrektur- 


1.29 


1.38 


1.43 


1.43 


1.39 


1.44 


1.38 


1.30 


faktor 


±0.15 


±0.20 


±0.19 


±0.18 


±0.19 


±0.13 


±0.13 


±0.17 


Reinheit in % 


51.9 


57.7 


53.4 


57.3 


63.0 


66.0 


69.6 


74.6 



Tabelle 5.4: Migrationen und Korrekturfaktoren bei Lepto mit 
Standardschnitten und zusätzlichem 10° < Oj^t < 145° Schnitt, in Bins 
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Wir haben neben unserem Lepto-MEPS-Monte-Carlo noch Simulationen von Ereignissen, 
die mit Ariadne bzw. mit Herwig erzeugt wurden. Die genauen Angaben zu diesen 
Datensätzen sind in Tabelle 5J. aufgeführt. 

Die Migrationen sind in Tabelle |5.5| für Ariadne und in Tabelle |5.6| für Herwig mit 
der gleichen Einteilung (siehe Tabelle |5^) unter Verwendung der Standardschnitte 
aufgeführt. 

Vergleichen wir die drei Tabellen, stellen wir fest, daß die Abhängigkeit der 
Raten bei den drei Monte- Carlos unterschiedlich ist. Bei dem mit Herwig generierten 
Datensatz ist sie am größten und beim Ariadne Datensatz am geringsten. 

Damit sind natürlich auch die Raten unterschiedlich. Die Abweichungen sind bei 
den niedrigen Eins jedoch etwas größer. 

Bei den Korrekturfaktoren wird die Abhängigkeit und die Abhängigkeit vom 
benutzten Monte-Carlo besonders deutlich. Im ersten Bin reichen die Werte von 
0.58 bis 1.25, im zweiten immerhin noch von 0.85 bis 1.31. In den letzten beiden 
Bins liegen die Werte zwischen 1.01 und 1.14 und stimmen im Rahmen der Fehler 
überein. Dies ist bereits ein deutlicher Hinweis darauf, daß die Probleme bei den 
Jet-Parton-Korrelationen in den niedrigen Bins wesentlich größer sind. 

Zusammenfassend können wir also sagen, daß unsere Berechnungen abhängig 
vom benutzen Monte-Carlo-Generator sind. Wir werden deshalb in Kapitel |5.8| 
verschiedene Verteilungen untersuchen, um Probleme bei unseren simulierten Ereignissen 
im Vergleich mit den Daten zu finden und die entsprechenden Phasenraumbereiche 
dann ausschließen zu können. 



5.6 Verschiedene Jetalgorithmen 

Wie bereits im letzten Kapitel angesprochen, müssen wir auch die Abhängigkeit 
der Werte vom verwendeten Jetalgorithmus näher betrachten. Dazu wurden auf 
den beiden Niveaus der Lepto Simulationen der JADE und der Cone Algorithmus 
angewandt. Aus Kapitel [4.1.2| wissen wir, daß bei der Berechnung der invarianten 
Masse im JADE Schema die Teilchenmassen vernachlässigt werden. Dies bedeutet 
jedoch, daß der Jetalgorithmus insgesamt nicht lorentz-invariant ist. Um diesen 
Effekt näher zu betrachten, wurde der JADE Algorithmus mit JADE Schema in drei 
unterschiedlichen Bezugssystemen angewendet. Dazu wurde neben dem Laborsystem, 
das unsere Standardeinstellung ist, das hadronische Schwerpunktsystem und das 
Breitsystem (siehe Kapitel 1.1.3 ) verwendet. Aufgrund des hohen Energieunterschiedes 



von Strahlproton und -elektron können jedoch Lorentz Transformationen mit extremen 
Werten in den Rotationsmatrizen auftreten Die Lorentz- Transformationen sind 
deshalb nicht unproblematisch und spezielle Tests mußten eingebaut werden, um 
nichtphysikalische Werte zu verhindern. Die einzelnen Partikel sind vor der Berechnung 
der Jets in das entsprechende System transformiert worden und nach der Anwendung 



'In diesen Fällen ist /3 = - « 1. 
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Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


1 


1+1 


2+1 


2 


1+1 


2+1 


3 


1+1 


2+1 


Clus- 


1+1 


5728 


454 


Clus- 


1+1 


4497 


400 


Clus- 


1+1 


4363 


391 


ter 


2+1 


305 


296 


ter 


2+1 


243 


262 


ter 


2+1 


227 


277 



Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


4 


1+1 


2+1 


5 


1+1 


2+1 


6 


1+1 


2+1 


Clus- 


1+1 


3997 


117 


Clus- 


1+1 


2552 


301 


Clus- 


1+1 


2875 


258 


ter 


2+1 


203 


282 


ter 


2+1 


154 


229 


ter 


2+1 


137 


263 



Bin 


Parton 


Bin 


Parton 


7 


1+1 


2+1 


8 


1+1 


2+1 


Clus- 


1+1 


1663 


169 


Clus- 


1+1 


673 


77 


ter 


2+1 


138 


214 


ter 


2+1 


75 


86 



Bin 


1 


2 


3 


4 


5 


6 


7 


8 


Rate Parton 


11.1 


12.3 


12.7 


14.3 


16.4 


14.7 


17.5 


17.9 


in % 


±0.4 


±0.5 


±0.5 


±0.5 


±0.7 


±0.6 


±0.9 


±1.3 


Rate Detektor 


8.9 


9.3 


9.6 


9.9 


11.8 


11.3 


16.1 


17.7 


in 'X 


±0. 1 


±0.1 


±0.5 


±0.5 


±0.6 


±0.6 


±0.8 


±1.3 


Korrekt ur- 


1.25 


1.31 


1.33 


1.44 


1.38 


1.30 


1.09 


1.01 


faktor 


±0.10 


±0.11 


±0.11 


±0.12 


±0.13 


±0.12 


±0.11 


±0.15 


Reinheit in % 


49.3 


51.9 


55.0 


58.1 


59.8 


65.8 


60.8 


53.4 



Tabelle 5.5: Migrationen und Korrekturfaktoren bei Ariadne mit 
Standardschnitten, in Bins 
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Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


1 


1+1 


2+1 


2 


1+1 


2+1 


3 


1+1 


2+1 


Clus- 


1+1 


4589 


56 


Clus- 


1+1 


3426 


94 


Clus- 


1+1 


3167 


108 


ter 


2+1 


129 


43 


ter 


2+1 


116 


74 


ter 


2+1 


122 


78 



Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


4 


1+1 


2+1 


5 


1+1 


2+1 


6 


1+1 


2+1 


Clus- 


1+1 


2950 


90 


Clus- 


1+1 


1879 


91 


Clus- 


1+1 


6223 


304 


ter 


2+1 


135 


109 


ter 


2+1 


94 


93 


ter 


2+1 


217 


491 



Bin 


Parton 


Bin 


Parton 


7 


1+1 


2+1 






1 + 1 


2+1 


Clus- 


1+1 


3570 


266 


Clus- 


1+1 


1237 


116 


ter 


2+1 


189 


522 


ter 


2+1 


86 


288 



Bin 


1 


2 


3 


4 


5 


6 


7 


8 


Rate Parton 


2.1 


4.5 


5.4 


6.1 


8.5 


11.0 


17.3 


23.4 


in % 


±0.3 


±0.4 


±0.4 


±0.5 


±0.7 


±0.4 


±0.6 


±1.1 


Rate Detektor 


3.6 


5.1 


5.8 


7.4 


8.7 


9.8 


15.6 


21.7 


in 'X 


±0.3 


±0. 1 


±0.1 


±0.5 


±0.7 


±0.1 


±0.6 


±1.0 


Korrekt ur- 


0.58 


0.88 


0.93 


0.82 


0.98 


1.12 


1.11 


1.08 


faktor 


±0.11 


±0.13 


±0.14 


±0.11 


±0.14 


±0.08 


±0.08 


±0.10 


Reinheit in % 


25.0 


38.9 


39.0 


44.7 


49.7 


69.4 


73.4 


77.0 



Tabelle 5.6: Migrationen und Korrekturfaktoren bei Herwig mit 
Standardschnitten, in Bins 
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Abkürzung 


Algorithmus 


Lorentzsystem 


Schnittparameter 


JL 


JADE 


Labor 


0.02 


JC 


JADE 


hadr . S chwerpunkt 


0.02 


JB 


JADE 


Breit 


0.02 


3L 


JADE 


Labor 


0.03 


GL 


ConeQ 


Labor 


0.7 


EL 


E Schema 


Labor 


0.02 



"Mit den Parametern -Et, min = 1-5 GeV und i?t,sum,min = 7.0 GeV 



Tabelle 5.7: Aufstellung der in den Abbildungen |5.7| bis |5.11| benutzten 
Jetalgorithmen. Genaueres zu den einzelnen Algorithmen ist in Kapitel [^. j| zu 
finden. 



des JADE Algorithmus wurden die Jets ins Laborsystem zurücktransformiert, um 
einen Vergleich der drei Versionen zu ermöglichen. 

Als weitere Variante wurde der JADE Algorithmus mit einem veränderten Schnittwert 
von Heut = 0.03 benutzt, um die Größenordnung der Abhängigkeit vom Auflösungs- 
parameter zu bestimmen (siehe Kapitel [4.2|) . 

Zum Schluß wurde auch der JADE- Algorithmus unter Verwendung der Massenberechnung 
und Rekombination durch das E Schema eingesetzt, das den Algorithmus lorentz- 
invariant macht. Frühere Untersuchungen der JADE Varianten haben jedoch schon 
gezeigt, daß gerade dieser Algorithmus große Korrekturfaktoren bewirkt ( ||Nis94b| 
Kapitel 10, Abbildung 10.19 und 10.20). 



Eine Aufstellung der verwendeten Algorithmen ist in Tabelle |5.7| wiedergegeben. 

In den Abbildungen |5.7| bis |5.11| sind die entsprechenden Werte graphisch für 
verschiedene Bins dargestellt. Die Ereignisse aus den ersten vier Bins wurden 
dabei zu einer Abbildung zusammengefaßt. Zum einfacheren Vergleich sind die 
Graphen in einer festen Skalierung gezeigt. Die Werte für alle acht Bins sind in 



den Tabellen q]8| bis ß.lO| zusammengefaßt. 

Die Abbildungen ermöglichen sowohl den Vergleich zwischen Daten und Monte- 
Garlo-Detektorniveau als auch die Bewertung der Korrekturfaktoren R und den 
Vergleich der korrigierten Raten 



Rl 



korr 
2+1 



R ■ R. 



Daten 
■2+1 



(5.15) 



pParton 

-"-2+1 pDaten 

pDetektor ' ""-2+1 
-"■2+1 



Zum letzten Punkt ist anzumerken, daß wir nur für den JADE Algorithmus mit 
Schnittwert Ucnt = 0.02 in den verschiedenen Systemen die gleiche korrigierte Rate 
erwarten. Bei Verwendung eines anderen Schnittparameters oder anderen Algorithmus 
verändern sich die von der Theorie vorhergesagten Koeffizienten zur Berechnung 
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JL • 

JC • • 

JB ■ 

3L • • ^ • 

CL • • 

EL 

0.2 0.2 1 2 0.2 

a) Daten b) MC Detektor c) Korrekt jrfakior d) korrigierte Rote 



Abbildung 5.7: Vergleich verschiedener Jetalgorithmen, Bin 1-4. Dargestellt 
sind auf der Ordinate die Jetalgorithmen JADE im Laborsystem (JL), JADE im 
hadronischen Schwerpunktsystem (JC), JADE im Breitsystem (JB), JADE mit 
Vcnt = 0.03 (3L), Cone Algorithmus (CL) und der JADE Algorithmus mit E-Schema 
(EL) (siehe auch Tabelle In den vier Craphen sind die Raten R2+1 für die 

Daten (a) und das Monte- Carlo-Detektorniveau (b) gezeigt. Der Korrekturfaktor ist 
in c dargestellt und die sich ergebende korrigierte Rate in d. 
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0.2 0.2 1 2 0.2 

a) Daten b) MC Detektor c) Korrekturfaktor d) korrigierte Rate 



Abbildung 5.8: Vergleich verschiedener Jetalgorithmen, Bin 5. Für nähere 



Erläuterungen siehe Abbildung \5. ?| und Tabelle 5.1 



JL * • — ^ 

JC * • - ^ 

JB * • ^ 

3L • • — 

GL • • - ^ 

EL 

0.2 0.2 1 2 0.2 

o) Daten b) MC Detektor c) Korrekturfaktor d) korrigierte Rote 

Abbildung 5.9: Vergleich verschiedener Jetalgorithmen, Bin 6. Für nähere 
Erläuterungen siehe Abbildung \5. ?| und Tabelle |5. T\ 
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0.2 0.2 1 2 0.2 

a) Daten b) MC Detektor c) Korrekturfaktor d) korrigierte Rate 



Abbildung 5.10: Vergleich verschiedener Jetalgorithmen, Bin 7. Für nähere 



Erläuterungen siehe Abbildung \5. ?| und Tabelle 5.1 





0.2 0.2 1 2 0.2 

o) Daten b) MC Detektor c) Korrekturfaktor d) korrigierte Rote 



Abbildung 5.11: Vergleich verschiedener Jetalgorithmen, Bin 8. Für nähere 



Erläuterungen siehe Abbildung \5. ?| und Tabelle W7i 
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1 


2 


3 


4 


5 


6 


7 


8 


Datenrate -R2+i,Daten in % 


JL 


6.9 


8.2 


8.8 


10.2 


13.6 


14.2 


21.2 


25.4 




±0.7 


±0.8 


±0.8 


±0.8 


±0.9 


±1.2 


±1.6 


±2.3 


JC 


25.4 


26.9 


27.8 


28.4 


32.3 


26.3 


37.6 


42.1 




±0.9 


±0.9 


±0.9 


±0.9 


±1.0 


±1.4 


±1.7 


±2.5 


JB 


4.7 


6.5 


7.1 


9.6 


14.0 


15.5 


27.4 


37.5 




±0.7 


±0.7 


±0.7 


±0.8 


±0.9 


±1.2 


±1.6 


±2.4 


3L 


2.5 


3.3 


3.4 


4.5 


6.8 


7.4 


12.8 


17.3 




±0.6 


±0.7 


±0.7 


±0.7 


±0.8 


±1.0 


±1.4 


±2.0 


CL 


0.3 


0.4 


0.4 


0.6 


0.8 


3.3 


6.3 


8.4 




±0.5 


±0.6 


±0.6 


±0.6 


±0.6 


±0.9 


±1.1 


±1.6 


EL 


67.6 


69.8 


69.8 


71.3 


73.1 


70.3 


71.3 


68.6 




±0.9 


±0.9 


±0.9 


±0.9 


±1.0 


±1.4 


±1.7 


±2.4 



Tabelle 5.8: Zusammenfassung des Vergleichs verschiedener Jetalgorithmen 

der starken Kopplungskonstante (siehe Kapitel ^.2| ). Da das von uns verwendete 
Programm zur Berechnung der theoretischen Werte PRO JET nur den JADE Algorithmus 
verwendet , ist auf einen Vergleich der Partonraten mit den theoretischen Werten hier 
verzichtet worden. 

Vergleichen wir zuerst die unterschiedlichen Systeme. Uns fällt auf, daß sich bei 
der Datenrate große Unterschiede ergeben. Obwohl bei allen Versionen die Datenrate 
mit steigt, ist die Stärke der Abhängigkeit unterschiedlich. Vergleichen wir dies 
mit der Rate auf Monte-Carlo-Detektorniveau, so ergibt sich ein ähnliches Bild. Die 
Raten selber sind jedoch um etwa 2% kleiner als bei den Daten. Auf Partonniveau 
sieht das Bild schon wesentlich besser aus. Die Werte sind nun wesentlich enger 
zusammen und die Abhängigkeit ist ebenfalls sehr ähnlich. Aus den letzten 
beiden Punkten ergeben sich zwangsläufig sehr unterschiedliche Korrekturfaktoren. 
Während die Korrekturfaktoren für das Labor- und das Breitsystem zwischen 0.88 
und 1.18 im Rahmen bleiben, sind sie für den JADE Algorithmus im hadronischen 
Schwerpunktsystem mit Werten zwischen 0.43 und 0.82 unakzeptabel. 

Betrachten wir nun den JADE Algorithmus mit dem größeren Schnittwert, so 
ist die Rate auf allen drei gemessenen Niveaus im Einklang mit unserer Erwartung 
kleiner als bei der Standardeinstellung. Die Reduktion beträgt in allen Fällen ungefähr 
50%, wobei das Verhältnis in den niedrigeren Bins auf Detektorlevel eher kleiner, 
auf Partonlevel etwas höher liegt. Für den Korrekturfaktor ergibt sich ein recht 
gleichmäßiger Verlauf mit Werten zwischen 1.14 und 1.30. Für die korrigierte Rate 
gilt somit ebenfalls eine Reduktion um die Hälfte gegenüber den Standardeinstellungen. 

Der Cone Algorithmus hat sehr niedrige Raten. Die Beschreibung der Daten 
durch das Monte-Carlo-Detektorniveau ist recht gut, aber es zeigen sich große Abweichungen 
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1 


2 


3 


4 


5 


6 


7 


8 


Monte-Carlo-Detektorrate -R2+1, Cluster in % 


JL 


5.7 


6.0 


6.5 


7.3 


9.2 


12.4 


19.2 


21.5 




±0.7 


±0.7 


±0.7 


±0.8 


±0.9 


±0.9 


±1.1 


±1.5 


JC 


13.2 


14.5 


15.5 


17.0 


20.4 


20.3 


30.2 


33.0 




±0.8 


±0.8 


±0.9 


±0.9 


±1.0 


±1.0 


±1.2 


±1.6 


JB 


3.9 


4.5 


5.4 


7.0 


9.8 


13.6 


25.1 


31.0 




±0.7 


±0.7 


±0.7 


±0.8 


±0.9 


±0.9 


±1.1 


±1.6 


3L 


2.0 


2.2 


2.5 


3.2 


4.4 


6.4 


10.9 


13.3 




±0.6 


±0.6 


±0.7 


±0.7 


±0.8 


±0.8 


±1.0 


±1.3 


CL 


0.7 


0.4 


0.6 


0.9 


0.8 


3.6 


5.5 


10.2 




±0.6 


±0.6 


±0.6 


±0.6 


±0.6 


±0.7 


±0.8 


±1.2 


EL 


60.8 


64.3 


63.8 


64.8 


66.5 


63.5 


68.8 


65.7 




±0.9 


±1.0 


±1.0 


±1.0 


±1.1 


±1.0 


±1.2 


±1.6 


Monte-Carlo-Partonrate -R2+i,Parton in % 


JL 


5.2 


5.4 


6.6 


8.1 


9.6 


14.6 


21.2 


23.6 




±0.7 


±0.7 


±0.7 


±0.8 


±0.8 


±0.9 


±1.1 


±1.4 


JC 


5.7 


6.5 


7.4 


9.2 


11.9 


15.9 


24.8 


25.7 




±0.7 


±0.7 


±0.7 


±0.8 


±0.9 


±0.9 


±1.1 


±1.4 


JB 


3.9 


4.5 


5.7 


7.1 


10.1 


15.0 


24.8 


27.3 




±0.7 


±0.7 


±0.7 


±0.8 


±0.9 


±0.9 


±1.1 


±1.4 


3L 


2.3 


2.7 


3.1 


4.1 


5.2 


7.4 


12.7 


15.2 




±0.6 


±0.6 


±0.7 


±0.7 


±0.8 


±0.8 


±1.0 


±1.2 


CL 


1.9 


1.4 


1.6 


2.1 


1.8 


8.7 


12.2 


16.7 




±0.6 


±0.6 


±0.6 


±0.6 


±0.7 


±0.8 


±0.9 


±1.3 


EL 


13.1 


13.8 


15.7 


19.3 


23.0 


33.4 


43.7 


45.1 




±0.8 


±0.8 


±0.8 


±0.9 


±1.0 


±1.0 


±1.2 


±1.5 



Tabelle 5.9: Zusammenfassung des Vergleichs verschiedener Jetalgorith 
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1 


2 


3 


4 


5 


6 


7 


8 


Korrekt urfaktc )r R 


JL 


0.91 


0.90 


1.01 


1.11 


1.04 


1.18 


1.10 


1.10 






±0.07 


±0.08 


±0.08 


±0.08 


±0.07 


±0 07 

—1— W . W 1 


±0 09 

—1— W . W iJ 


JC 


0.43 


0.45 


0.48 


0.54 


0.58 


0.78 


0.82 


0.78 




+0 03 


±0.03 


±0.03 


±0.03 


±0.04 


±0.04 


±0 04 

_l_ W . Wt: 


±0 06 

_l_ w . w w 


JB 


0.99 


1.00 


1.05 


1.02 


1.03 


1.11 


0.99 


0.88 




+0 08 

_i_w . wo 


±0.09 


±0.09 


±0.08 


±0.08 


±0.06 


±0 05 

_l_ W. WO 


±0 07 

_I_W.W i 


3L 


1.19 


1.23 


1.26 


1.30 


1.19 


1.16 


1.16 


1.14 




+0 1 3 

_I_W. ±kJ 


±0.15 


±0.15 


±0.14 


±0.13 


±0.10 


+0 09 

_l_ w . w u 


+0 1 3 

_1_ w . ± o 


CL 


2.72 


3.17 


2.59 


2.30 


2.14 


2.40 


2.21 


1.63 




±0.43 


±0.71 


±0.53 


±0.40 


±0.48 


±0.23 


±0.22 


±0.19 


EL 


0.22 


0.22 


0.25 


0.30 


0.35 


0.53 


0.64 


0.69 




±0.01 


±0.01 


±0.01 


±0.01 


±0.02 


±0.02 


±0.02 


±0.03 






korrigierte Rate R 


2+r,korr. lH % 






JL 


6.2 


7.4 


8.9 


11.3 


14.2 


16.7 


23.4 


28.0 




±1.1 


±1.3 


±1.4 


±1.6 


±2.0 


±2.3 


±3.0 


±4.7 


JC 


11.0 


12.1 


13.2 


15.3 


18.9 


20.5 


31.0 


32.8 




±1.2 


±1.4 


±1.4 


±1.6 


±1.9 


±2.2 


±2.9 


±4.3 


JB 


4.7 


6.5 


7.5 


9.7 


14.4 


17.1 


27.0 


33.1 




±1.0 


±1.3 


±1.4 


±1.5 


±2.0 


±2.2 


±3.0 


±4.5 


3L 


3.0 


4.0 


4.3 


5.8 


8.1 


8.6 


14.8 


19.8 




±1.0 


±1.2 


±1.2 


±1.4 


±1.7 


±1.8 


±2.6 


±4:A 


CL 


0.9 


1.3 


1.0 


1.4 


1.8 


8.0 


13.9 


13.7 




±0.8 


±1.0 


±0.8 


±0.9 


±1.1 


±2.1 


±3.2 


±3.9 


EL 


14.6 


15.0 


17.1 


21.2 


25.3 


37.0 


45.3 


47.2 




±1.0 


±1.0 


±1.1 


±1.2 


±1.4 


±1.9 


±2.4 


±3.6 



Tabelle 5.10: Zusammenfcissung des Vergleichs verschiedener 
Jetalgorithmen 
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zwischen den beiden Monte-Carlo-Niveaus. Dies bedingt auch die unakzeptablen 
Korrekturfaktoren mit Werten zwischen 1.63 und 3.17. Die sich daraus ergebende 
korrigierte Rate zeigt ein merkwürdiges Verhalten. Während die Rate in den 
unteren fünf Bins nahezu konstant ist, tritt ein Anstieg zwischen Bin fünf und acht 
um einen Faktor sieben bis acht auf. 

Beim E Schema ergeben sich sehr große Raten auf Daten und Detektorniveau. Die 
Ubereinstimmung hier ist gut. Die Raten auf Partonniveau sind um ein Vielfaches 
kleiner und die Abhängigkeit ist hier wesentlich stärker. Somit ergeben sich 
unakzeptable, abhängige Korrekturfaktoren zwischen 0.22 und 0.69. 

Zusammenfassend ergeben sich starke Unterschiede bei Verwendung anderer Algorithmen 
oder anderer Einstellparameter. Der Cone Algorithmus und das E Schema zeigen 
merkwürdige Abhängigkeiten insbesondere beim Korrekturfaktor. Eine endgül- 
tige Bewertung dieser Algorithmen kann aber erst nach einem Vergleich mit den 
Theorieraten erfolgen. Dies ist zur Zeit noch nicht möglich. Für die vier JADE 
Versionen zeigen sich zu große Abweichungen insbesondere bei Verwendung des Algorithmus 
im hadronischen Schwerpunktsystem. Die korrigierte Rate liefert jedoch akzeptable 
Werte. 



5.7 Untersuchungen auf dem Hadronniveau 



Der Korrekturfaktor beschreibt die Ratenänderung, die durch mehrere Effekte eintreten 
kann. Neben der Hadronisierung und Fragmentation können auch Detektoreffekte 
eine Rolle spielen. Um den Anteil dieser beiden Schritte an den Migrationen zu 
untersuchen, spalten wir unsere Betrachtung auf. Zur Untersuchung der durch die 



Hadronisierung verursachten Migrationen sind in Tabelle |5.11| die Migrationen vom 
Parton- zum Hadronniveau verzeichnet. Die Migrationen durch Detektoreffekte oder 



Fehl- Erkennung des gestreuten Elektrons sind in Tabelle |5.12| beim Ubergang von 
Hadron- auf Detektorniveau notiert. Bei diesem Schritt ist auch zu berücksichtigen, 
daß die Akzeptanz auf Parton- und Hadronniveau im gesamten Raum eins ist, während 
wir auf Detektorniveau Akzeptanzlücken haben. Daher benutzen wir nur auf Detektorniveau 
ein Pseudoteilchen als Näherung für den Protonrest. 



Vergleichen wir diese beiden Tabellen mit den Gesamtmigrationen in Tabelle |5.3| , 
so stellen wir fest, daß in beiden Teilschritten die Reinheit in den hohen Bins in 
etwa gleich ist. Gegenüber den Parton-Cluster-Migrationen hat sich die Reinheit nur 
leicht verbessert. In den niedrigen Bins ist die Reinheit beim Ubergang zwischen 
Parton- und Hadronniveau besser, ansonsten gilt das gleiche wie für die hohen Bins. 
Der Korrekturfaktor läßt ebenfalls nur den Schluß zu, daß die Probleme in beiden 
Bereichen gleich stark sind. 



Zum Vergleich der Raten ziehen wir die Partonrate von Tabelle |5.3| und die 
Hadron- und Detektorrate von Tabelle |5.12| heran, da zur Berechnung der Tabelle p.ll 



auch die Ereignisse verwendet wurden, die auf Detektorniveau die Standardschnitte 
nicht erfüllen. Die Hadronrate liegt in fünf Bins zwischen den beiden anderen Raten, 
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Bin 


Parton 


Bin 


Parton 


Bin 


Parton 




1 


1+1 


2+1 




2 


1+1 


2+1 




3 


1+1 


2+1 


Ha- 


1+1 


14825 


327 


Ha- 


1+1 


10839 


251 


Ha- 


1+1 


10095 


291 


dron 


2+1 


433 


508 


dron 


2+1 


325 


385 


dron 


2+1 


351 


442 



Bin 


Parton 


Bin 


Parton 


Bin 


Parton 




4 


1+1 


2+1 




5 


1 + 1 


2+1 


6 


1+1 


2+1 


Ha- 


1 + 1 


9151 


306 


Ha- 


1+1 


6262 


251 


Ha- 


1+1 


7067 


123 


dron 


2+1 


301 


528 


dron 


2+1 


230 


440 


dron 


2+1 


363 


849 



Bin 


Parton 


Bin 


Parton 




7 


1 + 1 


2+1 




8 


1 + 1 


2+1 


Ha- 


1+1 


3990 


349 


Ha- 


1+1 


1675 


128 


dron 


2+1 


301 


805 


dron 


2+1 


128 


429 



Bin 


1 


2 


3 


4 


5 


6 


7 


8 


Rate Parton 
in% 


5.2 
±0.2 


5.4 
±0.3 


6.6 
±0.3 


8.1 
±0.3 


9.6 
±0.4 


14.6 
±0.4 


21.2 
±0.6 


23.6 
±0.9 


Rate Hadron 

in 7i 


5.8 

±0.2 


6.0 

±0.3 


7.1 

±0.3 


8.1 

±0.3 


9.3 

±0.1 


13.9 

±0.1 


20.3 

±0.6 


23.6 

±0.9 


Korrekt ur- 
faktor 


0.89 
±0.06 


0.90 
±0.07 


0.92 
±0.07 


1.01 
±0.07 


1.03 
±0.08 


1.05 
±0.06 


1.04 
±0.06 


1.00 
±0.08 


Reinheit in % 


54.0 


54.2 


55.7 


63.7 


65.7 


70.0 


72.8 


77.0 



Tabelle 5.11: Migrationen und Korrekturfaktoren bei Lepto mit 
Standardschnitten, Vergleich von Parton- und Hadronniveau in Bins 
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Bin 


Hadron 


Bin 


Hadron 


Bin 


Hadron 


1 


1+1 


2+1 


2 


1+1 


2+1 


3 


1+1 


2+1 


Clus- 


1+1 


12395 


389 


Clus- 


1+1 


9358 


308 


Clus- 


1+1 


8934 


327 


ter 


2+1 


419 


358 


ter 


2+1 


324 


279 


ter 


2+1 


298 


325 



Bin 


Hadron 


Bin 


Hadron 


Bin 


Hadron 


4 


1+1 


2+1 


5 


1+1 


2+1 


6 


1+1 


2+1 


Clus- 


1+1 


8128 


325 


Clus- 


1+1 


5183 


237 


Clus- 


1+1 


5968 


114 


ter 


2+1 


282 


368 


ter 


2+1 


242 


316 


ter 


2+1 


270 


590 



Bin 


Hadron 


Bin 


Hadron 


7 


1+1 


2+1 






1 + 1 


2+1 


Clus- 


1+1 


3368 


297 


Clus- 


1+1 


1304 


137 


ter 


2+1 


223 


616 


ter 


2+1 


90 


295 



Bin 


1 


2 


3 


4 


5 


6 


7 


8 


Rate Hadron 


5.5 


5.7 


6.6 


7.6 


8.8 


13.9 


20.3 


23.7 


in% 


±0.2 


±0.3 


±0.3 


±0.3 


±0.4 


±0.5 


±0.6 


±1.0 


Rate Detektor 


5.7 


5.9 


6.3 


7.1 


8.9 


11.9 


18.6 


21.1 


in 'X 


±0.2 


±0.3 


±0.3 


±0.3 


±0.1 


±0.1 


±0.6 


±1.0 


Korrekt ur- 


0.96 


0.97 


1.05 


1.07 


0.99 


1.17 


1.09 


1.12 


faktor 


±0.07 


±0.08 


±0.09 


±0.08 


±0.09 


±0.08 


±0.07 


±0.10 


Reinheit in % 


46.1 


46.3 


52.2 


56.6 


56.6 


68.6 


73.4 


76.6 



Tabelle 5.12: Migrationen und Korrekturfaktoren bei Lepto mit 
Standardschnitten, Vergleich von Hadron- und Detektorniveau in Bins 



70 KAPITEL 5. ERGEBNISSE DER UNTERSUCHUNGEN 



Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


1 


1+1 


2+1 




2 


1+1 


2+1 




3 


1+1 


2+1 


Ha- 


1 + 1 


1890 


37 


Ha- 


1+1 


8191 


443 


Ha- 


1+1 


2164 


382 


dron 


2+1 


39 


19 


dron 


2+1 


409 


942 


dron 


2+1 


318 


1037 



Hadronanzahl im Ereignis 


kleiner 25 


25 bis 50 


größer > 50 


Rate Parton 


2.8 


13.9 


36.4 


in % 


±0.4 


±0.4 


±0.8 


Rate Hadron 


2.9 


13.5 


34.7 


in % 


±0.4 


±0.4 


±0.8 


Korrektur- 


0.97 


1.03 


1.05 


faktor 


±0.26 


±0.06 


±0.05 


Reinheit in % 


32.8 


69.7 


76.5 



Tabelle 5.13: Migrationen und Korrekturfaktoren bei Lepto mit 
Standardschnitten, in Bins bezüglich der Hadronmultiplizität 

in zwei Bins entsprechen sie den Partonraten und in nur einem Bin ist die Rate 
höher als die auf Parton- und Detektorniveau. Die Ursache der Migrationen scheint 
somit ein Effekt zu sein, der sowohl beim Ubergang vom Parton- zum Hadronniveau 
alsauch bei Berücksichtigung der Detektoreffekte auftritt. Um dies zu beweisen, sind 
jedoch weitere Studien und ein besseres Verständnis der Hadronisierung und der 
Detektoreffekte nötig. 

Um den Einfluß der Hadronanzahl zu untersuchen, habe ich die Ereignisse in Bins 
entsprechend ihrer Hadronanzahl geteilt. Die sich daraus ergebenden Migrationen 
sind in Tabelle |5.13| verzeichnet. Wir erkennen, daß die Reinheit mit der Hadronanzahl 



stark ansteigt. Unsere Probleme scheinen also nicht durch die von Niveau zu Niveau 
steigende Partikelanzahl ausgelöst zu werden. Der Korrekturfaktor ist im Rahmen 
der Fehler in allen Bereichen gleich. Der Grund für das Ansteigen der Raten mit der 
Multiplizität ist unbekannt und deutet die Wichtigkeit von weiteren Untersuchungen 
der Hadronisierungseffekte an. Diese würden jedoch den Rahmen der Diplomarbeit 
sprengen. 
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5.8 Einführung zusätzlicher Schnitte 
5.8.1 Jetwinkel öjct 

Aufgrund der Hinweise von Kapitel |5.4| schauen wir uns nun die Winkelverteilung der 
harten Jets an. In Abbildung ^.12| a sind die Verteilungen für das Parton- und das 
Detektorniveau für alle Ereignisse und in Bild b die inklusiven Verteilungen für die 
2+1 Ereignisse alleine gezeigt. Beide Verteilungen sind auf gleiche Flächen normiert. 

Wir erkennen schon hier, daß im Vorwärts- und im Rückwärtsbereich die Beschreibungen 
nicht gut übereinstimmen. Deshalb sind die beiden Bereiche für alle Ereignisse in 
c und d nochmal vergrößert dargestellt. Dazu wurde die gleiche Normierung wie in 
Bild a benutzt. 

Die Beschreibung der Bereiche kleiner 70 Grad und größer 150 Grad zeigt Abweichungen 
zwischen den Monte-Carlo-Niveaus. Im Bereich dazwischen stimmt die Verteilung der 
2+1 Ereignisse gut und bei allen Ereignissen unter Berücksichtigung der dann leicht 
geänderten Normierung einigermaßen überein. Es scheint jedoch eine Verschiebung 
von großen Winkeln auf Partonniveau zu kleineren auf Detektorniveau vorzuliegen. 

Betrachten wir jedoch zuerst den Randbereich bei hohen Winkeln (Abbildung 
5.12| d). Das LAr Kalorimeter endet bei 153 Grad und der auf Hadronen empfindlichste 



Bereich bereits etwas eher. Dies hat jedoch auch schon Auswirkungen bei leicht 
kleineren Winkeln, da dann der Jet einseitig abgeschnitten wird und sich der Jetwinkel 
leicht zu kleineren Werten verschiebt. Die Abweichungen hier sind also verständlich 
und es empfiehlt sich ein Schnitt auf Jetwinkel kleiner als 150 Grad. 

Im Bereich kleiner 20 Grad liegt die Situation etwas anders. Hier übersteigt die 
Detektorverteilung die auf Partonniveau. Der Detektoranfang liegt bei vier Grad 
und dort ist auch ein entsprechender Abfall auf beiden Niveaus zu erkennen. Ein 
anderer Effekt erhöht jedoch die Anzahl Jets unter kleinen Winkeln. Der Protonrest 
verläßt den Detektor ungesehen durch das Strahlrohr. Dieser Rest kann durch 
Abstrahlung oder Hadronisierung in mehrere Hadronen zerfallen, die ähnlich den 
harten Jets einen gewissen Bereich im Detektor überdecken. Diese Hadronen fallen 
dann unter kleinen Winkeln in das Kalorimeter und hinterlassen dort ihre Energie. 
Wird ein harter Jet unter kleinen Winkeln gestreut, können sich diese überlagern und 
es ergeben sich Verfälschungen in der Energie und dem Winkel der Jets. Dies ist in 



Abbildung |5.13| zu erkennen. Dort ist die Energie des Protonrestes auf Partonniveau 
gegen die Energie des Jets aufgetragen, der aus dem Pseudoteilchen und eventuell 
zusätzlichen, rekonstruierten Clustern gebildet wurde. Die Größe der Boxen ist 
wieder logarithmisch. In den Teilbildern wurden jeweils unterschiedliche Jetwinkel 
Bereiche für den unteren harten Jet ausgewählt. 

Im Bild links oben sind die Ereignisse gezeigt, deren unterster harter Jet im 
Bereich zwischen sieben Grad und zehn Grad liegt. Wir erkennen hier eine sehr breite 
Verteilung, die auf eine mangelnde Rekonstruktion schließen läßt. In der Abbildung 
daneben ist der Bereich zwischen zehn und zwölf Grad gezeigt. Hier zeigt sich eine 
nur leicht verbesserte Korrelation. Erst in den beiden unteren Graphen ist eine 
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0.035 




Abbildung 5.12: Vergleich der O^et Verteilungen auf Parton- und 
Detektorniveau a) alle Ereignisse, b) 2+1 Ereignisse, c) und d) sind 
Vergrößerungen von a. Alle Plots sind mit Standardschnitten berechnet worden 
und auf gleiche Flächen normiert. Dargestellt sind Parton- (durchgezogen) und 
Detektorniveau (gestrichelt). 
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Abbildung 5.13: Vergleich der Protonrestenergien für Ereignisse mit Jets im 
Bereich kleiner öjct- a) 7° < öjet < 10°, b) 10° < öjet < 12°, c) 12° < öjet < 15°, d) 
15° < öjet < 20°. Die Größe der Boxen ist logarithmisch. 
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Abbildung 5.14: Vergleich der Jetwinkel zwischen Parton- und 
Detektorniveau im Bereich 20° < öjet < 145°. a) alle Ereignisse, b) 2+1 
Ereignisse. 



zufriedenstellende Korrelation erreicht. Die Bereiche sind hier zwölf bis 15 und 15 
bis 20 Grad. Dies legt einen Winkelschnitt von 15 Grad nahe. 

Bei Anwendung des Jetwinkelschnittes werden solche Ereignisse zur totalen Anzahl 
gezählt, die einen Winkel in dem angegebenen Bereich haben, zur 2+1 Jetanzahl 
jedoch nur dann, wenn beide harten Jets die Bedingung erfüllen. Dadurch kann 
jedoch die Situation eintreten, daß in einem Ereignis zwar zwei harte Jets gefunden 
werden, es aber als 1 + 1 Ereignis gehandhabt wird, da nur ein Jet den Winkelschnitt 
erfüllt. Dieses Vorgehen ist aber dem kompletten Verwerfen solcher Ereignisse vorzuziehen, 
da auch Ereignisse, deren zweiter Jet z.B. durch Detektorakzeptanzen erst garnicht 
detektiert wurde, beim totalen Wirkungsquerschnitt berücksichtigt werden. 

Die Abhängigkeit der Rate von der Größe des unteren Schnittes ist in den Tabellen 
5.14| und ß.l5| für die oberen Bins gezeigt. Wir erkennen ein Abnehmen der Raten 
und ein Ansteigen der Korrekturfaktoren. 

Betrachten wir nun den mittleren Jetwinkelbereich. In Abbildung ß.l4| ist die 
Korrelation des Jetwinkels auf Parton- gegen den auf Detektorniveau gezeigt. In 
Bild a sind alle Ereignisse verzeichnet. Die Ubereinstimmung ist gut zu sehen, es 
erfolgen jedoch auch viele Fehl-Erkennungen. Bei den 2+1 Ereignissen, die in Bild b 
gezeigt sind, ist die Korrelation ebenfalls gut, die Statistik ist jedoch begrenzt. 

In Anbetracht der schlechten Beschreibung der Jetwinkel im Bereich kleiner und 
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Bin 6 




o 







2 


4 


6 


8 


10 


Parton- 


14.6 


14.6 


14.1 


13.0 


11.6 


10.0 


rate 


[%] 


±0.4 


±0.4 


±0.4 


±0.4 


±0.3 


±0.3 


Detektor- 


12.4 


12.4 


11.9 


10.4 


8.5 


6.7 


rat 


[</ 


] 


±0. 1 


±0.1 


±0.1 


±0.1 


±0.1 


±0.3 


Korrektur- 


1.19 


1.19 


1.19 


1.26 


1.38 


1.51 


faktor 


±0.07 


±0.07 


±0.07 


±0.08 


±0.09 


±0.11 


Ojet 


o 




12 


14 


16 


18 


20 




Parton- 


8.9 


7.6 


6.5 


5.6 


4.9 




rate 


[%] 


±0.3 


±0.3 


±0.3 


±0.2 


±0.2 




Detektor- 


5.5 


4.5 


3.8 


3.0 


2.7 




rate 


[%] 


±0.3 


±0.3 


±0.3 


±0.2 


±0.2 




Korrektur- 


1.63 


1.72 


1.75 


1.89 


1.86 




faktor 


±0.13 


±0.16 


±0.17 


±0.21 


±0.22 




Bin 7 












2 


4 


6 


8 


10 


Parton- 


21.2 


21.2 


20.1 


17.4 


14.5 


12.3 


ratc 


[%] 


±0.6 


±0.6 


±0.5 


±0.5 


±0.5 


±0.4 


Detektor- 


19.1 


19.1 


17.8 


14.2 


10.9 


8.5 


rate 


[%] 


±0.6 


±0.6 


±0.6 


±0.6 


±0.5 


±0.5 


Korrektur- 


1.11 


1.11 


1.13 


1.23 


1.34 


1.45 


faktor 


±0.06 


±0.06 


±0.07 


±0.08 


±0.10 


±0.12 


Ojct 






12 


14 


16 


18 


20 




Parton- 


10.2 


8.6 


7.2 


6.0 


5.2 




rate 


[%] 


±0.4 


±0.4 


±0.4 


±0.3 


±0.3 




Detektor- 


6.6 


5.3 


4.3 


3.8 


3.1 




rate 


[%] 


±0.4 


±0.4 


±0.4 


±0.3 


±0.3 




Korrekt ur- 


1.56 


1.65 


1.68 


1.60 


1.70 




faktor 


±0.15 


±0.18 


±0.20 


±0.21 


±0.25 





Tabelle 5.14: Raten- und Korrekturfaktorbetrachtung bei 9 jet Schnitt, in 
Bins 
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Bin 8 


6'jct 


o 







2 


4 


6 


8 


10 


Parton- 


23.6 


23.6 


22.4 


20.3 


17.0 


13.9 


rate 


[%] 


±0.9 


±0.9 


±0.9 


±0.8 


±0.8 


±0.7 


Detektor- 


21.8 


21.8 


21.3 


18.3 


14.0 


10.9 


rat c 


%] 


±1.0 


±1.0 


±1.0 


±0.9 


±0.8 


±0.8 


Korrektur- 


1.09 


1.09 


1.06 


1.12 


1.22 


1.28 


faktor 


±0.09 


±0.09 


±0.09 


±0.10 


±0.12 


±0.15 


^jet 


o 




12 


14 


16 


18 


20 




Parton- 


11.5 


9.4 


8.0 


7.1 


6.6 




rate 


[%] 


±0.7 


±0.6 


±0.6 


±0.6 


±0.6 




Detektor- 


8.2 


6.7 


5.4 


4.5 


4.2 




rate 


[%] 


±0.7 


±0.6 


±0.6 


±0.6 


±0.6 




Korrektur- 


1.40 


1.41 


1.50 


1.57 


1.57 




faktor 


±0.19 


±0.22 


±0.26 


±0.31 


±0.33 





Tabelle 5.15: Raten- und Korrekturfaktorbetrachtung bei 9jet Schnitt, in 
Bins 



großer Winkel, ist ein Schnitt in dieser Größe sinnvoll. Wir erreichen dadurch 
eine gute Beschreibung im mittleren Winkclbcrcich. Dies ist aufgrund der guten 
Korrelation dort erlaubt. Der obere Schnitt sollte einen Wert von 145 bis 150 Grad 
haben. Bei der unteren Grenze ist die Angabe eines Wertes komplizierter. Ein 
zu tiefer Wert behält einen zu großen Teil der Problemzone bei. Um eine gute 
Ubereinstimmung der beiden Monte- Carlo-Niveaus zu erreichen, wäre ein Schnitt 
von mindestens 20 Grad angebracht. Dabei verliert man jedoch viel Statistik, so daß 
sich der entsprechende Fehler vergrößert, und der Korrekturfaktor steigt, was eine 
Ratenmessung problematisch werden läßt. Der Schnitt sollte deshalb im Bereich 
zwischen zehn und 20 Grad je nach Anforderungen an die Messung liegen. Ist eine 
hohe Statistik erwünscht, so sollten schwächere Grenzen gewählt werden, wobei sich 
der systematische Fehler dann erhöht. Wir wählen die Grenzen 15 und 145 Grad. 



5.8.2 Skalenvariable Zp 

Eine weitere Größe, die ich untersucht habe, ist die Verteilung in der Skalenvariable 
Zp. Sie ist definiert für 2±1 Jet Ereignisse durch den Viererimpuls des einlaufenden 
Protons P und die Viererimpulse der beiden harten Jets Pj 

'P,^ - p (5-16) 

p • (Pl ± P2) 

= i(l-cos(ö*)) (5.17) 



5.8. EINFÜHRUNG ZUSÄTZLICHER SCHNITTE 



77 



0.03 
0.025 

0.02 
0.015 

0.01 
0.005 



Partor 

Detektor 

Daten 







0.05 



0.15 



0.2 



0.25 

a) 



0.3 



0.35 



0.4 



0.45 



0.5 

z„ 



Abbildung 5.15: Vergleich der Zp Verteilungen auf Parton- und 
Detektorniveau und bei den Daten nach Standardschnitten. Dargestellt 
sind Parton- (durchgezogen) und Detektorniveau (gestrichelt) der Lepto-Monte- 
Carlos und die Daten (gepunktet). Die Verteilungen sind auf gleiche Flächen 
normiert. 



Der Winkel 6* ist der Winkel zwischen den beiden harten Jets im hadronischen 
Schwerpunktsystem. Es gelten die beiden Beziehungen 



1 ~ ^p,2 

180 - ei 



(5.18) 
(5.19) 



Abbildung ß.l5| zeigt die Verteilung des jeweils minimalen 



mm z. 



für Parton- und Detektorniveau und für die Daten. Die Verteilungen sind auf 
gleiche Flächen normiert. Wir erkennen eine gute Ubereinstimmung. Ein Schnitt 
in dieser Größe wird keine Änderung des Korrekturfaktors hervorrufen, da Parton- 
und Detektorniveau in gleicher Weise getroffen werden. 

Um jedoch einen Einfluß des Schnittes auf die Migrationen und die Reinheit zu 
untersuchen, berechnen wir die Raten in Abhängigkeit eines minimalen Zp Schnittes. 
Tabelle p.l6| zeigt das Ergebnis für die Bins drei, sechs und acht für drei unterschiedliche 
Schnittwerte. Diese Werte können mit denen in Tabelle |5.3| verglichen werden. Für 
die einzelnen Schnitte gilt, daß die Gesamtereigniszahl gleich bleibt, wobei jedoch 
einige der 2+1 Ereignisse durch die Schnitte nicht mehr als 2+1 Ereignisse gezählt 
werden. Auch hier ist dieses Verfahren sinnvoller als ein Verwerfen der Ereignisse, 
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da Ereignisse, deren zweiter Jet nicht gefunden wurde, ebenfalls zur Gesamtanzahl 
gerechnet werden. Dadurch werden die Raten kleiner. Wir erkennen, daß der Einfluß 
des Schnittes auf den Korrekturfaktor im Rahmen der Fehler wie erwartet keine 
Änderung bewirkt, wobei das niedrige Bin etwas empfindlicher reagiert. Die 
Reinheit zeigt auch keine erkennbare Veränderung. 

Ein Schnitt in Zp erscheint aus diesen Gründen nicht nötig, aber auch nicht 
verboten. 



5.8.3 Transversaler Impuls pl 

Der transversale Impuls der Jets im hadronischen Schwerpunktsystem ist ebenfalls 
untersucht worden. Wir erwarten hier nur geringe Impulswerte für 1+1 Jet Ereignisse, 
da der Protonrest den Detektor weiterhin in 2;- Richtung verläßt. Bei 2+1 Jet Ereignissen 
sollten die transversalen Impulse der beiden harten Jets die gleiche Größe haben, da 
die transversale Impulssumme der harten Jets ebenfalls klein sein sollte. 



Die zugehörigen Verteilungen sind in Abbildung ^.16| dargestellt. Im Bild links 
oben ist die Verteilung für alle Ereignisse in den drei Niveaus, daneben für die 2+1 
Ereignisse gezeigt. Wir erkennen, daß der mittlere pl Wert für die 2+1 Verteilungen 
klar von Null abweicht, während die Verteilung für alle Ereignisse bei Null einen 
Peak zeigt. Im unteren Teil ist die Korrelation der beiden p^ Werte der harten Jets 
für 2+1 Ereignisse abgebildet. Die Größe der Boxen ist wieder logarithmisch. 

Die Korrelationen auf Partonniveau sind gut. Auf Detektorniveau und bei den 
Daten zeigt die Gegenüberstellung schlechtere Wert. Dies liegt an der Ungenauigkeit 
des transversalen Impulses im Laborsystem, die sich auch hier zeigen, da die Lo- 
rentztransformation fast ausschließlich in z Richtung stattfindet. Im Rahmen der in 
Kapitel ß.2.4| betrachteten Impulsverteilung ist die Korrelation auch hier gut. 



Vergleichen wir die Niveaus in Abbildung |5.16| a, so fällt auf, daß auf Detektorniveau 



und bei den Daten weniger Jets mit einem transversalen Impuls kleiner als 1 GeV 
gefunden werden. Dies ist aufgrund der begrenzten Energieaufiösung erklärbar. 
Insbesondere ist die Berechnung der Lorentztransformation auf die Genauigkeit der 
Rekonstruktion des hadronischen Endzustandes angewiesen. Dies verursacht wie 
erwähnt auch die Verbreiterung der Korrelationen in den Bildern c bis e, die in der 
gleichen Größenordnung liegen. Bei größeren transversalen Impulsen stimmen jedoch 
alle drei Verteilungen unter Beachtung der geänderten Normierung einigermaßen 
überein. Dies gilt ebenso für die 2+1 Ereignisse, da deren transversale Impulse, wie 
bereits erwähnt, größer sind. Die Beschreibung der Daten gelingt auch hier durch 
das Detektorniveau sehr gut. 

Bei einem Schnitt in den transversalen Impuls in der Größe von 2 GeV erwarten 
wir also den Wegfall vieler 1+1 Jet Ereignisse, jedoch nur weniger mit zwei harten 



Jets. Dies wird in Tabelle ß.l7| an den Migrationsmatrizen deutlich. Die Raten 
sind stark angestiegen. Beim Vergleich mit Tabelle |5.^ stellen wir einen Anstieg des 
Korrekturfaktorwertes fest. Die Variation ist mit 1.26 bis 1.40 kleiner und alle Werte 
sind im Rahmen ihrer Fehler mit dem Mittelwert 1.33 verträglich. Die Reinheit der 



5.8. EINFÜHRUNG ZUSÄTZLICHER SCHNITTE 



79 



Bin 3 


Parton 


Bin 3 


Parton 


Bin 3 


Parton 


Schnitt 1 


1 + 1 


2+1 


Schnitt 2 


1 + 1 


2+1 


Schnitt 3 


1 + 1 


2+1 


Clus- 


1+1 


9330 


188 


Clus- 


1+1 


9563 


123 


Clus- 


1+1 


9668 


85 


ter 


2-1-1 


168 


198 


ter 


2+1 


91 


107 


ter 


2+1 


60 


71 


Bin 6 


Parton 


Bin 6 


Parton 


Bin 6 


Parton 


Schnitt 1 


1+1 


2+1 


Schnitt 2 


1+1 


2+1 


Schnitt 3 


1+1 


2+1 


Clus- 


1+1 


6367 


306 


Clus- 


1+1 


6678 


195 


Clus- 


1+1 


6852 


140 


ter 


2+1 


177 


392 


ter 


2+1 


118 


251 


ter 


2+1 


67 


183 


Bin 8 


Parton 


Bin 8 


Parton 


Bin 8 


Parton 


Schnitt 1 


1+1 


2+1 


Schnitt 2 


1+1 


2+1 


Schnitt 3 


1+1 


2+1 


Clus- 


1+1 


1431 


95 


Clus- 


1+1 


1536 


82 


Clus- 


1+1 


1596 


69 


ter 


2+1 


71 


229 


ter 


2+1 


52 


156 


ter 


2+1 


44 


117 



Bin 


3 


6 


8 


Schnitt 


1 


2 


3 


1 


2 


3 


1 


2 


3 


Parton 


3.9 


2.3 


1.6 


9.6 


6.2 


4.5 


17.7 


13.0 


10.2 


in % 


±0.2 


±0.2 


±0.2 


±0.1 


±0.3 


±0.3 


±0.9 


±0.8 


±0.8 


Detektor 


3.7 


2.0 


1.3 


7.9 


5.1 


3.5 


16.4 


11.4 


8.8 


in % 


±0.2 


±0.2 


±0.2 


±0.4 


±0.3 


±0.3 


±0.9 


±0.8 


±0.7 


Korrek- 


1.05 


1.16 


1.19 


1.23 


1.21 


1.29 


1.08 


1.14 


1.16 


turf a. 


±0.11 


±0.16 


±0.20 


±0.10 


±0.12 


±0.16 


±0.12 


±0.15 


±0.17 


Reinheit 


54.1 


54.0 


54.2 


68.9 


68.0 


73.2 


76.3 


75.0 


72.7 


in% 





















Tabelle 5.16: Migrationen und Korrekturfaktoren mit verschiedenen Zp 
Schnitten Exemplarisch werden drei Bins mit drei verschiedenen Zp Schnitten 
gezeigt. Schnitt 1 ist hier Zp > 0.05, Schnitt 2 Zp> 0.10 und Schnitt 3 Zp > 0.15. 
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0.7 ^ 

0.6 - 

0.5 - 

0.4 - 
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a) 
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D 



: D 
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0.1 
0.08 
0.06 
0.04 
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Abbildung 5.16: Vergleich der Verteilungen auf Parton- und 

Detektorniveau und bei den Daten nach Standairdschnitten. a) Verteilung 
für alle Ereignisse, dargestellt sind Parton- (durchgezogen) und Detektorniveau 

(gestrichelt) und die Daten (gepunktet), b) Verteilung für 2+1 Jet Ereignisse, 
c,d,e) Korrelation des transversalen Impulses für beide harte Jets mit logarithmischer 
Boxengröße hei 2+1 Ereignissen für Parton- (c), Detektorniveau (d) und Daten (e) 
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2+1 Jetklasse ist in allen Bins deutlich gestiegen. Dieser Schnitt verbessert auch 
aufgrund der geänderten Normierung die Verteilungen der transversalen Impulse. 

Die Qualität des Schnittes trügt jedoch. Dies sieht man an der Ratenänderung in 
Abhängigkeit der Größe des Schnittes. Wir erkennen hier sehr große Korrekturfaktoren. 
Der Unterschied zwischen den Detektorraten bei den Migrationsmatrizen und bei der 
Berechnung in der Tabelle ^.19| liegt darin begründet, daß bei den Migrationsmatrizen 
nur die Ereignisse gezeigt werden, die die Schnitte auf beiden Niveaus erfüllen. Viele 
Ereignisse haben aber, wie wir in Abbildung ß.l6| a gesehen haben, einen verschwindenden 
transversalen Impuls auf Partonniveau, jedoch einen höheren auf Detektorniveau. 
Eine Stabilisierung auf hohem Niveau tritt erst ein, wenn die Verteilung der transversalen 
Impulse auf beiden Niveaus übereinstimmt. 

Ein Schnitt in dieser Größe ist gefährlich, da die Verteilung in keinem Bereich gut 
beschrieben wird. 



5.8.4 Hadronische Schwerpunktenergie s 

Die bereits in Kapitel |3.1.1j in Gleichung |3.5| eingeführte Variable s kann in der Praxis 
für 2+1 Jet Ereignisse durch zwei unterschiedliche Methoden berechnet werden. Die 
eine berechnet die invariante Masse des harten hadronischen Endzustandes und die 
andere benutzt die Jetwinkel im hadronischen Schwerpunktsystem. Während die 
zweite Berechnung eindeutig vom Jetalgorithmus abhängig ist, trifft dies für die erste 
Methode nur in abgeschwächter Form zu. Das Problem bei der Berechnung über 
den hadronischen Endzustand ist, daß nur die Cluster verwendet werden dürfen, die 
zu den harten Jets gehören, nicht jedoch die es Protonrestes. Dazu ist jedoch ein 
Jetalgorithmus nötig. 

Die Berechnung kann dann durch die Formeln 



^had 



E Pj (5-20) 

*\hadr. Partikel / 

gjet ^ w^^{v*i+V*2) (5.21) 



s- 



mit den Pseudorapiditäten rj* der beiden harten Jets im hadronischen Schwerpunktsystem 
erfolgen. 

In Abbildung p.l7| sind die Werte gegenübergestellt, die mit den unterschiedlichen 
Methoden berechnet wurden. Aufgetragen wurde die Wurzel von s. Dabei ist oben 
für Daten, Parton- und Detektorniveau die Korrelation dargestellt und in Bild d der 
Absolutbetrag der Differenz aufgetragen. Wir erkennen, daß die Verteilung in allen 
Niveaus ähnlich ist, nur im Bereich kleiner Unerschiede macht sich die zusätzliche 
Begrenzung durch die Energieauflösung bemerkbar. Die Korrelationen zwischen 
den beiden Berechnungen ist jedoch nicht sehr gut. Die Skalierung der Boxen ist 
logarithmisch. 

Die zugehörigen Migrationen wurden berechnet und in Tabelle |5.20| exemplarisch 
für die Bins drei, sechs und acht mit drei verschiedenen Schnitten aufgetragen. 
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Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


1 


1+1 


2+1 


2 


1+1 


2+1 


3 


1+1 


2+1 


Clus- 


1+1 


607 


221 


Clus- 


1+1 


419 


172 


Clus- 


1+1 


453 


204 


ter 


2+1 


108 


326 


ter 


2+1 


72 


247 


ter 


2+1 


78 


269 



Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


4 


1+1 


2+1 


5 


1 + 1 


2+1 


6 


1+1 


2+1 


Clus- 


1+1 


101 


230 


Clus- 


1 + 1 


320 


176 


Clus- 


1+1 


811 


369 


ter 


2+1 


87 


330 


ter 


2+1 


61 


274 


ter 


2+1 


105 


556 



Bin 


Parton 


Bin 


Parton 


7 


1+1 


2+1 






1+1 


2+1 


Clus- 


1+1 


543 


283 


Clus- 


1+1 


239 


127 


ter 


2+1 


97 


560 


ter 


2+1 


33 


274 



Bin 


1 


2 


3 


4 


5 


6 


7 


8 


Rate Parton 


43.3 


46.0 


47.1 


53.3 


54.2 


50.2 


56.8 


59.6 


in% 


±1.4 


±1.7 


±1.6 


±1.6 


±1.8 


±1.2 


±1.3 


±1.9 


Rate Detektor 


34.4 


35.1 


34.6 


39.7 


40.3 


35.9 


44.3 


45.6 


in 'X 


±1.1 


±1.6 


±1.6 


±1.6 


±1.8 


±1.2 


±1.3 


±2.0 


Korrektur- 


1.26 


1.31 


1.36 


1.34 


1.34 


1.40 


1.28 


1.31 


faktor 


±0.09 


±0.11 


±0.11 


±0.09 


±0.10 


±0.08 


±0.07 


±0.10 


Reinheit in % 


75.1 


77.4 


77.5 


79.1 


81.8 


84.1 


85.2 


89.3 



Tabelle 5.17: Migrationen und Korrekturfaktoren bei Lepto mit 
Standardschnitten und zusätzlichem > 2 GeV Schnitt, in Q"^ Bins 
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Bin 6 


Ptrain [GeV] 





1 


2 


3 


4 


5 


Parton- 


14.6 


32.2 


46.9 


60.6 


68.9 


73.3 


rate [%] 


±0.4 


±0.7 


±1.0 


±1.1 


±1.2 


±1.2 


Detektor- 


12.4 


14.2 


18.5 


23.7 


27.6 


31.4 


rate ['X] 


±0.1 


±0.5 


±0.6 


±0.8 


±1.0 


±1.2 


Korrektur- 


1.19 


2.28 


2.55 


2.56 


2.50 


2.34 


faktor 


±0.07 


±0.12 


±0.13 


±0.13 


±0.13 


±0.12 


Plrain [GeV] 


6 


7 


8 


9 


10 




Parton- 


75.9 


77.3 


77.9 


75.5 


74.3 




rate [%] 


±1.3 


±1.4 


±1.5 


±1.7 


±1.9 




Detektor- 


34.5 


38.3 


41.0 


42.6 


45.1 




rate [%] 


±1.4 


±1.7 


±1.9 


±2.2 


±2.6 




Korrektur- 


2.21 


2.02 


1.90 


1.77 


1.65 




faktor 


±0.12 


±0.12 


±0.12 


±0.13 


±0.13 




Bin 7 


Kmin [GeV] 





1 


2 


3 


4 


5 


Parton- 


21.2 


37.6 


53.6 


65.5 


72.3 


77.4 


rate [%] 


±0.6 


±0.9 


±1.1 


±1.2 


±1.3 


±1.3 


Detektor- 


19.1 


20.3 


23.5 


26.8 


29.4 


30.7 


rate [%] 


±0.6 


±0.7 


±0.8 


±0.9 


±1.1 


±1.2 


Korrektur- 


1.11 


1.86 


2.29 


2.45 


2.46 


2.52 


faktor 


±0.06 


±0.10 


±0.12 


±0.12 


±0.13 


±0.14 


Kmin [GeV] 


6 


7 


8 


9 


10 




Parton- 


78.1 


78.6 


77.7 


77.9 


76.9 




rate [%] 


±1.4 


±1.5 


±1.7 


±1.8 


±2.0 




Detektor- 


32.0 


31.6 


31.7 


32.0 


33.1 




rate [%] 


±1.4 


±1.6 


±1.8 


±2.0 


±2.3 




Korrektur- 


2.45 


2.49 


2.46 


2.44 


2.33 




faktor 


±0.15 


±0.17 


±0.19 


±0.21 


±0.22 





Tabelle 5.18: Raten- und Korrekturfciktorbetrachtung bei Schnitt, in 
Bins 
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Abbildung 5.17: Vergleich der \/§ 

Verteilungen auf PEU"ton- und Detektorniveau und bei den Daten nach 
Standctrdschnitten a,b,c) Korrelation der Jet- und der Summations-Methode mit 

logarithmischer Boxenskalierung für Parton- (a), Detektorniveau (h) und Daten (c). 
d) Absolutbetrag der Differenz der Methoden. Dargestellt sind Parton- (durchgezogen) 
und Detektorniveau (gestrichelt) und die Daten (gepunktet). 
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Bin 8 







1 


2 


3 


4 


5 


Parton- 


23.6 


39.4 


54.7 


66.1 


73.9 


77.8 


rate [%] 


±0.9 


±1.3 


±1.6 


±1.7 


±1.8 


±1.8 


Detektor- 


21.8 


22.7 


25.2 


27.4 


29.2 


31.2 


rate [%] 


±1.0 


±1.0 


±1.2 


±1.3 


±1.5 


±1.7 


Korrektur- 


1.09 


1.74 


2.17 


2.42 


2.54 


2.50 


faktor 


±0.09 


±0.13 


±0.16 


±0.18 


±0.18 


±0.19 


Kmin [GeV] 


6 


7 


8 


9 


10 




Parton- 


81.1 


81.9 


79.0 


79.4 


79.5 




rate [%] 


±1.8 


±1.9 


±2.1 


±2.3 


±2.5 




Detektor- 


32.3 


34.0 


32.4 


31.4 


32.0 




rate [%] 


±1.9 


±2.1 


±2.2 


±2.4 


±2.6 




Korrektur- 


2.51 


2.41 


2.44 


2.53 


2.49 




faktor 


±0.20 


±0.20 


±0.23 


±0.26 


±0.28 





Tabelle 5.19: Raten- und Korrekturfaktorbetrachtung bei pl Schnitt, in 
Bins 



Der höchste Schnitt 3 ist hier der loseste und gibt deshalb die höchste Statistik. Der 
Korrekturfaktor ändert sich nur im Bin acht, dort vergrößert er sich jedoch stark. 
Die Reinheit verschlechtert sich in allen Bins mit stärkerem Schnitt. 
Ein Schnitt in dieser Variable ist nicht empfehlenswert. 



5.8.5 Rapiditätsdifferenz Ar] 

Ein weiterer Schnitt, den es zu untersuchen lohnt, bezieht sich auf die räumliche 
Verteilung der Jets. Schon anschaulich ist klar, daß Verwechslungen verstärkt auftreten, 
wenn die Jets nahe beieinander liegen und sich so nach der Hadronisierung teilweise 
durchdringen. Um diese Ereignisse zu unterdrücken ist ein Schnitt sinnvoll, der 
einen minimalen Rapiditätsunterschied zwischen den beiden harten Jets verlangt. 



In Abbildung 5.18| ist deshalb die Verteilung des Rapiditätsunterschiedes auf den 



beiden Monte-Carlo- Niveaus und bei den Daten aufgetragen. Wir erkennen eine gute 
Beschreibung der Daten durch den Detektorlevel des Monte-Carlos. Die Partonjets 
sind jedoch im Durchschnitt etwas dichter beieinander und die Verteilung ist leicht 
verschoben. 



In Tabelle |5.21| sind die Migrationen für einen Mindestabstand von Ar] > 2 
eingetragen. Im Vergleich zu Tabelle |5.3| erkennen wir eine Erniedrigung des Korrekturfaktors 
und eine Verschlechterung der Reinheit. Auffällig ist weiterhin die Änderung der 
Abhängigkeit. Die Raten schwanken zwar auf beiden Niveaus und steigen leicht an, 
jedoch ist die Steigung minimal und im achten Bin ist eine starke Reduktion zu sehen. 
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Bin 3 


Parton 


Bin 3 


Parton 


Bin 3 


Parton 


Schnitt 1 


1+1 


2+1 


Schnitt 2 


1+1 


2+1 


Schnitt 3 


1+1 


2+1 


Clus- 


1+1 


9507 


147 


Clus- 


1+1 


9233 


235 


Clus- 


1+1 


9043 


294 


ter 


2+1 


157 


73 


ter 


2+1 


251 


165 


ter 


2+1 


314 


233 



Bin 6 


Parton 


Bin 6 


Parton 


Bin 6 


Parton 


Schnitt 1 


1 + 1 


2+1 


Schnitt 2 


1 + 1 


2+1 


Schnitt 3 


1 + 1 


2+1 


Clus- 


1+1 


6674 


261 


Clus- 


1+1 


6352 


342 


Clus- 


1+1 


6125 


404 


ter 


2+1 


184 


123 


ter 


2+1 


261 


287 


ter 


2+1 


288 


425 



Bin 8 


Parton 


Bin 8 


Parton 


Bin 8 


Parton 


Schnitt 1 


1+1 


2+1 


Schnitt 2 


1 + 1 


2+1 


Schnitt 3 


1 + 1 


2+1 


Clus- 


1+1 


1540 


149 


Clus- 


1+1 


1412 


171 


Clus- 


1+1 


1361 


153 


ter 


2+1 


56 


81 


ter 


2+1 


77 


166 


ter 


2+1 


85 


227 



Bin 


3 


6 


8 


Schnitt 


1 


2 


3 


1 


2 


3 


1 


2 


3 


Parton 


2.2 


4.0 


5.3 


5.3 


8.7 


11.4 


12.6 


18.5 


20.8 


in % 


±0.2 


±0.2 


±0.3 


±0.3 


±0.4 


±0.4 


±0.8 


±1.0 


±1.0 


Detektor 


2.3 


4.2 


5.5 


4.2 


7.6 


9.8 


7.5 


13.3 


17.1 


in% 


±0.2 


±0.3 


±0.3 


±0.3 


±0.4 


±0.4 


±0.7 


±0.8 


±0.9 


Korrek- 


0.96 


0.96 


0.96 


1.25 


1.15 


1.16 


1.68 


1.39 


1.22 


turfa. 


±0.13 


±0.10 


±0.09 


±0.14 


±0.10 


±0.08 


±0.25 


±0.16 


±0.12 


Reinheit 


31.7 


39.7 


42.6 


40.1 


52.4 


59.6 


59.1 


68.3 


72.8 


in % 





















Tabelle 5.20: Migrationen und Korrekturfaktoren mit verschiedenen Av^ 
Schnitten Exemplarisch werden drei Bins mit drei verschiedenen Av^ Schnitten 
gezeigt. Schnitt 1 ist hier < 5 GeV, Schnitt 2 A\/§ < 10 GeV und Schnitt 3 
Av^ < 15 GeV. 
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Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


1 


1+1 


2+1 


2 


1+1 


2+1 


3 


1+1 


2+1 


Clus- 


1+1 


12705 


260 


Clus- 


1+1 


9593 


220 


Clus- 


1+1 


9208 


231 


ter 


2+1 


420 


176 


ter 


2+1 


327 


129 


ter 


2+1 


323 


122 



Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


4 


1+1 


2+1 


5 


1+1 


2+1 


6 


1+1 


2+1 


Clus- 


1+1 


8 100 


276 


Clus- 


1+1 


5755 


169 


Clus- 


1+1 


6555 


276 


ter 


2+1 


294 


133 


ter 


2+1 


229 


125 


ter 


2+1 


230 


181 



Bin 


Parton 


Bin 


Parton 


7 


1+1 


2+1 






1+1 


2+1 


Clus- 


1+1 


4083 


152 


Clus- 


1+1 


1731 


35 


ter 


2+1 


154 


115 


ter 


2+1 


32 


28 



Bin 


1 


2 


3 


4 


5 


6 


7 


8 


Rate Parton 


3.2 


3.4 


3.6 


4.5 


4.7 


6.3 


5.9 


3.5 


in % 


±0.2 


±0.2 


±0.2 


±0.3 


±0.3 


±0.3 


±0.4 


±0.5 


Rate Detektor 


4.4 


4.4 


4.5 


4.7 


5.6 


5.7 


6.0 


3.3 


in 'X 


±0.2 


±0.3 


±0.3 


±0.3 


±0.3 


±0.3 


±0.1 


±0.5 


Korrekt ur- 


0.73 


0.77 


0.79 


0.96 


0.83 


1.11 


0.99 


1.05 


faktor 


±0.07 


±0.08 


±0.08 


±0.10 


±0.10 


±0.11 


±0.12 


±0.27 


Reinheit in % 


29.5 


28.3 


27.4 


31.1 


35.3 


44.0 


42.8 


46.7 



Tabelle 5.21: Migrationen und Korrekturfaktoren bei Lepto mit 
Standardschnitten und zusätzlichem Arj > 2 Schnitt, in Bins 
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Partor 

Detektor 

Daten 



0.5 



1 .5 



2.5 



3.5 



Abbildung 5.18: Vergleich der Ar] Verteilungen Die Verteilungen sind auf gleiche 
Flächen normiert. Dargestellt sind Parton- (durchgezogen) und Detektorniveau 
(gestrichelt) und die Daten (gepunktet). 

Ein Schnitt in dieser Variable ist für Jetanalysen ungeeignet. 



5.8.6 Quarkwinkel 9. 



q 



Wie wir schon in Kapitel |5.8.1| an Abbildung ß.l2^ und b gesehen haben, treten 



große Jetwinkel hauptsächlich bei 1+1 Ereignissen auf, wohingegen nahezu alle 2+1 
Ereignisse einen Jet unter kleinen Winkeln haben. Durch einen Schnitt ließe sich 
somit das Verhältnis von 2+1 zu 1+1 Jet Ereignissen und damit die Rate verändern. 
Besser als der Jetwinkel eignet sich hierfür der Winkel des nach dem harten Stoß 
auslaufenden Quarks im Quark Parton Modell. Dieser Winkel läßt sich wegen der 
relativ einfachen Reaktionskinematik des Born Graphen aus der Messung des gestreuten 
Elektrons durch die folgende Gleichung herleiten : 

Die Werte für und i/b werden dabei wie in Kapitel |5.2.1| erläutert berechnet. 

Abbildung [5.19| zeigt die Verteilung des Quarkwinkels für alle und für die 2+1 
Ereignisse in beiden Bereiche getrennt. Die Verteilungen zeigen eine unbefriedigende 
Beschreibung der Daten. Nur im Bereich < 100 GeV ist die Beschreibung gut. 
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Das gleiche gilt für den Unterschied zwischen den Monte-Carlo-Niveaus. Ein Schnitt 
auf einen maximalen Quarkwinkel verwirft im hohen Bereich wesentlich mehr 1+1 
Ereignisse als solche mit 2+1 Jets. Wir erwarten daher eine Ratensteigerung. Im 

< 100 GeV Bereich gilt die gleiche Aussage. Die Verteilungen sind jedoch steiler, 
so daß ein Schnitt hier die Statistik stark reduziert. 

Ein Schnitt auf einen maximalen Quarkwinkel verwirft mehr 1+1 Ereignisse als 
solche mit 2+1 Jets. Wir erwarten somit eine Ratensteigerung. Der Effekt ist bei 
einer Grenze von 100 Grad am deutlichsten zu sehen. Die Migrationen, die sich aus 
der Anwendung dieses Schnittes ergeben, sond in Tabelle |5.22| gezeigt. Ein Vergleich 
mit Tabelle ^]3| zeigt die Bestätigung des oben gesagten. Die Raten sind auf beiden 
Niveaus gestiegen, wobei es in den verschiedenen Bins jedoch starke Unterschiede 
gibt. Die Reinheit ist in etwa konstant geblieben und der Korrekturfaktor hat sich 
nur in den drei unteren Bins verändert. In diesen Bins macht sich die Reduktion der 
Statistik stark bemerkbar. Der relative Fehler ist von ca. 5% ohne öquark Schnitt auf 
ca. 20% mit Schnitt angestiegen. 

Unter Berücksichtigung dieser Ergebnisse ist ein Schnitt im Quarkwinkel nur in 
den hohen Bins sinnvoll. Er verbessert die Migrationen jedoch nicht. 



5.9 Alternative Ratenberechnung 

Im letzten Abschnitt dieses Kapitels untersuchen wir einen alternativen Ansatz die 
Jetrate zu bestimmen. Dieser verzichtet auf einen Jetalgorithmus und macht sich die 
Tatsache zu nutze, daß die Energieverteilung im Breitsystem je nach Anzahl der Jets 
unterschiedlich ist[[Hoe95|. Hierzu machen wir uns klar, wie sich 2+1 Ereignisse im 



Breitsystem verhalten. Dabei tritt das Problem auf, daß s ungleich Null ist (siehe 
Kapitel p.8.4D und somit nach Gleichung ß.ll| die Bjorken Skalenvariable xb nicht 



dem Partonimpulsbruchteil Xi/p entspricht. Deshalb wird das einlaufende Parton 
nach der Transformation nicht mehr durch das Austauschteilchen wie an einer harten 
Wand reflektiert (siehe Kapitel [4.1.3| und Abbildung ^]^) . Es verbleibt statt dessen 



ein Restimpuls. Im Spezialfall Xi/p = 2xb ist der Impuls des Photons entgegengesetzt 
gleich groß wie der Impuls des Partons. Die Jets laufen in entgegengesetzte Richtungen 
aus. Ist Xi/p noch größer, werden sogar beide Jets in die Photon Hemisphäre gelenkt. 
Dies wird in Abbildung p.20| verdeutlicht. 



Legen wir nun zwei Bereiche im Breitsystem fest. Die Grenze wird durch einen 
Winkel 9^ relativ zur einlaufenden Partonrichtung festgelegt. Wir benutzen im 
folgenden 90 Grad zur Trennung der sogenannten Current Hemisphäre, d.i. die Seite 
mit dem einlaufenden Parton, und der Target Hemisphäre, d.i. die Seite mit einlaufendem 
Photon. Wir erwarten in der Current Hemisphäre die gesamte Schwerpunktenergie 
|(5 im Falle eines 1+1 Ereignisses und bei 2+1 Jet Ereignissen mit Xß ~ a^j/p- Bei 
Xi/p » 2a;ß erwarten wir keine Energie in der Current Hemisphäre. Abbildung 



5.21 zeigt die Verteilung auf den beiden Monte-Carlo-Niveaus und bei den Daten. 



Auf Partonniveau entspricht die Verteilung den Erwartungen. Das Detektorniveau 
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Abbildung 5.19: Vergleich der ÖQuark Verteilungen a) zeigt alle Ereignisse mit 
> 100 GeV, b) die 2+1 Ereignisse, c) und d) sind die analogen Graphen für die 
Bins mit < 100 GeV Dargestellt sind Parton- (durchgezogen) und Hadronniveau 
(gestrichelt) und die Daten (gepunktet). Die Verteilungen sind jeweils auf gleiche 
Flächen normiert. 
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Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


1 


1+1 


2+1 


2 


1+1 


2+1 


3 


1 + 1 


2+1 


Clus- 


1+1 


198 


5 


Clus- 


1+1 


668 


17 


Clus- 


1+1 


1183 


27 


ter 


2+1 


14 


10 


ter 


2+1 


29 


20 


ter 


2+1 


55 


44 



Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


4 


1+1 


2+1 


5 


1+1 


2+1 


6 


1+1 


2+1 


Clus- 


1+1 


1726 


89 


Clus- 


1+1 


2165 


121 


Clus- 


1+1 


1526 


187 


ter 


2+1 


83 


100 


ter 


2+1 


133 


148 


ter 


2+1 


136 


261 



Bin 


Parton 


Bin 


Parton 


7 


1+1 


2+1 


8 


1 + 1 


2+1 


Clus- 


1+1 


2187 


253 


Clus- 


1+1 


1069 


137 


ter 


2+1 


178 


490 


ter 


2+1 


92 


260 



Bin 


1 


2 


3 


4 


5 


6 


7 


8 


Rate Parton 


6.6 


5.0 


5.4 


9.5 


10.5 


21.2 


23.9 


25.5 


in% 


±1.7 


±0.9 


±0.7 


±0.7 


±0.7 


±0.9 


±0.8 


±1.2 


Rate Detektor 


10.6 


6.7 


7.6 


9.2 


10.9 


18.8 


21.5 


22.6 


in % 


±2.1 


±1.0 


±0.8 


±0.7 


±0.7 


±0.9 


±0.8 


±1.1 


Korrekt ur- 


0.62 


0.76 


0.72 


1.03 


0.96 


1.13 


1.11 


1.13 


faktor 


±0.28 


±0.23 


±0.16 


±0.15 


±0.11 


±0.10 


±0.08 


±0.11 


Reinheit in % 


41.7 


40.8 


44.4 


54.6 


52.7 


65.7 


73.4 


73.9 



Tabelle 5.22: Migrationen und Korrekturfaktoren bei Lepto mit 
StandEirdschnitten und zusätzlichem ÖQuark < 100° Schnitt, in Q'^ Bins 




Abbildung 5.20: 2+1 Jetereignisse im Breitsystem a) bei Xi/p — 2xb, b) bei 
Xi/p » 2xb- 
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Abbildung 5.21: Vergleich der Energieverteilung in der Current Hemisphäre 
im Breitsystem Dargestellt sind Parton- (durchgezogen) und Detektorniveau 
(gestrichelt) des Lepto-Monte-Carlos und die Daten (gepunktet). Die Verteilungen 
sind auf gleiche Flächen normiert. 



beschreibt die Daten einigermaßen, die Verteilungen sind aber sehr breit. Wir 
benutzen nun einen Abschneideparameter e = 0.25 und definieren Ereignisse mit 
-E'cur < f Q als 2+1 Jet Ereignisse und die mit größerer Current Energie als 1+1 Jet 
Ereignisse. 

Die Migrationen in Tabelle |5.23| zeigen ein merkwürdiges Verhalten. Die 
Raten fallen mit steigendem Q^. Ob dies durch die Theorie beschrieben wird, ist 
unklar, da die vorhandenen Theorieprogramme die Partonjets erzeugen, nicht jedoch 
die Cluster. Ein Vergleich in nächst-zu-führender Ordnung ist somit nicht möglich. 
Auch die Reinheit der Jetklassen ist in den unteren vier Bins besser als beim 
JADE Algorithmus, während die oberen Bins eine schlechtere Reinheit zeigen. Der 
Korrekturfaktor zeigt ein uneinheitliches Verhalten. 

Ein Vergleich zwischen der Jetklassifizierung über die Current Energie mit dem 
JADE Algorithmus ist in Tabelle |5.24| dargestellt. Wie wir schon aufgrund des völlig 
unterschiedlichen Verhaltens erwarten konnten, ist die Korrelation sehr schlecht. 

Diese Jetanzahlbestimmung hat den Vorteil, daß sie die Klassifizierung nicht über 
einen Algorithmus durchführt, sondern durch Berechnung nur eines Wertes. Der 
Nachteil dieses Verfahrens, neben den oben angesprochenden, ist, daß die Jets selber 
nicht bestimmt werden. Es sind somit keine weiteren Schnitte auf z.B. Jetwinkel 
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1474 


2221 


ter 


2+1 


1187 


1479 


ter 


2+1 


1053 


1248 


Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


4 


1+1 


2+1 




5 


1+1 


2+1 


6 


1+1 


2+1 


Clus- 


1+1 


6691 


535 


Clus- 


1+1 


5117 


252 


Clus- 


1+1 


6648 


165 


ter 


2+1 


932 


945 


ter 


2+1 


499 


410 


ter 


2+1 


165 


264 


Bin 


Parton 


Bin 


Parton 










7 


1+1 


2+1 






1 + 1 


2+1 










Clus- 


1+1 


4308 


38 


Clus- 


1+1 


1803 


5 










ter 


2+1 


77 


81 


ter 


2+1 


11 


7 











Bin 


1 


2 


3 


4 


5 


6 


7 


8 


Rate Parton 


25.9 


23.1 


20.1 


16.3 


10.5 


5.9 


2.6 


0.7 


in% 


±0.4 


±0.5 


±0.5 


±0.4 


±0.4 


±0.3 


±0.3 


±0.2 


Rate Detektor 


27.2 


26.0 


23.3 


20.6 


14.5 


5.9 


3.5 


1.0 


in % 


±0.4 


±0.5 


±0.5 


±0.5 


±0.5 


±0.3 


±0.3 


±0.3 


Korrektur- 


0.95 


0.89 


0.86 


0.79 


0.73 


1.00 


0.75 


0.67 


faktor 


±0.03 


±0.04 


±0.04 


±0.04 


±0.05 


±0.10 


±0.13 


±0.35 


Reinheit in % 


60.1 


55.5 


54.2 


50.3 


45.1 


61.5 


51.3 


38.9 



Tabelle 5.23: Migrationen und Korrekturfaktoren bei Lepto in Bins 
mit Standcirdschnitten bei Anzahlbestimmung durch die Energie in der 
Current Hemisphäre 



Parton- 


alternativ 


Detektor- 


alternativ 


niveau 


1+1 


2+1 


niveau 


1+1 


2+1 


JADE 


1+1 


56980 


9442 


JADE 


1+1 


49991 


10659 




2+1 


3766 


2969 




2+1 


3379 


2529 



Tabelle 5.24: Vergleich vom JADE Algorithmus mit der Ratenbestimmung 
durch die Energie in der Current Hemisphäre 



KAPITEL 5. ERGEBNISSE DER UNTERSUCHUNGEN 



Kapitel 6 

Zusammenfassung und Ausblick 



6.1 Zusammenfassung 

Aus den im vorigen Kapitel gezeigten Untersuchungen werden die unten genannten 
Schnitte gefolgert : 



die Standardschnitte für die Kinematik (Kapitel |5. 2. 2| ), die Rekonstruktionsqualität 



(Kapitel |5]23) und die Elektronfinder (Kapitel ^2^. 



und der Jetwinkelschnitt 15° < öjct < 145° 



Mit diesen Schnitten ergeben sich die in Tabelle ß.l| aufgetragenen Migrationen. 



Im Vergleich zur Tabelle |5.3| zeigen sich hier die kleineren Raten und der erhöhte 
Korrekturfaktor. Der ungewichtete Mittelwert des Korrekturfaktors ist 1.683 ± 
0.191. Im Rahmen ihrer Fehler sind alle Werte mit dem mittleren Korrekturfaktor 
verträglich. Auch die Schwankungen in der Reinheit sind kleiner. Die Reinheit selber 
hat sich durch den Schnitt erhöht. 

Vergleichen wir nun die in Kapitel ^ untersuchten Verteilungen noch einmal 



unter Berücksichtigung des Jetwinkelschnittes in den Abbildungen ^]T] und |6]^. Bei 
der öjet Verteilung haben wir im Vorwärtsbereich nicht den gesamten ungenügend 
beschriebenen Bereich durch unseren Schnitt beseitigt. Die Gründe für die Wahl 
dieser Grenzen sind schon in Kapitel ß.8.1| besprochen worden. Die Verteilung für die 
2+1 Raten stimmt nun jedoch auf allen drei Niveau hervorragend überein. (vergleiche 
Abbildung |5.12|) . Für die Verteilung des transversalen Impulses im hadronischen 
Schwerpunktsystem gilt die gleiche Aussage (vergleiche Abbildung |5.16|) . Die Gesamtverteilun^ 



stimmt aus Gründen der Energiekalibration nur für den Vergleich zwischen Monte- 
Carlo-Detektorniveau und Daten. Auch hier hat sich die Verteilung der 2+1 Jet 
Ereignisse wesentlich verbessert. Die Verschiebung zwischen Parton- und Detektorniveau 
wurde also durch die schlecht beschriebene Jetwinkelverteilung ausgelöst. 

Die Quarkwinkelverteilungen waren vor dem Schnitt unbefriedigend beschrieben 
(siehe Abbildung |5.19|) . Dies hat sich auch mit der Einschränkung des Jetwinkelbereiches 
nicht geändert (Abbildung |0| a und b). Das Aussehen der Zp Verteilung hat sich 
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KAPITEL 6. ZUSAMMENFASSUNG UND AUSBLICK 



Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


1 


1+1 


2+1 


2 


1+1 


2+1 


3 


1+1 


2+1 


Clus- 


1+1 


6716 


128 


Clus- 


1+1 


6343 


74 


Clus- 


1+1 


6584 


114 


ter 


2+1 


61 


81 


ter 


2+1 


30 


52 


ter 


2+1 


34 


52 



Bin 


Parton 


Bin 


Parton 


Bin 


Parton 


4 


1+1 


2+1 


5 


1+1 


2+1 


6 


1+1 


2+1 


Clus- 


1+1 


6617 


121 


Clus- 


1+1 


5351 


91 


Clus- 


1+1 


5113 


210 


ter 


2+1 


40 


56 


ter 


2+1 


26 


43 


ter 


2+1 


66 


176 



Bin 


Parton 


Bin 


Parton 


7 


1 + 1 


2+1 




i 


1+1 


2+1 


Clus- 


1+1 


3954 


172 


Clus- 


1+1 


1508 


62 


ter 


2+1 


57 


137 


ter 


2+1 


21 


74 



Bin 


1 


2 


3 


4 


5 


6 


7 


8 


Rate Parton 


3.0 


1.9 


2.4 


2.6 


2.4 


7.1 


7.2 


8.2 


in% 


±0.3 


±0.2 


±0.2 


±0.2 


±0.3 


±0.4 


±0.4 


±0.7 


Rate Detektor 


2.0 


1.3 


1.3 


1.4 


1.3 


4.1 


4.5 


5.7 


in 'X 


±0.2 


±0.2 


±0.2 


±0.2 


±0.2 


±0.3 


±0.1 


±0.6 


Korrektur- 


1.47 


1.54 


1.93 


1.84 


1.94 


1.72 


1.59 


1.43 


faktor 


±0.23 


±0.31 


±0.36 


±0.33 


±0.40 


±0.19 


±0.20 


±0.27 


Reinheit in % 


57.0 


63.4 


60.5 


58.3 


62.3 


72.7 


70.6 


77.9 



Tabelle 6.1: Migrationen und Korrekturfaktoren bei Lepto mit 
Standardschnitten und zusätzlichem 15° < Oj^t < 145° Schnitt, in Bins 
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5 10 15 5 10 15 

Pi'/GeV Pi'/GeV 
c) d) 



Abbildung 6.1: Vergleich der wichtigsten Verteilungen mit abschließenden 

Schnitten a) Jetwinkelverteilung für alle Ereignisse, b) inklusive Jetwinkelverteilung 
für 2+1 Jet Ereignisse, c) und d) transversale Impulsverteilung für das hadronische 
Schwerpunktsystem für alle bzw. nur für 2+1 Jet Ereignisse. Dargestellt sind Parton- 
( durchgezogen), Detektorniveau (gestrichelt) und Daten (gepunktet) 
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KAPITEL 6. ZUSAMMENFASSUNG UND AUSBLICK 





1 


2 


3 


4 


5 


6 


7 


8 


Datenrate -R2+1, Daten in % 


JL 


1.2 
±0.6 


1.3 
±0.6 


1.2 
±0.6 


1.5 
±0.6 


1.7 
±0.7 


4.7 
±1.0 


5.5 
±1.1 


5.7 
±1.5 


JC 


0.8 
±0.6 


0.8 
±0.6 


0.6 

±0.6 


1.0 

±0.6 


1.2 
±0.7 


3.2 
±0.9 


4.2 

±1.1 


5.3 
±1.5 


JB 


0.8 
±0.6 


1.0 
±0.6 


1.0 
±0.6 


1.3 
±0.6 


1.8 
±0.7 


5.2 
±1.0 


7.0 
±1.2 


9.5 
±1.8 


3L 


0.5 
±0.6 


0.5 
±0.6 


0.3 
±0.6 


0.6 
±0.6 


0.7 
±0.6 


2.1 
±0.8 


2.9 
±1.0 


3.5 
±1.3 


CL 


17.0 
±2.8 


14.8 
±2.7 


9.2 
±2.0 


8.8 
±1.6 


4.2 
±1.1 


4.6 
±1.0 


5.7 
±1.1 


8.0 
±1.7 


EL 


0.5 
±0.6 


0.4 
±0.6 


0.5 
±0.6 


0.8 
±0.6 


0.7 
±0.6 


2.1 
±0.8 


2.3 
±0.9 


4.9 
±1.5 



Tabelle 6.2: Zusammenfassung des Vergleichs verschiedener Jetalgorithmen 

mit 15° < öjet < 145° Schnitt 

geändert. Im untersten Bin sind nun wesentlich weniger Einträge. In diesem Bin 
ist die Ubereinstimmung zwischen Parton- und Detektorniveau nicht mehr gegeben. 
Dies kann in Anbetracht der Stärke der Steigung jedoch ein reiner Migrationseffekt 



sein (vergleiche Abbildung ß.l5|) . Die Verteilung des Rapiditätsunterschiedes hat sich 
am stärksten geändert. Dies ist auch verständlich, da durch den Jetwinkelschnitt die 
Jets unter extremen Winkeln wegfallen. Daher finden wir nur noch solche Ereignisse, 
deren Jets wesentlich dichter zusammen liegen. 

Betrachten wir nun noch einmal die Abhängigkeit des Ergebnisses vom Jet- 
algorithmus. Der Jetwinkelschnitt wurde dabei immer im Laborsystem angewandt, 
d.h. die Partikel wurden erst in das entsprechende System transformiert, danach 
der Jetalgorithmus angewandt und danach wurden die Jets zurücktransformiert. 
Anschließend wurde der Jetwinkelschnitt durchgeführt. Das Ergebnis ist in den 
Abbildungen ^]3| bis |6.7| und den Tabellen |6.2| bis ^]4| gezeigt. Ein Vergleich mit 
den entsprechenden Abbildungen und Tabellen im Kapitel |5.4| verdeutlicht, daß 
wir durch den Schnitt im Jetwinkel wesentlich kleinere Raten und einen erhöhten 
Korrekturfaktor finden. Die Werte bei den verschiedenen Jetalgorithmen zeigen 
nun jedoch kleinere Abweichungen. Die Abhängigkeit vom Jetalgorithmus und vom 
Lorentzsystem ist deutlich kleiner geworden. Sicherlich kann noch nicht von einer 
Übereinstimmung der verschiedenen Algorithmen gesprochen werden. 
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0.2 r 
0.175 
0.15 

0.125 

0.1 r 
0.075 ^ 

0.05 
0.02b ^ 






0.18 
0.1 6 ^ 
0.14 
0.12 ^ 

0.1 
0.08 ^ 
0.06 ^ 
0.04 
0.02 ^ 






50 



100 



a) 



150 




0.2 0.3 0.4 0.5 

c) 



0.225 
0.2 
0.1 75 
0.15 
0.125 
0.1 
0.07b 
0.05 
0.025 




150 



2 3 4 5 

Ar] 

d) 



Abbildung 6.2: Vergleich der wichtigsten Verteilungen mit abschließenden 
Schnitten a) und b) Quarkwinkelverteilung für alle bzw. 2+1 Jet Ereignisse, c) Zp 
Verteilung und d) Rapiditätsdiff'erenz der beiden harten Jets. Dargestellt sind Parton- 
(durchgezogen), Detektorniveau (gestrichelt) und Daten (gepunktet) 
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KAPITEL 6. ZUSAMMENFASSUNG UND AUSBLICK 



JL • • 
JC • • 
JB • • 

3L • • • • 

EL • • 

I I I I I I I I I ^1 I I I I I I I I I ^1 I I I I I I I ^1 I I I I I I I I L 

0.1 0.1 1 2 0.1 

a) Daten b) MC Detektor c) Korrektjrfoktor d) korrigierte Rote 



Abbildung 6.3: Vergleich verschiedener Jetalgorithmen, Bin 1-4. Dargestellt 
sind auf der Ordinate die Jetalgorithmen JADE im Laborsystem (JL), JADE im 
hadronischen Schwerpunktsystem (JC), JADE im Breitsystem (JB), JADE mit 
Vcnt = 0.03 (3L), Cone Algorithmus (CL) und der E Algorithmus (EL) (siehe auch 
Tabelle In den vier Craphen sind die Raten für die Daten (a) und das Monte- 
Carlo-Detektorniveau (b). Der Korrekturfaktor ist in c gezeigt und die sich ergebende 
korrigierte Rate in d. 
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0.1 0.1 1 2 0.1 

a) Daten b) MC Detektor c) Korrekturfaktor d) korrigierte Rate 



Abbildung 6.4: Vergleich verschiedener Jetalgorithmen, Bin 5. Für nähere 



Erläuterungen siehe Abbildung |g.5| und Tabelle 5.1 



3L ' • • 

EL * • 

0.1 0.1 1 2 0.1 

o) Daten b) MC Detektor c) Korrekturfaktor d) korrigierte Rote 

Abbildung 6.5: Vergleich verschiedener Jetalgorithmen, Bin 6. Für nähere 
Erläuterungen siehe Abbildung \6.^ und Tabelle |5. T\ 
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KAPITEL 6. ZUSAMMENFASSUNG UND AUSBLICK 





0.1 0.1 1 2 0.1 

a) Daten b) MC Detektor c) Korrekturfaktor d) korrigierte Rate 



Abbildung 6.6: Vergleich verschiedener Jetalgorithmen, Bin 7. Für nähere 



Erläuterungen siehe Abbildung |g.5| und Tabelle 5.1 





0.1 0.1 1 2 0.1 

o) Daten b) MC Detektor c) Korrekturfaktor d) korrigierte Rote 



Abbildung 6.7: Vergleich verschiedener Jetalgorithmen, Bin 8. Für nähere 
Erläuterungen siehe Abbildung \6.^ und Tabelle |5. T\ 
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1 


2 


3 


4 


5 


6 


7 


8 


Monte-Carlo-Detektorrate i?2+i, Cluster in % 


JL 


1.8 
±0.6 


1.1 
±0.6 


1.2 
±0.6 


1.3 
±0.6 


1.3 
±0.6 


3.9 
±0.7 


4.7 

±0.8 


5.6 
±1.1 


JC 


0.8 
±0.6 


0.7 
±0.6 


0.8 
±0.6 


0.9 
±0.6 


1.0 
±0.6 


2.6 
±0.7 


3.3 
±0.8 


4.4 
±1.0 


JB 


1.3 

±0.6 


0.8 

±0.6 


0.9 

±0.6 


1.2 

±0.6 


1.5 
±0.7 


4.6 
±0.8 


6.4 
±0.9 


8.4 
±1.2 


3L 


0.6 

±0.6 


0.4 
±0.6 


0.3 
±0.6 


0.5 
±0.6 


0.7 
±0.6 


1.7 
±0.7 


2.4 
±0.7 


3.1 
±0.9 


GL 


16.5 
±2.1 


11.3 
±2.2 


9.6 
±1.8 


11.2 
±1.8 


4.3 
±1.1 


4.9 
±0.8 


4.9 
±0.8 


8.9 
±1.2 


EL 


1.6 

±0.6 


1.2 

±0.6 


1.0 

±0.6 


1.1 
±0.6 


1.0 
±0.6 


2.7 
±0.7 


2.6 
±0.7 


5.4 
±1.1 


Montc-Carlo-Partonrate -R2+i,Parton in % 


JL 


2.9 
±0.7 


1.9 
±0.6 


2.5 
±0.7 


2.8 
±0.7 


2.8 
±0.7 


7.4 
±0.8 


7.5 
±0.9 


8.6 
±1.1 


JC 


1.5 
±0.6 


1.2 
±0.6 


1.3 
±0.6 


1.6 
±0.6 


2.1 
±0.7 


4.9 
±0.8 


6.1 
±0.8 


7.1 
±1.1 


JB 


5.1 
±0.9 


3.0 
±0.8 


3.6 

±0.8 


3.5 

±0.8 


3.7 

±0.8 


8.8 
±0.9 


9.5 
±0.9 


10.7 
±1.2 


3L 


1.1 
±0.6 


0.8 
±0.6 


0.9 
±0.6 


1.1 
±0.6 


1.1 
±0.6 


3.5 
±0.7 


4.3 
±0.8 


5.0 
±1.0 


CL 


24.2 
±1.8 


19.5 
±1.9 


16.4 
±1.7 


14.2 
±1.5 


4.4 
±0.9 


10.4 
±0.9 


10.7 
±0.9 


14.4 
±1.2 


EL 


3.1 
±0.7 


1.9 
±0.6 


2.4 
±0.7 


2.8 
±0.7 


2.8 
±0.7 


8.0 
±0.8 


8.3 
±0.9 


10.1 
±1.2 



Tabelle 6.3: Zusammenfassung des Vergleichs verschiedener Jetalgorith 

mit 15° < öiet < 145° Schnitt 
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KAPITEL 6. ZUSAMMENFASSUNG UND AUSBLICK 





1 


2 


3 


4 


5 


6 


7 


8 


Korrekturfaktor R 


JL 


1.63 
±0.22 


1.77 
±0.33 


2.06 
±0.36 


2.11 
±0.35 


2.17 
±0.41 


1.92 
±0.20 


1.60 
±0.19 


1.54 
±0.27 


JC 


1.94 
±0.39 


1.70 
±0.43 


1.63 
±0.38 


1.71 
±0.36 


2.06 
±0.47 


1.88 
±0.24 


1.85 
±0.26 


1.62 
±0.32 


JB 


3.80 

±0.66 


3.62 

±0.83 


4.10 

±0.86 


2.88 

±0.54 


2.45 

±0.45 


1.92 

±0.20 


1.48 

±0.15 


1.27 

±0.19 


3L 


1.82 
±0.43 


2.24 
±0.70 


2.82 
±0.88 


2.09 
±0.55 


1.63 
±0.45 


2.02 
±0.32 


1.83 
±0.30 


1.62 
±0.38 


GL 


1.47 
±0.22 


1.73 
±0.39 


1.70 
±0.37 


1.27 
±0.23 


1.02 
±0.26 


2.13 
±0.21 


2.19 
±0.24 


1.62 
±0.22 


EL 


1.98 
±0.28 


1.60 
±0.29 


2.34 
±0.43 


2.58 
±0.45 


2.86 
±0.59 


2.94 

±0.35 


3.22 
±0.47 


1.87 
±0.33 


korrigierte Rate i?2+i,korr. in % 


JL 


2.0 
±1.0 


2.3 
±1.2 


2.5 
±1.2 


3.1 
±1.3 


3.6 
±1.5 


9.0 
±2.3 


8.9 
±2.5 


8.8 
±3.5 


JC 


1.5 
±1.0 


1.3 
±1.0 


1.0 
±0.9 


1.7 
±1.1 


2.5 
±1.4 


6.0 
±2.0 


7.8 
±2.6 


8.6 
±3.8 


JB 


3.0 

±1.4 


3.6 
±1.8 


4.1 
±1.9 


3.6 
±1.5 


4.3 
±1.7 


10.0 
±2.5 


10.4 
±2.6 


12.1 

±3.8 


3L 


0.9 
±0.8 


1.1 
±1.0 


0.8 
±0.9 


1.2 
±1.0 


1.2 
±1.0 


4.3 
±1.8 


5.3 
±2.2 


5.7 
±3.1 


CL 


25.0 
±7.6 


25.6 
±9.9 


15.7 
±6.4 


11.2 
±3.9 


4.3 
±2.1 


9.8 
±2.6 


12.6 
±3.2 


12.9 
±4.1 


EL 


0.9 
±0.8 


0.6 
±0.7 


1.2 
±0.9 


1.9 
±1.1 


2.1 
±1.2 


6.3 
±2.2 


7.4 
±2.9 


9.1 
±3.9 



Tabelle 6.4: Zusammenfassung des Vergleichs verschiedener Jetalgorithmen 

mit 15° < öiet < 145° Schnitt 
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6.2 Ausblick 

Zusammenfassend konnte deutlich gemacht werden, daß die Betrachtung von Jet- 
Parton-Korrelationen einem Problem in einem multi-dimensionalen Phasenraum entspricht. 
Es war nicht möglich einen Schnitt in einer einzelnen kinematischen Größe zu finden, 
der eine Beseitigung der Migrationen bewirkt, d.h. der eine Verkleinerung des Korrekturfaktors 
und gleichzeitig eine Verbesserung der Reinheit der 2+1 Jetereignisklasse erreicht. 
Immerhin konnte die Reduktion der Abhängigkeit von Rekombinationsschema, Jetalgorithmus 
und Bezugssystem durch Einführung eines Jetwinkelschnittes gezeigt werden. 

Weitere Arbeit ist nötig um neue Verteilungen zu überprüfen. Zu den dringendsten 
Problemen gehören sicherlich die Erweiterung der Untersuchungen auf das Theorieniveau. 
Diese wird besonders interessant, wenn es gilt, verschiedene Programme zu vergleichen. 
Neben dem von uns benutzten Programm PRO JET | |Gra93 | sind zur Zeit zwei weitere 



Programme in der Test- (MEPJET, |[MZ95 |) bzw. Entwicklungsphase. Diese neuen 



Programme werden eine höhere Flexibilität in Bezug auf die Auswahl des Jetalgorithmus 
haben. 

Daneben sind jedoch auch neue Möglichkeiten entstanden durch den Austausch 
des rückwärtigen Kalorimeters BEMC durch das SPACAL |[BW95|| , welches nun auch 



einen hadronischen Teil besitzt. Sobald die Energiekalibration dort verstanden ist, 
kann der obere Winkelschnitt vielleicht fallen gelassen oder zumindest verändert 



werden. Wie wir aber an der Jetwinkelverteilung in Abbildung |5.12| gesehen haben, 
sind dort nur wenige Jets von 2+1 Jetereignissen. Außerdem hat sich die Luminosität 
bei der Datennahme 1995 wieder erhöht, so daß wir mehr Ereignisse und damit 
geringere statistische Fehler erwarten dürfen. 

Das immer besser werdende Verständnis von Monte-Carlo-Simulationen und die 
Untersuchung der Hadronisierungseffekte werden sicher auch weiterhin neue Aspekte 
zum Verständnis der Jet-Parton-Korrelation beitragen. Wünschenswert wäre jedoch 
ein Monte-Carlo-Ereignisgenerator, der auf den Partonschaueransatz zur Beschreibung 
der nächst-zu-führenden Ordnung Korrekturen verzichtet und sie statt dessen bereits 
in den Matrixelementen beinhaltet. Dem stehen jedoch nicht nur software-technische 
Probleme entgegen. 

Die in Kapitel |4.3| erwähnte Messung der Gluondichte des Protons durch eine 
Jetanalyse stellt durch die zusätzliche Xß Einteilung neue Anforderungen an das 
Verständnis der Jet-Parton-Korrelationen. Hier ist eine Untersuchung der Größen 
auf eine Xß Abhängigkeit von Nöten. 

Wie wir aus dieser sicher unvollständigen Aufzählung erkennen können, bleiben 
noch viele Fragen offen, die die Neugier des Physikers anregen. 
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